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CHAPTER 1 
Object of the Study 

In recent years, the study of the attitude dynamics of a space- 
craft considered as a partly rigid, partly elastic or articulated body 

F 1— 1 1 

has become of increasing importance*’ . At first, such work did not 
present such a degree of urgency, as many investigations concentrated 
on rotational and librational dynamics of essentially rigid spacecraft, 

F 1-21 

as is apparent from the reviews of D.B. De Bra*" and R.E. Rober- 

r 1 _ 3 1^1 

son^ * * Any elastic body effects are conspicuously absent of 

V,V. Beletskii^s classic book on the "Motion of an artificial satellite 
about its center of mass" who writes at the outset that "the discussion 
is confined to problems which fall within the scope of the dynamics 
of rigid bodies". 

Satellites became increasingly "elastic", as booms were extended 
tens and hundreds of meters from the central body^^ ^ ^ or as 

r X“ 8 1 

large ^olar panels or manned toroiddl space stations are considered^ ^ 

Three methods are most commonly used in the study of the dynamics of the 

elastic spacecraft: discretization by modeling the continuous system by 

[1-9] 

finite elements; modal representation; and the Likin ’s’’ method of 

hybrid coordinates. 

The present work uses the modal approach. It is a study of the 
relevant equations and parameters in the dynamical analysis of the 
attitudes motion of a spin-stabilized spacecraft having flexible appen- 
dages. It is principally aimed at developing working tools, such as 
stability diagrams, tables or simulation analyses by means of computer 
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programs. These programs are of low time-consumption, and their use 
is quite easy to learn. As such, it is hoped that they v/ill prove 
valuable to the engineer engaged in the design of spin-stabilized 
elastic spacecraft. 
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CHAPTER 2 

A Study of Modal Shapes and Eip;en£requencies of 
Flexible Appendages on a Spin-Stabilized Satellite 


2,1 Introduction 

In order to study the dynamics of the spin-stabilized satellite with 
flexible appendages, by the methods of generalized dynamics, the con- 
tinuum of the elastic parts should be represented by generalized coor- 
dinates q. (1 = 1,2,...). The q. are ^fun- tions of time describing 
^i i 

'"k 

the amplitude of the non-dimensional displacements, , of boom k at 

abscissa ^ terms of modes 


Wk(0|= S q^(t) «>^(0 (2.1-1) 

w, will be (in the assumption of small displacements) along y for equa- 

1C 

torial displacements (E) and along z for meridional displacements (M) 
(See Fig. 2.1) . 


n, ^ are the geometric coordinates x,y,z non-dimensionalized 
by £, undeflected length of the boom. Co ~ non-dimensional 

radius of the central hub. 

The system of mode shapes, adopted here are the modes of the 

[ 2 - 1 ] 

rotating structure corresponding to the boom's Etkin number 

A = and non-dimensional radius Co = ^ * P is the (uniform) 

El s ^ 
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lineal density of the boom, in units of mass/length. . E is the boom's 

Young modulus, in units of force/unit area, I is the geometric moment 

4 

of inertia of the boom's cross section, in units of length , and w is 
the spin rate, in rad. /sec. Thus X is non-dimensional. Finally, 
is the radius of the central hub, at which distance the elastic boom 
is assumed to be cantilevered. As will be seen, these significantly 
depart in shape and frequency from those of the non-rotating structure 
corresponding to X = 0 and Ko =0. 

In the following, it is assumed that only antisymmetric motions 


are considered, or that the motion of the CM away from the origin is 

[ 2-21 

negligible. The latter amounts, as has been shown by F. Vigneron 
to assuming that the central mass is sufficiently large for terms 
of order 

P^(^) [ w dx]^ 

^ boom 

to be neglected in comparison with terms like 


P 


f 

dx 

-'boom 


Typically, for the ALOUETTE and ISIS satellites. Ref . [2-2] gives the 

pZ 

values; ^ = 0.005 to 0.01 (copper-beryllium booms) 


2.2 Equations of Motion; equatorial vibrations 
2.2.1 Basic equation 

We shall first consider motions in the "equatorial" plane of the 

satellite, i.e. (x,y) or (S,n)* These were the first type of vibrations 

[2-31 

considered by this author and J.E. Rakowski 
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Any section of boom located at p , of abscissa x, is in rotational 
equilibrium under the action of (Fig. 2.1). 

- bending moment from the left, which for pure flexure in the equa- 
torial plane, is 


dM 


el 


= 1 


( 2 . 2 - 1 ) 




in which w(x) is the assumed small displacement of the boom element 
in the y-direction, and 1^ is the unit vector along the z-direction. 

- the moment about p of inertia forces 

dF . = -p dx. r^ (2.2-2) 

in 1 Q 

imparted by the particles of the boom to the right of p, i.e. 
having abscissa between x and 
Therefore, 



In terms of their components, we have 



Thus 


K 

a 


-T ~r f \ 




/ 


' ^o+^ ' 

L 0 J 


w(y,) 

w(x,)+ w., 


I 

wYx) 

. 0 
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Also 

- £ 02 _ W(X{) - 

K-- 1 ()t ) t w A = Wj (Xo+X|) + wfx,) 

1 w(*,) + «x w(x/j-«j (x„+x,j 

i| ' 

/ l-w(x,)u)^ W(x,)^w^(x,^x,) 

(“x - (‘*^j+“z)(^o+x,) -2(0, w(X|) 

( w (X|)- 10.^ w(x,) ) - w(y,) + f(^^ + 0)^ Wj ) fx„ +x, 

(w<U^-W^)(Xo+Xi]+ (iO;(+ W^tojj wlXi] +2 (Ox w(X|) 
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Under the assumption of small displacements and transverse angular 
rates, terms of order w^, neglected, and r^ re- 

duces to 


- 0 J^W(x,) - w(x,j 

(w(x,)-wt W(x,)) + 4 (Xo^X.) 

K + X,j 


(2.2-3) 


Finally, along 1 


[(^ -;t ) aT 1 ^ 

L\ j,y ft J 


X, - X 


(jj^ w fx I j — u)-2_ ^Xq+ X, j 
-2£0j, w(xJ 


wfX|) - w(x) 

w(x,)-«i w(x,) + 

(xo + X|] 


and neglecting quantities of smaller order 


Pa'^p)^ y = (x,-x)[w(x,)-^^ix,)+wJx„^J]4flx(xJ-^^^^^ 




w(Xi) 

With the same notation as above, let Ci= ” » hj = — — , ct = ; this 

becomes 




(^O 
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With the abbreviated notation ri 


I } 

k times 




, we obtain 


If) -OC 

'a 


'i 




Taking the derivative of (2.2-4) with respect to C, and using Leibnizes 
formula, being the integrand, 

J I (i ^ ij - ( - "2 fio + ij ] 4^ 


Finally 

(2.2-5) 


The non-dimensionalization is completed by introducing the non-dimen- 

r2-11 

sional Etkin's number *■ 

^ 2 : 

where is the first cantilever frequency of the non-rotating boom. 

It is to be stressed that X is a constant only if o)^ , the satellite spin- 
rate, may be considered such. Equation (2.2-5) is rewritten in the form 


X = c( fa). 




41 


fa).r .1 




j - = - « 0)^ ( 



( 2 . 2 - 6 ) 
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So far, quantities which have been neglected were of order of 
smallness, or smaller. Now otm itself is of order e^, i.e. with 


dd) = 0) dt 
^ z 




0)^ 


if the product X the percentage change of (a^ per unit angle of rota- 
tion is very much smaller than quantities assumed to be of order e. 
Assuming that such is the case, we are then left with the homogeneous 
Equation (2.2-6) with a r.h. side equal to zero. 


- \ y ^ y 


= o 


(2.2-7) 


2.2.2 Solution of the basic equation 


d d 


Using separation of variables, with ' = ^ 


“C = 0) t 
2 » 




Hence 


I 

i 

j 

yielding 

f • - 


7 ■ 

Co^ 

where w. 

is t 


( 2 . 2 - 8 ) 


4)1 




ciated with (^o» ^)* This equation is in agreement with that obtained 

[ 2-11 

by Etkins and Hughes'” , in the special case Co = 
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Determination of w, (or w.) from Equation (2) proceeds as follows 

j J 

Equation (2) is linear, with ^ varying coefficients. Thus any linear 
combination of solutions of (2) is a solution of (2). 

Let^^ . be the solution satisfying the b.c. 

3,j 

: o o / 0 (2.2-9) 

f 

and/O^ j be the solution satisfying 


C = 0 ? 0 0 1 (2.2-10) 

Therefore, the desired solution, which satisfies the "built-in, free" 


boundary conditions 


conaxtions ... 

4- 4^y 4 <^) 

^ 7 ^ 7 

5-0 n o 


( 2 . 2 - 11 ) 


is of the form 


C ^^3 ; + X ; 


0 0 


( 2 . 2 - 12 ) 


(2.2-13) 


with C„, C, unknown. (2.2-11) is automatically satisfied by (2.2-13). 
3 4 


Expressing (2.2-12) 

i?l. 


5 

' y ^ "7 •* 1 =/ 

^ ^ ^7 -It./ 


(2.2-14) 


(2.2-15) 
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order to be satisfied for non-zero values of C^, should 


be such that the determinant 


i >4 1,J ‘4 <4 Jj , , 


(2.2-16) 


The successive eigenf reguencies , , are determined to any prescribed 

accuracy by iteration * j’^2 j determined by numerical integra- 

tion of differential equation (2), subject to b.c. (fl) and (/O.) 
respectively. 

The modal shapes (^) , which as expected are defined only to an 
arbitrary multiplicative constant, are determined, once is known, 
as 

1 (i) = C f i - J . 1 (2.2-17) 

<)■ 

2.23 Orthogonality of the mode shapes 

It is now proven, that given X, 2: 0, the modes <f) . , are 

J 

orthogonal, i.e. 

*.>•!' =« //I 


/ / 2 ^ 

„ « ” ^ ^ d7f 

Let oO be the operator 


12 ^ J ^ 




Hi,). 


X § 

K k 


(2.2-18) 
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Then, from multiplying by <|) and <j> . respectively, and subitracting 



- ^ li i ( '' i - ^ 

+ > ^ 1.J f $A~H f - 


Integrate with respect to from ^ = 0 (root) to C~ 1 (tip) > 



Terms corresponding to (1) and (3) will cancel in the difference. Terms 
(2) and (4) will cancel the terms resulting from the last term in in 
(2.2-19). 
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in the l.h. side of (2.2-19). Finally 0 



Again, terms (7) and ( 8 ) will cancel in the difference. Terms (5) 
and ( 6 ) will cancel the terms resulting from the last term in * 
in the l.h. side of (2.2-19). We are left with 







( 2 . 2 - 20 ) 


The modal mass, m. 

J 


(X ,^o)> is defined, for j=k, as 


in 


. " f 




( 2 . 2 - 21 ) 


in which normalized to correspond to a unit deflection at 
the boom^s tip,C = 1. The following quantity, to appear later, is 
also of interest 


. = [ 








( 2 . 2 - 22 ) 


with = Co + C, varying between Co (root) and Co + 1 (tip). It is 
readily determined when the modal shape, (C) > is known. 

Also, for later use, two identities are given here, which are ob- 
tained by multiplying Eq. (2.2-19) written for <J)., by <1), , and integrating 

J K. 


over the boom, 
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2.3 Equations of motion: meridional vibrations 

The developments in the case of motions in the (x,z) plane, of 
a boom located along axis +x in its undeflected position, or "meri- 
dional" vibrations, closely parallels those for equatorial vibrations* 
given in Section 2.2. In the following, only those terms which depart 
from the ones in Section 2.2 will be given in detail. 

2.3.1 Basic equation 

The equation expressing the equilibrium, at any section "x" of 
the boom, at point P, between the flexure moment from the left and the 
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moment , about P, of inertia forces imparted by the particles Q 
of the boom to the right of P (i«e. those having an abscissa x, 
between x and S. ^ reads 



(2.3-1) 


Now w(x) is an elastic displacement parallel to z. 
relevant quantities, 

(Xo + Xi)^ + 'z 

Thus 







Computing the 


a 


u^(xo+x,)-0;^w(x,J 
iHx,) - U)j ()^o+'^l) 


Also 


/t 


a 


ta A /c 










‘.V 


W(x,) +2U^, w(x,) - to, ) ^ Ay A>, w(x, ) 

-tlx w(x,) -2U^ +fVX|) (o^Uy fuj 410 

W (X|j - W(Xi)(ux ) + (xo + X|)(wx(02- 0 j J 
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Again, under the assumption of small displacements and transverse 
angular rates, terms of order w , are neglected, r^ 

reduces to 


X y X y 


a 


^2. 

wG,) +(mX()(wx(02-Wj) 


Along 1^, 


(vir>rt] . 
' a pv al 


x,-x 


W(y,)-W(x) 


Substituting into (2.3-1), and non-dimensionalizing 

H 


-H 

or 

%i " ~ J / ‘I h< '7 J *^’2 ( 1) ( % 

+ Uri)%\^^i ( 2 - 3 - 2 ) 

Comparing (2,3-2) to (2.2-4), it is seen that terms (b) and (c) in 
(2.3-2) differ in the following way from the corresponding ones in 
(2.2-4) 

♦ « 

(b) here has a factor (co w - u) ) instead of w 

X z y z 

2 

(c) here has a factor instead of nj- 

Therefore, with these changes, the equation analogous to (2.2-6) which 
describes the meridional vibrations should be 

'1 % [ I 1 G+4J 


(2*3-3^ 
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So far, quantities neglected have been of order of smallness, or 
smaller. Now, in order for the r.h. side of (2.3-3) to be of order 
we should have 


Ut cl<h 

very small compared to quantities assumed to be of order E. If such 


is the case, we are left with homogeneous equation 


% 




(2.3-4) differs from (2.2^7) only in that term - An of (2.2-7) is 
not present. 

2.3.2 Solution of the basic equation 

After separation of variables and non-dimensionalizing time by 


T = (i) t , the solution to (2.3-4) will be 
z 

n, = ^*,(4)T,(t) 

j 3 J 
sin sin 

in which T = m t = m.t, and satisfies the differential 
^ cos ^ cos ^ J 


equation 




(2.3-5) 


- % f • =0 

0 i 


As expected, this equation is the same as that obtained in (2.2 8) 
for equatorial vibrations provided the substitution of 



in (2.3-5) is made for 



in [2.2-8] 


(2.3-6) 
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Therefore, the method outlined in Section (2.2.2) to solve for can 
be adopted and followed without any other modification than that speci- 
fied by (2.3-6). In fact, program SEARCH DP, which obtains the first 
three eigenvalues 

given a pair (X, ^p), iteratively solves an equation such as (2.3-5), 

_ A ip 4 f'-V] ^ ^ Q 

^ ^ (2.3-7) . 

- 0 

in which the coefficient "COEF" is determined as follows: 



l I 


COEF = ~X(l+w|) COEFF = -Xw| 

2.3.3 Orthogonality of the mode shapes 

Modes d) (^) (i = 1,2...) for meridional vibrations can be proven 
J 

to be orthogonal, as in Section (2.2v3), since Equation (2.2-19) holds 
equally well in the present case. Thus 


i CoEF 


7 


boom 


<|> ,(|) dC = 0 j?^k 
J 


(2.3-8) 
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and we define, for case M, 


“l.j dif J 

2 

<l>. dC > 0 

i 

boom 

(2,3-9) 

®2,j dIf J 

d? 

boom d 

(2.3-10) 

with Co + C» 




With the substitution + (o? in (2.2-8) in (2,3-5), the follow- 

® j J 

ing relations, valid for meridional vibrations, are deduced straight- 
forwardly from Equations (2.2-23) and (2.2-24) 


for i^k 





(2.3-11) 


and for ,j=k 



(2.3-12) 


It should be noted here that for the same pair of values (A, Co)» 
if (COEF)j is the value to be given to COEF in (2.3-7), in order for 
the determinant (2.2-16) to vanish, then 


(COEF). ^ = (COEF). = COEF 
J,E j ,M 

„a,5o) + 1 = „(x.5o) 

J »E 3 


(2.3-13) 


whereas the modal shapes determined from (2.3-7) with the value 


(COEF) . of COEF have to be the same in cases E and M 

e )= ^ ^ ) 


In (2.3-13), if it is found more convenient 'to non-dimensional ize 
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by a quantity proportional to the 1st eigenfrequency of the non-ro- 
tating cantilever boom, namely 


then (2.3-13) becomes 


* 


U) 






(2.3-14) 




Nil 


lO 




as illustrated in some examples of Section (2.8) 


2. A Program determining the modal frequencies for equatorial or 
meridional vibrations: SEARCH DP. 

Program SEARCH DP, listed at 'the end of the present chapter, is 
written in FORTRAN V and implements the developments of Section 2.2 and 

2.3. 

The calculations are carried out in double precision, which 
suffices for values of X up to about 5,000. For higher values of 
X, an arbitrary N-precision, scheme had to be used: this is described 
in Section 2.7, 

2.4.1 Description of the program 

Number of statements (including cortim.^at cards) : about 270 

Input: - 1 card giving Q = E or M? ; X; Co format (Al, F6.5, G5.4) 

Output: 1) - A heading, specifying "Equatorial case" or "Meridional 


2) - The values of a "frequency" number" defined as/x 
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Lines giving the value of determinant of Equation (2.2-16), 
called here FE34 ; the value ofVcofr , the value of index U, 
number of trials in y before converging to the root of <^9 (lo^ ) = 0 


- Lines labeled KKK number of iterations, giving the successive 
values of the determinant as y is changed to obtain convergence 
of the determinant to zero. The iteration stops when 
differs from \x^ by less than 10 


*- A statement that "MU converged" giving the value of FE34 and y. 


- A*print-out of FE34, y, X, 


and MATFRQ, 


defined as 



3) 

2.4.2 


- The value of the step in y, DLT, and the value of the order 
of the eigenvalue, j or NOR 
same for j = 2, 3, in that order. 

Schematic flow chart: 

The following flow chart schematically describes the main con- 


trol flow in SEARCH DP. 
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2.4.3 Coinments 

\ 2“ A 1 

a) It has been numerically determined^ ^ that 100 steps across the 
boom’s length would suffice, over the range of X and Co investi- 
gated, to obtain eigenvalues agreeing up to the 5^^ digit with 
those obtained with 200 steps across the boom's length. The 

”l00-steps” are therefore incorporated as a "fixed" feature 
in program SEARCH DP . 

b) A method of linear interpolation is used for finding the roots 

of (to ) = 0. The iteration on y (or equivalently the eigenvalue 
J 

to be ) stops when two successive values of y, in the itera- 
tion process, agree to at least 0.1%, 

c) The integration method is a simple Runge-Kutta with fixed step, 
having a per step error of the order of Ax^. 

d) Using double-precision arithmetic, the number of significant 

digits retained in the two terms in £) , in Equation 2,2-16, 

does not suffice for values of X larger than about 5,000, and 

an arbitrary precision package ("NP" - package, N > 0 integer) 
had to be developed and is described in Section 2.7. 

2.4.4 Listing and sample output 

A listing and a sample output of program SEARCH DP are given at 


the end of this chapter. 
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2.5 Program Determining the Modal Shaptes <1)^ and "Masses” j » ^2 j 
MODE 

MODE is a Fortran-V, double precision program determining the 
modal shapes, normalized to unit deflection at the boom's tip, 

j = 1,2,3 

which are solutions of Equation 2.3-7, in which 

rj. is the jth eigenvalue determined by SEARCH DP 


COEF = (COEF) . ^ = (COEF) . ^ 
J,E J,M 


-2 






2.5.1 Description of program MODE 

Number of statements (including comment cards) : 158 

Input: - 1 card giving IE - E or M? ; j ; A ; i 

y = COEF (to be used in Equation 2,3-7) 
in (Al, II, 3G12.6 format) 

Output: 1) - A heading, specifying "Equatorial Case" or "Meri- 
dional Case" 

2) - The values of y^ = COEF^ (as obtained from SEARCH DP) , 
A, Kot j (1, 2 or 3) 


- The values of m 


IJ 


boom 


$ d^; m = 

^ 'boom ^ 


!2u_. 


m 


m 


l|j l,j 


which are of interest in the dynamical 
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simulation of the evolution in time of the space- 
craft angular rates (o) , m , w ) and modal coor- 

X y z 

dinates 

- The deflection a function of I, the 

station index, varying from 1=1 (at the root) 

I = 101 (at the tip), in steps of 2. 

2,5.2 Schematic flow diagram. The main control flow in MODE is 
as follows; 
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2.5.3 Comments 

a) The number of steps of integration, experimentally determined to 
give values of \i agreeing up to the 5th digit when solving the 
step size, was found to be 100. As in 2.4,3 and SEARCH DP, the 
the 100 steps are a fixed feature incorporated in the program. 

b) The method of integration is Runge-Kutta with fixed' step. . 

c) The calculations are carried out in double-precision, which 

should suffice for values" of X of up to 10,000. The data y., 

J 

however, might have had to be determined with the use of "NP" 
arbitrary precision package. 

2.5.4 Listing and sample output. 

A listing and a sample output of program MODE are given at the 
end of this chapter. 

2.6 Parametric Study of Eigenf requencies and Modal Shapes as a 

Function of X (Etkin^s Number) and (Non-Dimensional Radius 
of the Hub) 

Given the design parameters X and the study of the eigen- 

frequencies m. , (which normalized to m , are noted cj., and to 
* El V 2 

^NR ~ ^ — 4^ * noted ) will be made easier by using several 

“nr 

programs described hereunder. 

2.6.1 Preliminary Comment 

First of all, it should be emphasized here that there is no 

point in comparing mode shapes (j> . for "E" and , since they are 

the same solutions to Equation (2.3-7), for COEF, = COEF = COEF , 

3 jjE j ,M 

once j has been chosen and X and have been given. Any slight 
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numerical departure, such as described in Ref. [2-4], 2-5] could only 

result from the inaccuracy is determining the eigenfrequency (0.1% 

relative accuracy on p , in program SEARCH DP). Only the eigenfre- 

quencies m. _ , ^ corresponding to these modal shapes will be different 
j ,E J 

2.6.2 Program computing dynamical parameters, given X,Ko' PARAM. 

Program PARAM, written in FORTRAN-V, will permit to get a quick 
look at various relevant dynamical parameters, given Q E or M, X 
and ^6, namely 


(m^/m^)^ and j = 1>2,3 

(m?/m ) . 

2 1 3 

and also the sum over one, two, three modes 

I i 

a quantity to be used later in this work. It will also plot the 
mode shapes (up to j = 3) in the computer printout. 

The data entered are 

w. (j = 1,2,3) obtained from SEARCH DP, 
J case M (NDS = 0) 

X 

The program basically computes (O and 
the relevant integrals, , m etc... as defined before. 
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A listing and a sample output of program PARAII is given at the 
end of this chapter, 

2.7 Arbitrary Precision Package: MP (for use on OS) and P (N- 

Precision Package), in Fortran. 

2.7.1 Motivation 

[2-4] 

An earlier version of SEARCH DP had been written to alle- 

viate a problem of numerical stability at large values of X (higher 
than about 5,000). This version used on IBM-library multiple precision 
(MP package). It was found, however, that this package was unavailable 
in a TSS environment. Therefore, an arbitrary precision package (NP) 
was written in Fortran V, and used for finding the eigenvalues 
at values of X, and the accuracy of determinant in Equation (2.2-16) 
will be critically atffected when taking differences of very large 
numbers . 

MPAP (Multiple Precision Arithmetic Package) is present in the 
Internal Library of the IBM-360. The routine calls on specialized 
subroutines to perform floating point calculations with precision 
to be specified by the programmer (typically, here, quadruple pre- 
cision was required) , 

MP-SEARCH, as used in Ref. [2-4], and MP-MODE,^^^ thus basically 
MPAP versions of SEARCH and MODE. Their one disadvantage, As expected ,15 
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an increased running time, of the order of 1»5 minutes for eigenvalue 
(IBM/360). For this reason, it is important that the eigenvalues or 
modal quantities of MP-MODE obtained for high X be stored for later 
use in the simulation (Option MGIV = 1 in program FLEXAT , see Chapter 
5), and that interpolation be used whenever possible. 

2.7.3 Multiple precision in TSS: OT-package 

Written in FORTRAN for case of conversion to any machine, N-PRES 
is a multiple-precision arithmetic system for scientific calculation 
It may be used on any machine which stores one integer per work, where a 
v 7 ord is i 31 bits long. 

2.7. 3.1 Short description of the program 

2. 7. 3. 1.1 Representable numbers. 

Let N, M be integers 


2 5 N < 16 

All numbers in the program are considered floating point constants of 
+ N precision, expressed in scientific notation. Thus, for N - 3, or 


precision 4N =• 12, we could have 
. 371246875003*10**8371 

The exponent must always be an integer, positive, negative or zero 
and less than or equal to 4 decimal digits long. Thus a number such 

+ d, drw • . * d, 


as 


- "12 


‘60 


* 10 ±^ 1 ^ 2 ^ 3^4 


2.7. 3. 1.2 Internal Storage (Multiple-point, floating) 

The mantissa is stored 4 digits to a word, in N digits. 
The exponent takes up the N+1 tion (Any 1 < N ^ 16) 
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Example : for M = 3, N = 2 

8.4326 X 10^ = .84326000 x 10^ 

6.0 X 10^ = ,6 X 10^ 

represented as 


NUMBER 1 


NUMBER 2 


8 

4 

3 

_2 

6 

0 

0 

0 


EXPCH Q 

0 

5 


6 

0 

0 

Jl 

0 

0 

0 

o’ 

0 

0 

2 


Word 1 


Word 2 


All operations are designed to handle such units, called W-CONS (for 
N constant) . 

2. 7. 3, 1.3 Quick guide to operations and subroutines 


Name 


Subroutine Function (all operating with N cons) 


IN IT 

INPUT 

OUTPUT 

CIN 

CNI 

CFN 

CNF 

NABS 

NPm 

NSCL 


Initialize the N-precision system 

Input 

Output 

Convert integer to N-CON 
Convert N-CON to integer 
Convert floating point to N-CON 
Convert N-CON to floating point 
Mem (Add) ^ ABS[Mem(Add) ] 

Mem (Add) [Mem(Add) ]**P 

with P a parameter to NPWR 
Mem (Add) [Mem(Add) ] *10**S 

with S a parameter to NSCL 
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Name 

Subroutine function (all operating with 

NCMPR 

if Mem(Add 1) = Mem(Add 2) , A = B 


COPY 

if Mem(Add 1) > Mem(Add 2) , A > B 
if Mem (Add 1) < Mem (add 2), A < B 
with A, B, parameters to NCMPR 
Mem (Add 2) ^ Mem (Add 1) 


RENORM , 

SHIFT ^ Internal use only 


PUNCH 

Output to punch 


IMUL 

Mem (Add) Mem (Add) *I 


FDIV 

I = 1 integer 1 < limit 
Mem (Add ) Mem (Add ) /F 


MADD 

F = floating point 
Mem (Add 3) ^ Mem (Add 1) + Mem (Add 

2) 

MSUB 

Mem(Add 3) Mem (Add 1) - Mem (Add 

2) 

MMUL 

Mem (Add 3) Mem (Add 1) * Mem (Add 

2) 

MDIV 

Mem (Add 3) Mem (Add 1) / Mem (Add 

2) 

2,7. 3.2 

Some examples of N— precision programming 


2. 7. 3. 2 

Square root 



A. Algorithm ; Newton-Raphso.n 

\ 

Let B = '/k , with old B = 1 

If Abs(01d B — B) > B*10** limit 
Old B = B 


Else done 
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B. Fortran Progarm; 

Limit = - 12 
Read (5,1) A 

1 Format (F10.2) 

Old B = 1. 

2 B = (A/Old B + Old B)/2 
X = ABS(01d B~B) 

Y == B*10**Limit 
If (X, LE, Y) G0 T0 3 
Old B = B 
GO TO 2 

3 WRITE (6,4)B,A 

4 FORMAT , F10.2, ' IS SQUARE ROOT OF, F10.2) 
STOP 

END 


C . N-Precision Program 

IMPLICIT INTEGER (A-2) 
CALL INIT(1,4) 

CON V = -12 
A = 1 
B = 2 

Old = 3 
TWO = 4 
X = 5 
Y = 6 
.HALF = 7 
CALL Input (A) 

Call NSCL(A,1) 

CALL C IN (TWO, 2) 


Comments 
(all N-cons.) 

(16 digits of precision, N=4) 
(limit) ’ 

(Allocation of variable 
names to N-con addresses) 


(N-con at address TWO contains 
the value 2) 
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Call CIN(01d B,l) 

Call CFN(HALF,.5) 
Call Output (HALF) 

Call Output (A) 

Q = 1 

Call to IV (A, 01dB,B) 
Call MADD(B,01d B,B) 
Call MMUL(B,HALF,B) 
Call Output (B) 

Call HM 

Call MSUB(B,01d B,X) 
Call NABS (X) ' 

Call COPY (B,Y) 

Call NSCL (Y,CONV. ) 
Call NCMPR(X,Y,I,J) 


If(I. LE, J) GO TO 2x 

Call COPY(B, Old B) 

GO TO 1 
CONTINUE 
Call HM 

(Write out results) 

2. 7. 3. 2. 2 Conversion of a stateirient 
Consider the FORTRAN statement 

IF(FE34aDECID) 52, 51, 50 - 
The N-PREC. translation would be 
CALL MMUL(DECID, FE34, TEMP) 
CALL NCMPR (TEtlP , ZERO , I , J ) 
If(I,LT. J) GO TO 52 
If(J.EQ. J) GO TO 51 


(N-con. at address 'Old B' 
contains the value 1) 

(Half contains 0.5) 

(Conversion OK; print 
and check) 

(Print input number) 

(Iteration Counter) 

(B - A/Old, B) 

(B = B+Old B) 

(B = B*.5) 

(write partial answer) 

(How many subroutines called 
so far. Print it out.) 

(X = B-Old B) 

(X = Abs X) 

(Y - (B)) 

(Y - (Y)*10*:«c CONV ) 

(Result : 

If X > Y, I > J 
X < Y, I < J 
X = Y,' I = J) 

(IF (ABS (B-Old B).LE. 
B*10**CONV) qq 2 
(Old B = (B)) 

(Done! ) 

(How many calls) 

from SEARCH J)P 
of SEARCH DP: 


50 CONTINUE 



2.8 Results from programs SEARCH DP, MP and NP 

The frequencies o). (normalized to o) =1) for j = 2, are given 
J s 

for case M. Those for case E are immediately obtained from 




-.2 


Q m 


Also given below is the quantity , which will be of special 

importance in Chapters 4 and 5. The first non-dimensional frequency 


/ pit VEI 


is also represented, for cases E and M, and various 


values of on Fig. 2.2. 
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CASK M - FIRST NONDI>IENSIONAL NATURAL FREQUENCY 



o 

o 

• 

o 

0.10 

0.25 

0.50 

0 

3.681 

3.703 

3.734 

3.787 

5 

1-913 

1.953 

2.013 

2.107 

10 

1.555 

1.605 

1.675 

1.788 

20 

1.339 

1.395 

1.476 

1.601 

30 

1.256 

1.316 

1,401 

1.531 

50 

1,183 

1.246 

1.335 

1.469 

100 

1,120 

1.186 

1.278 

1.417 

200 

1.081 

1,148 

1.242 

1.385 

. 500 

1.050 

1.118 

1.214 

1.358 

1000 

1.034 

1.104 

1,201 

1.346 

3000 

1,021 

1.091 

- 1.188 

1.339 

7000 

1.016 

1.087 

1.184 

1.329 

10000 

1.013 

1.083 

1.181 

1.327 





NOTE; 
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CASE M 


1 

5 

10 

20 

30 

50 

100 

200 

500 

1000 

3000 

7000 

10000 


- SECOND NONDIMENSIONAL NATURAL FREQUENCY r j 


0.00 0.10 


22.18 

22.20 

10.236 

10.276 

7. '*19 

7.476 

5.596 

5*624 

4.760 

4.849 

4.023 

4.128 

3.364 

3.486 

2.976 

3.113 

2.707 

2.855 

2.603 

2.754 

2.520 

2.673 

2.490 

2,644 

2.482 

2.635 


CM 

• 

0 

0.50 

22.23 

22.78 

10.339 

10.447 

7.561 

7.703 

5.736 

5.921 

4.981 

5.191 

4.281 

4.523 

3.665 

3.941 

3.308 

3.606 

3.060 

3.373 

2.964 

3.282 

2.886 

3.195 

2.857 

3.178 

2.849 

3.171 



NOTE: 
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NONDIMENSIONAL DYI^.'-JIICAL PARAJIIOTERS 


S\m OVER 3 MODES 



0.0 

0.10 

0.23 

10 

0.3328 

0 ,4A01 

0.6336 

100 

0,3329 

0.4404 

0.6343 

1000 

0.3332 

0 . 44l 3 

0.6366 


A = 0.3333 

A = 0.4433 

A = 0,64.38 



0:00 

(ONE mode) 

0.10 0.25 

0.50 

0 

0.3233 

0.4190 

0,5810 

0.9250 

10 

0.3249 

0.4212 

0.5893 

0.9325 

20 

0.3260 

0.423! 

0.5929 

0.9400 

30 

0.3268 

0.4246 

0.5957 

0.9457 

50 

0.3280 

0.4268 

0.5997 

0.9530 

100 

■ 0.3297 .. 

0.43QI 

0.6056 

0.9652 

200 

0.351 ! 

0.4330 

0.6111 

0.9757 

500 

0.3323 

0.4357 

0.6165 

0.9861 

1000 

0,3328 

0.4371 

0.6194 

0.9916 

3000 

0.3331 

0.4385 

0,6224 

0.9950 

10000 

0.3332 

0.4392 

0.6241 

■ 1.001 


A = 

A = 

A = 

A = 


0.3333 

0.4433 

0,6458 

1 .0833 
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PROGKAM LISTING 
AND 

SAMPLE OUTPUT 
SEARCH DP 



1. 


t- i IIA I I u I . 


t i ri S * PR O G f .'^1 ’ p 1 0 S "Th E 'f 'ii <’S f “ T HR K £ u 1 0 E V A L UE S ( MU ) FOR T H F ROT A T I H G 
BOOM If^ EITHER EOUAfORIAL OR MEk luION, FLEXURE • - - - 

...l.T... W I LL COMRU 1 E. TMEGE ...accurately,. F0K....V AUUES...QE....LAMDA...UP_jrO 

APPROXIMA I ELY BOOu 


OoUciLt^ FRECiSiOH P(^) (4) » M ( <4 ) ,L(4) fE5{101) #E31P(101) »E32P(101) ► 

1FE34 t UEC I D f I VAL f H Y r f EMOM4 r lS. tR3 f ESMR4 t f EPRV r .. — 

iE3.3P(10l) »E34P(ini) rEA(iOl. ) rF 4 lR(ini) fE42P(101> rF43P(lnl) » 

3E44P(irjl).f AfBfCrEf A0rb0.rCu».L0f MUl_r..MU»EP.SBrEP5C»LASrLASS» 

4 ULTf UPrD/RJrSiiOfSl 

E A L M f f'4 M f L A M f M A TE R 0 ... j 

Integer IrOfUrZ/ MiHT»lMTERrwrKflUNORrKKKrQ 
UATA/LM/1H.V1 

If iPUT_..DAT.A„ 1 S...E_.OK_..i^ ?.L AMBU.A., PS I «ZEao. 

,la„ READlSr2UrbHL) = ;i:2)0»LAMf S.IO. ... .. - 

20 F0RMAT(A1 »F6.5f GS.4) 

lF(Q.EG*LM)Mi)S~0 

i F ( NpS * Eu • I ) ‘.vR i TE < P ' 1 : 

lF(HUS*Ew*u) .vRITE(bf 17) 

1.6 FORMAT ( ♦ EOUaTORIAu CASE,* //) 

17.: FOKMmK’ Mc.R1D10.mAL CASE’//) 

’ .WRUEtb^Zl) LaMiSIO , 

21 Forma K lH » ’LAM-» » F12.6» 3x > *SiO-’ »D9.3/) 

R = 0 

FE34~0. 

MU—l^OD-A 

’ ■ ~ EPSB=1.01J-I4 



feprv-6 . 

EPSd^lO * **~I4 ■ ■ 


. .SET ..NOP T, = ,_.l . .FOR ., KE y Ea SL.r)_.I n.TEGR A T I ,oN_...( J, IP ..._T.0._ RO.OX)___ 
NOPT^i ' 


NDS =. UIi<ECTiOi'j SWITCH 

wiiEn nDs ~ seakch For cOuAroKiAL Poors f 
when NQS ~ i)r SEAkCtI FOK MER IL)I AM. Ru01 5 


NOR^_ i 

'n-l 

olt-1* „ _ „ _ 

IA = 0 

NlNf-inO _ 

InTER-NINT+I 

-(^.NFjiiQ-SOR I I.LAm.) 

f } R I T E ( I'l » 6() ) A i-iF R 0 

►n Forma! (ih t *nfrq=^sort eaM-* fFio*5) 

>9’ Si::0. 

T M AL— F E 3.4 

Jj"u 

clear arrays 

Do"3i I-liT” 
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1K-. 

llt>r 


K U ) =0 . - 

L(I)=U. 

Mil ) = U • 

"^^31. P(I)=U. 

Do 1 l = lri01 

E34PU)=0* 

E44R( i )=0. 

E33P(i)=0. 

E3;>PU ) ^0, 

E4^^P(1):^0, 

E31P(1)=0* 

E4lP(i)=0. - 

E4(X)-0. 

1 1)^-0* ... 

H“1 ,/hUOA l "(iMlNt ) 

~C SET Tf^IilAL CONofTi'o^ ON THE S3 A>si'[) S4 SOLUfior'S 
D-3 

ft I.p_^ D E Q D 0 TO ; I. 

EO-U • 

Bu=U,v , „> — 

IF (nOPT^GT .U) GO TO 12 

E32PU)-1. • 

BOr:E32PU) 

__ _ Gu,j;0 ,.J3 

' i? E3(i)=i. 

E0=E5P1) - 

13 CO^U. 

Au~u . ; : 

GO TO 3 

2_ AO-0* _ 

co'=o.’ 

IF (N0PT*GT ,0)GO TO 14 

E43P( i)=l . " 

C0-C.40P (1 ) _ 

GO to’ is ' 

14 E4dP(l)-l* , 

'''~Au=L4ipTiT 

15 Bo=0*^ ' 

EO-U • 

3 A~A0 _ _ . 

B:;B0 

_C=C0 

Y-EO ' • 

■■ C "beg i i i ' RUNG A‘ kuTT A I ,4 TEGK' aT! ON 

N=1 

Ni\|E:(4 

S 1 = t NN - 1 . ) ^ 

S R ( I ) -H+A 

+ B _____ „ - 

■■■■ 

_ MUl-l.+MU + MU __ __ 

■'if(nds*cQ*'i)''go' To •air*'’ 
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MUl-MUl-l. , . . 

40 P ( I ) = ( ( 1 .~Si*Sl+?.*bIU* U .-ST i )/2.+b-(Sl + SI0) ♦/^+^^u■3 -t^E) *LAm^H 

1 F ( isiOPT . GT . 0 )... P..( I ) = ( .( -S I *S I +2.. +S I *U . +S 1 0 ) )./2 . tB 

1+ ( 1 .-Sl+SiO) *A + .viUl*c:)=*^LAM.KH 

Sl=(NiM-l . )*rl, 

1 = 1 + 1 

- IF ( 1 .bT.3)_.bO -TO. 6 — - 

1™/+ Zri-i 

-Vdpfi £=E0tK(Z)/2. 

A=Al>+L(Z)/2. 

im* B::B0+tA(Z)/2. 

I3i+ C=C0+P(Z)/2. 

l|p 5 1 = S I + H / 2, . 

l.|p+ GO TO 5 

1 ^4.+ 6 I F ( I . GT ..!F).„G0„X0.„7 

t=EO+K(3) 

X K * A=A0+L(,5) 

iw+ B-BO+Fi ( ^) 

X5o^ C=C0+P.(5) i,„ 

1#-^+ S1=SI+H 

aM>L+ -..-...-Go . „T.O.,..S 

XHlr 7 IF(U*E0.4) GOTO 9 


„..S 1 =n.n *H 

Z=N+ 1 

: £S(Z)=F3(.si) + (k ( 1 ) +2*+k (2) +2*+K (o) +■< Un/6* .. 
Ej»1P(Z)=Eb1PIi'J) + (l( 1 ) +2.+L (2) +?.+i-(o)+L (4) )/6. 

...E32.P ( Z ) =ES2»Z( r^^ ) t ( .'1 (1 1+2 • + ( 2.) +2 . 

E^3PI/)=Eb: 5F( l\!) + (P( 1 )+?.*P (2) +?.+F( o) + P(4 ) )/6. 

. E.^4P(Z)=LAr-'i^( ( (SIO + 1. ) ++2- (Sl+SiOi *»2) *E32P(2) /2. 

i-(SI+Sl0)+E.51P(Z)+MUl*F5(Z) ) 

IF <HOPT*GT. U ) E54P(Z)=LAM+ ( ( (siu + i. ..> + + ( 1 • -SI +SlO) *.+ 2) 

i+E:52PlZ) /2. + ( 1 .-5X+blO)^F:=^lP(Z)+Mbl + E,i(Z) ) 

„.EFFp ^ 

A^E31P(im + 1) 

3 = FB2P(N+.l ) .. . 

C-Eo3P ( N+ 1 ) 

Eo.“E 

AO=A 

...Bo=b ; 

CO=C 

. N::N+.1 _ 

IF (lM.LT, INTER i GO TO 4 - . - - 

EviOM3=Eo2F( J.NTEH) 

£SHtU = E03P (INTER) 

_IF (NOPT.G r .u) EmOm3=Es3( INTER) 

i fTimoht . G r . 0 )' YsVir:3=e31p TinVe^^ ) 

0 = 4 _ __ _ __ 

'"go’ to B ’ ' ' 

SI=NN^H 

Z-N + i" 

E_4 ( Z ) =E4 ( i\l ) + ( K_( 1 ) >2 • + 1\ (2 )+2.+K(0 ).+^ ( 4 ) ) /6 • . _ 

*’"E4 1 P U ) =E4 ip ( N ) + ( L ( X ) 42. ( 2 ) +2 . +L ( 0 ) +T ( 4 )T/67 

E42P(Z)=Ea 2P(N) + (i^i(i )+P.*i-i(2)+2.4i-.(o) +M(4 ) )/6. 

E4 3P ( Z ) =E45P ( ) + ( P ( 1 ) +2 * +P ( 2 ) +2 . '^P ( 5 ) +P ( 4 ) ) /6 

E44P(Z) =LaM* ( ( (SIO + X. )+^2-(Sl+Siq) + 42l*E42P (Z)/2.. _ 

X-(Si + blO) +E4iP( Z) +?^Ul*E4 (?) ) 

. , I P .( iM 0 P 1 • G.r . ,u ) E 4 4 P ( Z ) = L. /A M + ( ( ( S T O + 1 • ) + * 2 - .(. I S 1 + S XO.).* +j?J. 
1^E42P(Z) /2.+( 1 .-Sl+SIO)4E4lP(Z)+MUl+E4(Z) ) 
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i:: 

I7v* 
* 
•*: 
:♦ 
■* 

* 







1^0 ♦ 




7* 
199T 









2! X L > ^ 


•T' 









E=F:4{r,i + l) 

Aze4ip{N+i ) ■ ' * 

B=t:42P((N+ij , 

C=E43P(N-^ll " 

JE0=t 

A0=A 

130 ::b __ 

co~c 

N = 

IP(^'j»LT♦Ii^J^c,R) GO TO 4 

E mom 4 - e 4 2 i M r EH ) 

ESHK4~E43P ( INTER ) 

iFCNOPT^^Gr.U) EM0 m4^:E4 (INTER) 

IF (NOPT.GT ,U) ESHR4~E41P ( (NTER) ’ 

C 

C RUNG A KUTTA' fTniSH'ej 

C NOW RcGiN linear IN rERPULATtON 

C FE:34 is the V/\LUE of IHL determinant ' (~sY a^ 

C ..... • _ 

FE.44=EM0M.3 + ESHR4"'ESriRo*EM0M4 

lF (K*ci'o« l ) G0_ ro_3.l ___ 

lF(U.c,0.1) GO TO gO' 

lF(FEo4x0t-:ClD) 61 r5U 

!3n DECID-FE34 ' ' ■" 

LaS^MU _ _ 

LASS-LAS 

_ „ R ^ r L ( 6 » A 6 ) F E 6 4 ^ M L) > U 

■ 8s Format c ih f ♦ fe34 * f oi? • 6^ 5x » • mi jr » > * o » sx ? * Ti 3) 

MU=MU+DLT „ „ ^ 

U-U+1 

Go to 99 _ _ _____ 

'sp UP^^^MU’ 

OivN^LAS 

~~’HY-FE64- 

TVAL=FE64 

SI IF (aBs(FE64) .lE’.EPSg) GO TO 53 

IF (ABS (OEClu) .L£.cPSB> GO TO 4^?_ 

'Vkl lEIGf Brt)RKK»FEo4»MU 

.BA Format (IM r »KkK » r I3f 3 a» »F|-34’ »0l2.G»3Xf »MU» f D1P.6) 

R=1 ’ ' ' " 

IF (FE64*0ECID) 55»61 » 56 

”55 up^^Mu ; 

MU-D Wi'i- ( DvvN-UP ) 'I'DLC i 0/ ( DEC I D-FE 64 ) 

KKK = Kt\K + l 

EPSC=aBS(mBS(mU)~aBS(LASS) ) 

IF (lPSC.lt «iviU + lL) . ^*-4 ) GO TO To ~ 

_ L A S S — i 'I U 

Go TO 58 ' 


56 OwN-MU 

0ECID=:FE3i4’ • 

a"v» Ml jLDWH- ( DWN-UP ) ^DlC 1 D/ ( DEc TD-HY ) 

KKK-KxKxl 

^ EpSC“ABS(/vBb(MU)“Ai3s (LASS) ) 

IF (EPSC .Lr .mU + 1u .>#=*“4 ) GO TO li) 

232+ LASS=MU 

" 58 FEl = AbSCFL-.PKV)‘-ABS(FE64 ) ■ ■ 

2 ^1^ *■ I F ( A 1.'} S(Fc,l) .GT .1 .Ou — 14 )GOTO 8 2 

236+ I A = 1 A + I 
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^"b+. 
23 /* 

I ;y.^ 


mr-* 

243* 


27 u ♦ 
Ml* 


Ir (lA.LI *h) GO TO 82 . - - 

WRITE l6fB3)FF-34 

83 H.O)?MA r ( IH... f.!.STUCK...0rL-.r.l:ll5...F£343.rD.iH •. 6J — 

IA-0 

FEPF<V-0. - - - “ 

GO TO B3 • 

... 82 FHPKV-Ft£34 _ - - --- - ■ ■ ■ 

Gu TO 99 

42 WR J TF ( h > 43 ) uFCID ^ 

43 Format (IH »*N0 good Dc-CID::* »D12.6) 

Go TO. 53 ... . . . • - 

in WR I 1 e ( 6 f 1 1) F F34 f MU 

11, Forma r ( IHO ^ * MU COiMVtROED .FE34=» r 024.1b »3Xr*MU-» »012.G) 

53 N/XTFRo — MU + bvvHT (LAM) 

fj\i IT £ ( 6 » 54 ) F £34 r MU f f UaTFRQ f. w — 

54 FORMAT (lH,)f »F£34* #Di2.6»5X» *Mu=* »ul2.b*bX» *lA|v'<-» rF1.2*6» 

15X ? ‘NATFRo-V^El^.orbX » ,Vw=» r 13) ■ - 

WRlT£(6r8f.) DoTrUOR 

^..8b._. Forma t (lHtM ,»DLT-’ > 09.3r.bX r..TN0R=*.?..i 3) 

IF (inIOk.lQ.3) go To b7 

R = 0 

KkK-1 

GO TO 99 

vvRiT£(b^iOu) 

lOo ... FORMAT (.Vl_.V.) ^ — - - - 

UOR-1 

DLT^lV 

... 

GOTO . lb - 

22 CALL EXIT 


F OF 1 1 H I V A C 1 1 0 b __ F 0 R J \i AN _ y _ C_0 m P 1 L A T_I_0 N • 


0 tOlMGNOSTlC'fr MESSAGE(S) 
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rttTATORXAL CAm£ 


Af 


iO.UUOflOO 


SlO= *100+000 


r 
c 

K] 

K 
KK 
K 

I 

K 
K- 


i 

^ ... 

^ tC 

I 


=f;QRT LAM- 
.. . 17u lAU+nOl 
.&07?bU+o00 


3*16228 

■ MU=„ 

MU- 


. luOOOO-uOb U-.. 1. 

.lOOuOO+UOl 2 


J 


. E £ 64 - ,, 2 3,9 b4 8 + 00 1.. 


!U •a.UOOQn.t.OOi 


2 

r £34 

.134973+000 

Mu 

* 120223 + 001 

..3 

Fti34 

•2589 76-0 01._ 

Mt.J 

.124479+001 

4 

F£34 

. 482325-002 

Mi,J 

.125280+001 

.5 

Ft 34 

*893275-003.. 

MU 

•12543o+U01 _ 

6 

F £34 

.1652o4-a03 

MU 

*125464+001 


:U 


OiJVtRGEij' FE34 'i:' • lh52~6’4 1 o”682VoTi 1^5'4 69+^^ 

L A rd To do do d 


c.*>4 .166264-003 


m‘U= . 125469 + 0 6 1 


NATF60= 


i: 

sfi 

£.'•4' 

at 

I 


.luu+udl " fsioR 

• .„5 3 7 « 2 2+0,0 1. 

' . 7d6‘^*32+»i01 
‘.1 11908 + 002 
^ • 1 26v.»38+0 02 
lUU^; 39 + 0 02 
'.165669 + 001 

p FE34-. . 396210 
3 F£ 34 -.798063 
*4 F£o4— ,159426 

FL34-, 317792 


- 1 

225469 + 00,1 

M ij- *325469 + 001 
MtJ- *425469+001 
i«4U~ *525469 + 001 

mU= .625469 + {)01 

tViU- *.725469 + 001 
+ 002 _VUt_ . 825469 + UOl 
+ doT'"" Fiu *♦ 737738+0 01 

-001 MiJ *740274 + 0 01 

— OUl ^'!iJ *74078^ + 001 

-002 Mu *740685+001 


uB.0 

£.^,4-.3l7 792-n02 


t 

:.'j4- 

I 

£o 4 
r, 5 4 


.vVERbL'u FE34-- . 317791B321 7501Q385-002 '-iU= *74090 d + u01 

MU- *740905+001 LAM= lO.OOuOuO 


WATF 60 — 


t 7>- 


-OH 

I 

:o4 

I 

vl\ 

I 

C( 

I 


.100+001 UOR- 2 

,162 6 66 +n U 2 Ml.J- * 84 0 9 0 5 

,4.U3 l99+nu2 MU- .94 090 5 

*.7i89bb + n02 * MU= *lU4u91 

.109/83+003 MU- .114091 

.l5lbbb + il03 MtJ- *124u9i 

.193483+003 MU- .134091 

.230505 + 003' ' MU= *l44o9j. 

,256265+003 MU- *154o91 

. 2b32bo + nd3~~ ""’mU- . 164u9l 

.243li54 + f]03 MU= *174091 

,186657 + 003 MU- .18M(j9i 

,849254+002 MU- *194091 

1 FL34-. 7086 ^5 + 002 Mt) * 

;> ,,Fj:' 34 .722558 + 001 Mii . 

8 - £34' , 5 1 9 3 8 2 + n U 0 ” Mi li * 

4 Ft.34 .368619-001 MU * 


+ 001 
+ 0 01 
+ 002 
+ 0 02 
+ 002 
+ 0 0 2 
+ U02 
+ 002 
+ 0 02 
+ 002 
+ 0 02 
+ 002 


U= 

u- 

u- 

u= 

u= 

uz’ 

u- 

u=’ 


1 

2 
'3 

4 

5 

6 

7 

8 

’"'9' 

10 

11 

12 


204091+U02 
1 9954 2 + 0 02 
199963+002 
199993+002 


u COtjVEKbEU FE34:: . 36861886235783458 /-OOl mU=: . 1 999^'»5 + 002 

.36hnl9-001 MU= .199995 + 002 uAM= lO.UOuOUO NATFr?o: 


. 396769+01 


*234295+02 


.632439+02 


lli'^iOlONAL CASE 


lU .OUOOOO .... . SlO= »10U + U0u 


S U^SOKf LAM:;, j;5.,162;?8„ 
^ .284:^9^ + ^)0i .viU- 

’ . .. 1.7 U. 1.7 9 + C\ 0..1 fvU lr..„ 

1 f- l54-«, 144094 + 001 

?. .207673+000 

3 FL34 .173932-001 

1 ,..._4 FL34 .140369-002 

S hL34 .112971-003 


lOOuOo-OOb U: 

, 1 .0 ,P 0 0 u + 0 0.1 Uj 

.200UOO+001 
... Mi» . 164150 + 0.01 
MU .169920+0 01 
Mu .160404+001. 
Mu .160442+OOi’ 


i d j\| V E R u £ U F O 4 = . 1 f2 9 7 d 8 5 7 1] 2 3’6 3 6 5 1,7'- 0 03 w'O = . T o d~4 4 6+001' 


E.>4 .112971-003 


MU= .160446+001 


lam- lu.oOoduo 


NATF^Qr ,507373+01 


■ , 1 0 UT- d 0 r NO R = ' ' i 

; - . 3 9 b 4 3 5 + ft 0 1 M I J 

”•641747+001 MU- 

117975 + 002 MU= 

”.123912+002 MU- 

”.o426o8+001 MU=: 

I FL64 .166603+001 
„ 2 F £ 34-, 48 37 14 + 000 

3 FL344, 17^3^35-.qU X 

4 FE34-. 623739-003 


NOR= 1 

M 2 6,0 4 4 b + 0 p ,1_„ U 

MU- .360+40+001 U 

Ml U— .460445 + 001 _ U 

MU- *560445+001 U 

MU— .660445+001 U 

603+001 Mu .760446+001 

714 + 0 0 0 Mu ..7437 73 + 00 1 

355-001 Mu .747490+001 
739-003 _ Mu .747623+001 


F£346-.623739425547498172”U03 MU- .747o26+001 


.623/39-003 


,viU- *747628 + 0 01 


lam= lo.ooonoo 


NATFpo= *236421+02 


.100+001 NOR- 

.164d05+O02 mu 

.407204+002 MU 

,72b7ou+no2 mU 

.110/35+1)03 iV|U 

.1d2709+003 mU 

.194097+003 MU 

.231.051+0 03 MU 

,255671+003 MU 

♦ 2530 4 8+0 0 3 MU 

,241r3J8 + nu3 MU 

. 183:>0o+n03 "mU 

* 795^£/5+('|(j2 ,Vi(J 

1 H £34-. 788630+01 

2 Ft-: 54 .732752+01 

3 F L 3 .5782 6 2 + 0 L 

4 Fl34 .450 579-0 L 


'847o26 + (j0l L 

'947o25+00i U 

104763+002 'U 
1 14 763 + 002 _ U 
124 763-^0 02 U 

134763+002 U 

144765+002 U 

154763+002 U 

164765+002 u 

174765+002 _ _ U 

184 765 + 002"’ ’u 

194763+002 U 
Mu .204763+002 
Mu .I9978of002 
Mu .200209+U02 
Mu *200242+ U02 


In 
u= li 

u- 12 


J®;0>vyER6EU FE54- .45057870441269187 /-OOl Ml.l- .200245 + 002 


MU= .200245+002 


LAM= in.oooouo 


NATFo{)= .633229+02 


. lOu+001 


PROGRAM LISTING 
AND 

SAMPLE OUTPUT 


M0DE 
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mode 

_THIS program calculates tHe mode shapes 

gTvEN'the ElGENVALurS,LAMBr)A»' ANrr'PSI^ZERo" ' 

reversed integration methOo only, case e or caSe m -specify on input 



implicit double PRfC I S I on ( a-H »0-Z ) 

■ integer DtZl ” ‘ ' 

double precision Mu*LAM,MUi 

‘ Fn U 0 L E P R E C 1 S I 0 N ' K ,l-\ 7% , P ) 

p I mens ION EG3 (ion ,Eg^ ( 10 n ,nBPT (1 n I ) .BPT < 101) 
UaT a/LE/7 he 


^ „ NpRfal FOR printed qU t PU j / NP R t = 0 F 0« PUNCHED OUTpUT 

1 “ Read ( ^ 1 1 eo ) ' NPR I 

IfiOEoRMAT(II) _ 

C ' a K E " S U R E ' 'm u'"'" i S t H E” c 0 R RE C t 6 N E ' F V, R*"' E l T HE R "'Yh E'^'ET R ' n'T A 's E V 




^ E A D ( s » 9 3)1 E t j z »„E A M , S I D » M U 





NN=FLOaT(N) 


5 j 3 ( f-JN- 1 . ) ♦HH 


. K < f l=H* A 


L ( I > =H*B 


_M (J V = H *C 


MU 1 “MU* Mu ♦ I • 


_I fr (NOS t EO • 1 ) GO 

TO Ho 

MU1=MU1- 1 • 


Pl3-(SI*Sl )/2* 


P1»B *(Pl*SI*( 

S 1 0* 1 , J } 

rSlU-Sll 

+ 1.0 

pi^ sE *MU 1 


..P.U > = h*L AM * (pH 

+ P 1 + p 3 ) 

S-I e t NN - { , ) « HH 


I*I ’ 


I F ( I *GT . 1 ) GOTO 

6 1 


E»eO^K ( z\) /z* 

A = AO + L ( ^ n /2 * . 
B=bO+H{ z I ) /2. 

C = CO'*-P ( Z I ) /2* 

s I sS I • 

GnjO^I 

1 F { I *GT •*+ ) GOTO 7 
E;bK t 3 ) +E0 
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A u A 0 + L ( 3 I 

B s B 0 M ( 3 > 

c'sco'^pi'ii ^ 

St =SI+H 

GOTO 5 ■■ • • ■'■■"■ 

si«nn*hh_ 

Ev/»EV^ ( K { n +K < 'O +2 .0. ( K ( 2 ) 4.K ( 3 ) ) ) /a • 

EVl =*EV1>(L ( n+L ) + 2 .0* (L (2 ) ( 3 ) ) j / A. 

E V 2 = E V 2 ♦ <’m ( I ) ♦ M ( H I + 2 . * ( M ( 2 } ♦ M ( 3 ) ) ) / 6 , 

EV3*EV3^-(P(1 P < 2 . • ^ f 2 ) *pi 3 >J | /6 . 

P I 3-5 1 os I /2 • 

Plc((Sl0>i,)*Si+Pi}*rv2 

p3= CS lO-S I + I • > *EV 1 - - --- - 

ptjaMU I *E V 

"e V M t P 1 P 3 + P VJ o LA M " "" " ^ 

E==EV 

'a-evT ■ ”■ ~ ~ 

B=>EV2 ___ _ _____ _ 

C = EV3 

E 0 = E _ 

A 0* A 

Bn = B _ 



RUNG A kUTTA finished 
N s N '♦' 1 

I F ( D.EQ *S >gOTO 70 
EG3(NI=EV 

G0TO7i ___ 

"EgV<N)=:'Ev''”’ ” 

1 F ( N*LT • I NTEH J goto H 

' I F ( o'. EQ . M ) goto 9 

EmoM3=£V 

rEset for D_*'-t Integration 

OaH _ 

cast d "1 H T. "c'. 

. A 0 = I • _ 

EOaO.' 




B ^ sO * 

CO-0 • 
A * I 

F=~o • 

B = 0 * 

c=c • 

E V = 0 . 



EV2«0. 
Ey3 = C • 
E y M = 0 • 






E V 1 = AO 
Na 1 

GOTO M 
EmoM 9 =EV 

ALF*EM0M 3/EM0Ht| 

'do *7 2 Loa 1 ♦ 1 □ r 
LLb 102 -L» 
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M " BRPTaI)=£G3(LB)-ALF.YGV(LB) 

__ ■ continue 

I ... , ^ ^ ^ 

DO 7 3 LCa 1 * I 0 1 : 

BPT (LC ) »8BPT I Lc ) /B eT 

1 7 3 CONTjNUE_ 

SMcO*ODq 

DO 216 1=2, l_Cj , 

SHsSH^YpT ( I )+rPT ( i-i I )/2**( (FLOAT ( 1 J-l •BIohh^S 

1 2 U CONTI NUE 

A02*SM 

SM = 0 • oog 

I on 218 I =2 , I 0 1 

SiJ = s N B FXU 1 + a^XA IjrYiAaeJijjij *_H L 

2 I ft continue 

AHI=SM _ 

C0R=AM2»AH2/AMi 

__ SqoT = AM2/AM1 

''''RITE<6,7h) AMitAM2,5NOT,COR 

« - -7*4 F OH ^ fT ( V ,M 1 0 !>• A ».3.\ » *_.M25.__!_i 0 IS .„6.^3 A M2 /Ml 

l3x,« M2**2/M1= *iDlfi*6l 

IF(nPP 1 I iSl I l8l . 182 .. 

I 1R20q'jtINUE 

Write I 6, 76) ( BPT< 1 ) , I , I = I , 10 I l2 ) ... 

Go TO 75 

f l P 1 _Pu.NCH_ 7 7 , ( BpT (Jl »,I.= I 1. 1 0 I I 2 ) 

5 continue 

7 6 Format f • bpt = ..’ . o 1 2 .6 . 3 X . • i = . • 1 13) ^ 

7 7 F 0 1-V r } A T < G 1 I • 5 ) 

II End 


NO D I AGNqS T I CS • 
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Case ■ e 


M 1 “ 

•?HG9H6^c00 


M2=“ • 28 287*1-00 r 

BPT = 

• 88r 1 7«-0 1 S 

1 3 

1 

bPT = 

♦ 1 1030 1 "00 1 

f s 

3 

bPT = 

. • 9 3nM9b"00 1 

1 V 

5 

bPT = 

•920902-001 

r « 

7 

BPT=^ 


1 3 

0 

bPT = 

•22b733*000 

I 3 

n 

bPT = 

• 3 0 7 B M 2 •*’ 0 0 0 

T ^ 

J 3, ^ 

BPT^' 

•3B2Q53-*“o00 

1 ° 

15 " 

bPT = 

*N6n03B^n00' 

? * 

17 

BPT = 

• 5 3 3 7 0 l 0 0 0 

I - 


BPT = 

• 6On2M6 + n0O 

1 ^ 

2 1 

BPT = 

• 657 1 90‘^‘OOD 

I = 

23 

bPT = 

•70241 7^000 

f » 

25 

6PT = 

• 73^203*000 

I ^ 

27 

BPT = 

•75 1 2H6*00Q 

I — 

29 

6PT = 

•752680*000 

1 = 

3 t 

bPT = 

• 73f}0B9*0 00 

! - 

33 

BPT = 

• 7Q750S*00D 

I = 

35 

BPT = 

•66 1 400*000 

I * 

37 

BPT = 

• 600672*000 

1 S 

39 

bPT = 

•52a613*000 

I ^ 

4 1 

bPT- 

•44P877+0Q0 

1 = 

4 3 

BPT = 

•345438+000 

1 = 

45 

bPT = 

*242535*000 

I = 

H7 

bPT = 

• 1 3 <4 6 2 1 + P 0 0 

1 - 

49 

BPT = 

•243010-001 

t ^ 

5 1 

BPT== 

-•857352-0 0 1 

A =. 

.,.5 3... 

BPT = 

-•192775*000 

1 = 

55 

BPT = 

- • 29ti 150*000 

I = 

5 7 

bPT = 

-•387303*000 

I = 

5P 

HPT- 

- *469844*0 00 

I = 

. '6 .1 , 

BPT = 

-•539613*000 

I 9 

63 

BPT “ 

-•594726*000 

. I f 

„ 6 5 

BPT = 

-•633625*000 

I ^ 

67 

BPT = 

-•655103*000 

1 ^ 

6 9_ 

BPT = 

- *658338*000 


7 1 

BPT = 

- *642905*000 

I = 

7 3 

bPT = 

-•604779*000 

I ^ 

75 

bPT~ 

- • 5 5 3 2 7 * 0 0 0 

I = 

77 

BPT = 

- *486292*000 

I = 

7? 

oPT = 

-•399758*000 

1 ~ 

8l_ ,, 

BPT = 

- • 29r 1 1 4*000 

I ° 

83 

bPT = 

-•182998*000 

1 = 

85 

BPT = 

-*562353-001 

I 3 

0 7 

BPT = 

• 802290-00 1 

1 * 

89 

’ bPT = 

• 224 4 1 1 *000 

T ^ 

9 1 

BPT = 

•374373*000 

I - 

9 3 

bPT = 

•52p3 1 2*000 

t = 

95 

BPT“ 

•684655*000 

I 3 

..9 7 

BPT = 

•84? 1 56*000 

I = 

9 9 

bPT = 

• 1 0n000*00 1 

r 3 

lo 1 



• ■ 






• i \ 3A28^000 


M2^*2/M 1 


•321H2H-0D2 



PROGRAId LISTING 


AND 

SAMPLE OUTPUT 


PARAM 



r |T c; f OF 3 u B 0 r ! -V E A iV A ^ 

1 N T (', a r p L ' j I' s 
t< p ,1 \ L A M 

0 0 U B L E RE C 1 5 I -i ■'lU ( 3 j , S i 0 

C 0 ' • H ^ 0 '"J E / i_ >v !-l I 3 I 0 

C 0 V M 0 ■ ; / T ■- O / '^ J 
C 0 0 i'i / F t V e / -"J A J ■( 2 » L ^ » A U ^ 1 

C 0 M 0 ''• / N I isj E / "! ;) ^ ‘C 

C. o ^' 0 N / G R A E / ^ L ^ r S 

R E A 0 1 '5 » 2 3 I ) -1 J ^ E » '- A j S 1 . P L a T s 

2 0 1 F ,"i & [-'■ A 1 1 t B , F { 2 «> ^ t G I 2 • 6 , A ^ 1 

r A 0 ^ G , 2 0 2 I ( .i J ( K ) , K = 1 * '-1 0 0 F ) 

2(12 F 0 -< ->! A T ( 3 0 I 2 , 6 ) 

CaU. RARaM 

Stop 

t N 0 


HI 


SuiA F 0 U i 1 "-i e p A R A 0 

S.I — OMIINF ^M«A1 is iJfslGNKn 10 COMPUTE AND OUTPuT PARflMrTERS 

‘ ' fii,o Thkik sums Throughput thC 

ILL PtOT OUT the P,RaT THREE 

U E 3 I E n • C A L t F 0 


:^?i ^'‘,2 /’111 k2sQ'^ A^E ii/ M I 

THEF]R3T TbREEHliDES* al 50 rAcriA 

THRii main JN5T AS CaSEM 


0 0 F S IF 
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1 $ cfvf^Rt'D In Common 

C 0 M 0 N / 0 E / L A M » S I 0 
C 0 M M 0 ''i / T 0 / -M u 
r n f t M 0 N / (- 1 V E / N A f i X 2 , 1 ♦ A u X 1 

C n • i M 0 N / N I a; E / M 0 D r 
C 0 M. 0 N / G R A K / P L 0 r S 

\< (T |_ N N , A 

Uj T e p c R I , 0 , !'J . Z , I'^I I N T t T M T E R , P L 0 T s ' ‘ 

D 1 >i r N S 1 y j 0 J T P U T ( 2 3 ) , / M U ( 3 ) 

0 0 !,j B L E P R i- C I 5 I 3 .vj P | ^ ) , fr ( R ) , M ( M ) , L ( H t , p 3 ( 1 0 I I , E 3 I P < 1 0 1 ) ] E 3 2 p' ( 1 0 1 ) , 

1 E M O !‘i 3 f E M 0 M , A L P A » B L T A , M H X 3 t 1 Q ) ) , M M X M ( j G ! ) » p T’ T { I T I » R ) » 

1 £ 3 3 P < 1 J ! ) » E 3 H ; M t G I ) , £ 4 { I r i 1 ) , £ h 1 p ( I 0 i" I i £ M ? P I I f| I ) t E R 3 p J Jr,! } > 

I E B P ( I 0 I ) , A , -3 , C , E , A n , h Cl ® C 0 , £ 0 , fA u { * M i ; ( 3 ) » S I 0 , 3 I 

•Double phecisiOn ah i ( m ) » am^! m ) , sm - . .. 


N A i . X 2 I s A Plot c q n t o i_ p a r a m e t e r-J 

N 0 3 =: 0 E 0 meridional C a B E / N i j S - I for equatorial C a S e 


N A ! ; X 2 = [ G I 

w R r T E ( 6 , 3 H 7 ) 0 T s 

EOP'* AT ( 1 H } , « P aR am 
N 0 P T = 1 
NL'S = 0 
M 0 0 F S = 3 
N I ^ T = i 0 0 
I N T E R = N I M r + 1 
D 0 3 R M 0 0 L = j , M Q {> E S 
S T = 0 0 
N = 1 

Vi = 1 , / eVq a T I f i [ ! j T 1 


P L 0 T = * » A A / / ) 


Dy 3 i I = 1 , H 
K ( t ) = 0 * 

E ( I ) =0. “ 

ii ( I ) = Q . 

R ( I ) ==0 • . . . . 

D 0 I I = 1 , I 0 I 
E ii p ( T ) s Q , 
r-M P ( T ) =0 , 

E 33 P n ) = u ~ 

£ H 3 f I > = ij • 

£ 3 :> p ( 1 ) - n , 

E M R P ( I ) = 0 , 

£3 1 P M ) =:;3 c 
£ M 1 P { I ) = 0 . 

Eh M ) = 0 • 

£31 1 ) = 0 « 

T H 1 S S E C T 1 ■) !-! computes T h E E I R S T H 0 0 E s hLA P E AND T h E U T I I E m 0 n E 5 H A P E 
PARa'-ETEKS Ml and m 2 FnR CASE n 


U s 3 

I F in • F 0 . H 0 r '3 

E 0 5! 0 » 

B 0 0 “ 



lFt^.'OPr.GT*Q) '^'0 TO I 2 
1 1 ? P ( 1 ) - I . 
ao=.-i-:32P( I }■ 

oQ TO I J 

I ? £3 ( ? ) = 1 • 

L Q = f' 3 t I ) 

1 3 c 0 = 0 • 

A 0 = {? * 

C:>0 TO 3 ■ 

2 A Q =; 0 • 

0 = 0 * 

I F ( N 0 P 1 . G r • 3 J c. 0 TO 1 H 
li S 3 f> ( I ) = 1 4 
C 0 = F ^ 3 I ) 

0 T 0 IS 

I ^ £ 4 1 P ( j ) = 1 • 

Ao=f:H ) 

1 s B 0 ^ n • 

£0 = 0 • 

3, A - A Q 

B = p 0 
C = CO 
£ - ( ■; 0 
iJs } 

H 1 = 1 

N f\J i: [m 

S t = ( r,j f'.; - I , ) * H 

5 K ( J ) =H * A 
L { 1 ) = H M3 
M ( 1 ) = H » C 

M !J i 3 1 • * M U ( M 0 D £ • • M 0 ( Si 0 -i F ) 

I F ( fv 0 S • F 0 • I ) G 0 TO 4 0 

M U 1 = M U 1 " I . 

1 r pv 0 P T . G I * 3 ) P ( I ) = f t - S 1 « 5 I + 2 . “ S M ( 1 . ♦ S I 0 ) ) / 2 , » a 

1 + ( 1 , “ S I + S 1 0 ) « A ■kM 0 F J I . A M * H 
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CHAPTER 3 

Application, to Some Problems of Satellite Dynaipics 


The present chapter considers the use that can be made of the 
results of the previous chapter in some problems of interest in 
satellite dynamics^ A first field of application is in studying 

the nutational divergence of a satellite equipped with flexible appen 
aages, but this is the topic of Chapter 4 and 5. We shall be con- 
sidering here some other problems , such as the simulation of free 
oscillations, thermal flutter and variation of the spin rate due to 
the booms motion. 

3.1 Simulation of free oscillations 


3.1.1 Generalities 

In a motion of type E (equatorial) or M (meridional), the free 
oscillations can be simulated in the following manner. Given N modes, 
•^.(^)> C = 1,...N, with associated frequencies w . , and given an ini- 

j ^ 1 1 . 

. . dim^ 

tial distribution of displacements and velocities (t - j 

3 




I'n 
'C t Jt - 0 


I,. 


the displacement is written as a sum of modes 
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1 fi) (c, ‘“j V 
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Then 


'*'■-14 I 






r 'll 

y hfO Jf 




As an example. Figures 3.1 and 3./ are meant to illustrate that 
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starting with an initial shape identical to the first mode at A = 10 ^ 

= 0.1 , with no initial velocities, the stationary wave which exists 
in this case cannot be maintained if A is changed to 100. Not only 
has the frequency changed appreciably (t , is the period of the first 
mode oscillation for A = 10, = 0.1), but the second mode is present 

to an appreciable extent. 

3.1.2 Application to Satellite UK-4 

From data received through NASA GSFC on satellite UK-4, we com- 
puted the eigenfrequencies and modal shapes for satellite UK-4. This 
satellite has the following physical characteristics: 


UK- 4 Computations 


P 


" 30 rpm; 15 rpra; 6 rpm. 

« 0,00058 lb mass /in 

*= % ~5 4 X 102 X .45359 kg/m 
= 1.036 xl0~^ kg/m 


i 
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2 —4 
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(3.141592)^ X 1,036 x 


-2 

10 


X 

-2 


(7.01)"* 

2,869 

X 7.01^*/ 


.869 


86.06 
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ISrpra 


?^^^30rpm 


86.06 


« 21.515 


'6rpra 


86.06 

25 


3.44 


Data for prof.rams: 


» 1; 3.44; 10; 16.8; 21.515; 50; 86.06; 100 

« .042 

determined from 


[(-- - 1)C“- - - 1 

^ • V • 


'xh 


giving 


zh 


1.0767 


Results (see graphs) 

Graph 1; Resonance on thermal flutter at /X = • 4 or X = 16 . , i . e .• at spin rat' 

to “ ^ * 1.35 rad/sec = 12.9 rpm 

s l,rot 


Graph 2 : Mode shapes 

(0 =6;15;30 r.p.m. 
s 
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SATELLITE UK4 r 

ATTITUDE STABILITY 

- 


• 


Table I; 

Case M 






w(rpm) 

6 

10.2 13.25 

15 

22.8 

30 

32.3 

X 

3.44 

10.0 16.8 

21.515 

50 

86.06 

100.0' 

A 

1.84 

3.162 4.1 

4.64 

7.07 

9.28 

10.0 


4.08 

4.98 5.76 

6.24 

8.56 

10.75 

11.48 


1,38 

• 1.685 1.945 

2.11 

2.895 

3.63 

3.88 


Table IX: Case 

w(rpm) 

E 

3.23 

6,0 ‘ 

10.2 

13.25 

15 

22.8 

30 

32.3 

X 

1.0 

3.44 

10.0 

16.0 

21,515 

50 

86.06 

100.0 

A 

1.0 

1,84 

3.162 

4.1 

4 . 64 

7.07 

9.28 

10.0 


3.55 

3.64 

3.85 

4.05 

CO 

rH 

4.82 

5.44 

5.64 


1.2 

1.23 

1.30 

1.37 

1.415 

1,63 

1.84 

1 . 905 

« /8.715 « 2.96 

Resonance on Themal Flutter at 

u) »= 0) “4. 

s n 

- 4 

0/2.96 

.00, X - 16.0 

= 1.35 (Hz) 

== 12.90 rpm. 





K > 1 

P 


No posigrade resonance 
No nutaCional instability 


Var. of spin rate for 10% de 
.57%- 30 R1"M 

.785% 15 RPM 
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Fig. 3.3 represents the first mode of vibration for the three 
values of the spin rate being contemplated. Centrifugal effects are 
noted as Etkin's number X is increased. 


3.2 Resonant thermal flutter 

3.2.1 Determination of resonant frequency 

[3-2] 

It has been shown by Etkins and Hughes that assuming a 

relatively simple model for the boom's thermal curvature (inde- 

pendent of due to the sun's heat input during the spinning motion, 
the steady-state oscillation of the booms would be described by 
'h s. - ^o) 

^ Q , • 

In order to find for which spin rate w the motion will diverge (have 

s 

an amplitude tending to infinity) , these authors solved equation 
for boundary conditions, 

E(0) = 0 E'(0) = 0 E"(0) = E"'(0) = 0 


and V iry X until very 'large values of ^(1) are observed. The analysis 


was limited to satellites of zero radius,. 

[3-1] 

An alternative approach was proposed , which is recalled 

here. If in Equation (2-2” S) » *^1 tend continuously to o) along 

the eigenfrequencies curves ^^(X, Co)- t’he spatial part of a solution 
to Equation (2-2“7 )> normalized to unity at the tip, satisfies b.c. 

In order to also admit boundary conditions 0,0 for the zeroth and first 
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derivatives at ^ 0 for the second the third derivatives at 

5 = 1, ^o) should be scaled up by an infinite factor, i.e,, the 

amplitude at the tip tends to infinity. Thus, resonance on thermal 
flutter will correspond to the intersection of the curve, for given 
Co 


(Oj f t 


iT 


I ('^) 


with the bisectrix of the first quadrant (Fig. 3.4) 


No thermal flutter resonance can occur for 

second or higher modes 
as is shown on Fig. 2 ♦ 2* 


3.2.2 Application to UK-4 

Using the above data for UK-4, the thermal flutter resonance 
point was found at (Fig. 3.5) 

X.llS.O 0.042) 

and for the physical characteristics of the satellite, this translates 
to 


(0 




^ 435'i^z = fZ. ^0 


a spin rate to be avoided for steady-state operation. 

3,3 Variation of the spin rate due to the free oscillations 
3.3.1 Method of calculation 

It is often of interest to satellite users to know what amount 
of spin rate variation can be expected, due to the vibrations of 
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the boom. The equatorial vibrations will cause a very slight varia- 
tion of the spin rate described by 

7 , T 

Li di 


where T is a torque due to the moment at the oot of t boom and to 
the she r force acting through the central hub radius. This is de- 
scribed in non-dimensional form by i 


dl*) 

i 'C 

f 


[ 



1 



1 



TT 

- k 

^ . 

ii.c 






p 

or, after integration, and with T -jr — 

hub 




L 




f % 




fY 


( * TjLF 

J r\ 


(3.3-1) 


in which x is the value of x maximizing the integral. This value 
can be obtained using a program such as SIH, which is listed at the 
end of this chapter. 

3.3.2 Application to UK-4 

Using the above data for satellite UK-4, the maximum variation 
of the spin rate for an assumed 10% deflection of the boom was de- 


0.57% 

at CO =30 
s 

r.p.m. 

0.755% 

at (0 =15 

s 

r .p.m. 


termined to be 
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8*0 


4*0 

3, 6-/5 


0«0 


/ 


SfirtUl7t U * rt/l^:>AM£W 7 flL NATV^dL ff^EQacNcy ^s- 

SAr£LLlTE Spi ^ ^4T£ 


f 


/ 


/ 


M 

t*r 

\ 


/ 


Case £ 


/ 


/ 


I = o-04a 


0 


/ 


^ l?^sonaftce on ih« final flt^tUr 


/ 


/ 




/ 


/ 


/ 






tZ*0 


FIG* 3-5. RESONANCE ON THEIiMAL FLUTTER: UK-4, 
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PROGRAM SIM 
LISTING 



I 

I 


L 


$DUP E730MK I8t i^^INT*3C0 

“^■"FOR " S IM ■■ ” ' ~ 

SUBROUTINE SIM(NDS) 

D0UBXF”P RPC' I SI GN ' MU VX N Cr; Y NO 

REAL MU 

integer SAM 
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S IM 


IT- 

I 


THI<; PRir.RAM COMPUfES THE DISPLACEMENT OF THE FREELY 
—TTOT^TTNG“BOOH''Ar'SPECIFlED STATIONS, STARTING FROM A 
C DISPLACEMENT INITIAL CO NO I II ON 


VIBRATING 


GIVEN 


C 

T 




DfE 

(lOl 


32P( 101) 

)t 


t 


r 


t'OaTfOX/T) N £7 M U7L AM S I OV RG A M A" » P K , N U • X N 0 , Y M 0 
DOUBLE PRECISION C0R,MEl,ME2 c^dni 1 FRlPdO 

IE33PI 101) ,E34P( l.Ol) ,E4( lOl 1 .EAlPdOl ) ,E42 P(l0l ) ,E 

— 3E44P ( 10 I ) » At BfCfCfAOfBOfCOtEOt 

4S10,SI 

— 0 ox BTE“^P RE C'l S I ON EDA(2)fEIG(2) 

REAL HtNNfLAM - 

■“■INTEGER 1*0, NtZf NINT, INTER 

DOUBLE P REC I SI ON E M0M3 ( 2 ) U MOM^ ( 2 ) * ESHR3 ( 2 _) * FS HR^ ( 2 ) 


t 


DOUBLt C. -tcj.. , TPAPI 

IRM0M3( 2),RMnM4( 2 ) .RSHR3I2) ,RSHR4(2) •^PTUOl 
4-rriFnId2;2).'ALFA(2), ' M«X3t 101,2),MMX4( 101,2) ,BE 

DIMENSION AA(2 ,2) ,BBB(2) ,CCC(2 ) ,GAMA(2 ), CC2 (2 ) 

" INTeGER PPP( 2) ,QQO( 2) 

double precision RTORK ,DT0RK,£0A2P l*EDA3P 1 

"“REAL DLTTtMOOE 

GAM A ( 1 ) =RGA MA — 

'GAMAT2Y=0'.0 

ElG(l)=MU - 


* 

2 ) , 

TA( 2) ,MUl 



EiG(2) =0. 
EDA( 1) = XN0 
'E0A(2)=YN0 

lENO^l 

fYAXr = 2'7 

_ STP = 25. 

NR=2 

THIS PROGRAM 


REVERSED I NT EGRAT ION ONLY 


300 

TT" NDS 


USES 
WRITE16,300) 

FORMAT(*OS1MULATI ON_^EN T_£ R E 0 * ) 

■^DT RECTIUN'' SW ITCH 
= 1 IN PLANE, = 0 OUT OF PLANE 
IB1G=1 TO READ BIGMQDE 1 ONLY 


C NDS 

SET 
I VAR'l 

■ IB1G=1 

1FIRST=1 


KK = l 
N=l 

NlNT=iOO 


I 

I 


INTER=:NlNT-»-l 

ANFRg=SQrtT { LAM) 

A^J_1!IL^*60)_ ANFRQ 
FO'RMA'r (Th *mfrq='sqrTTaM^^^ 

S 1 = 0 . 


CLEAR ARRAYS 


DC 31 I=l»4 


K( I ) = 0 . ^ 

__L( I ) = 0 . 

M{I)=oV '■ ^ 

_ 

DO 1 1 = 1,101 

__E^ 4 P( I ) = 0 . 

E 44 l>Tl ) =d, 

_E 33 P( I ) = 0 . 

ewTi) = o. ■ 

_E 32 P(I)= 0 . 
e 42 P{I)^ 0 : " 

E 31 P( I ) = 0 . 

E 4 iPTl )=d. 

J^<I)= 0 . 

E 3 ( I y =:0 . 

‘ 

IF(D,EQ. 4 ) go to 2 

E 3 ( U=X, 

E 0 =E 3 (iy “ “ 

_co=o. 

Ad=or" ' — 

G 0 ^T 0__3 

A 0 =t 41 PU) ' ' ' " — 

^ 80 = 0 , 

E 0 = 0 .”“” ■■ 

A= AO 

C=C 0 

e=eqT ■ ■ ' - 

N=l_ . 

.NN = ^4 

s I = i'm- 1 . ) ^ ^ — 

L(n=H <'8 - - - - 

M( ly 

MUl= I . + E IG(KK )«EIG(KK ) ' 

1 F(N 0 S,E Q - 1 ) „ J^O _ T 0 4 0 

MUi=MUl“~l,' " ■" " 

P( r )=( (-Sl^Sl+ 2 . *SI ♦( l.+SIQ) )/ 2 .*B 

♦ 1 1 • “SI + S 1 0 ) * A + MU I ♦£ ) ^L AM*H 



=c 

N=N+l _ 

IF (iSl,LT* INTEK) GO TO' 4 ^ 

EMOM3( KK) =E'3{ [ NTER) _ _ 

ELSM3{ KK )=E 32P( 1) 

RMOM3(KK )=E32P( INTER ) 

ESRR3rKK)=E33P(l ” 

RSHR3(KK)=E33P( INTER) 

00 28 1=1 » 101 
KMX3U ,KK)=E3 ( I ) 

0=4 ; 

GOTO 2 


si=nn*h 

E4fZ ) = E4(N ) + ( K( 1) + 2.=>K( 2) + 2,*K{ 3) ^K(4) ) /6. 
E41P(Z ) = E41P (N) + (L( l) + 2.’^L(2) + 2.*L(3)^“L(4)) 
~E42P ( Z ) =E42 P (N)^(M(l)+2.^iV(2)+2*^M(3)^M{4)) 
E43P( Z) =£43P (N) + ( PI l)+2.*P{2)+2.*P(3)+P(4) ) 
E44pt^ )=LAM«( ( ( SIO+ U ) 1. 

l*E42P( Z)/2*+( I .-S I+S I0)#E41P (Z )+MUl^E4(Z ) ) 
E=64(N<^1) 


/6. 

/ 6 . 

/6. 

-SI + SI 










t" 


73 DO 213 I=l .2 

“7 13 ~W R T T E { 6 , 2 1 4 ) A L F A ( I ) f I 
214 FORMAT (IH , • AL F A= • » 0 1 2 . 6 # 3 X f * KK= • » I 3 ) 

ITTTBT GI EO. T ) " GO“T 0 ‘109 * 

DO 104 KK = l ,2 _ 

-WM= 5 I ' " ■ 

00 104 LL=1»2 



THESE COEFFICIENTS FIT THE BOOM T^0_ INITIAL^ 


COFOIS (LLf KK )=BPT(MM,KK ) 
AA(Lr,KK) =COFDIS{LL,KK) 
MM=MM+50 


17^ V ERT HE" COE F F I CT ENT 'M AT R IX 




rriTOD'A 

c 

109 


CALL GJRI AA, 2» 8BB * CCC , PPP t QOQ * S 1 1 3 ) 
W“l 6 r r=l ,NR' 

$UMM = 0. „ 

'oo“r6‘o'”j'=r;N"R 

SUMM=SUMM+AA n 1 J ) *EOA u ) 

CC2iI ) =SUMM 

WR I rei 6, 162)CC2( I ) * I 

FORMAT ( IH » *CC2= » , E 12.6, 3X, *1= •, I 3) 

riT'soTiJ rroN~F'OR rt"H er e 



FF^OV 

FEE^O.. 

kT0RK = 0.' 

0T0RK = 0. 

t“INT-0; 

0LTT=1 .0/ (STP’J'GAMAIl ) .*3 . 1416 

WITF( 6,50)DLTT 

FORMAT! IH * *OLTT= * , F 10. 6/) __ 

SFAC = l ./EIG( 1 ) 

V^RI TE ( 6,131 ) SF AC _ , 

TORM’A'TI IH V* THE' NO.*" OF ROT CYC PER FUND VI B 

DO 122 I = 1,MAXT _ _ 

f=OL TT* FLOAT ( 1“1 ) 

DO 129 KK=1 ,2 . 

TFAC (KK )=CC2{KK ) ♦COS! GAMA( KK I*T) 

EDA2P1=0. 

ElDA"3Pl=0. 

^DAE?=0. _ . 

EDAE3=0. 

DO 1L9 KK=l ,2 

S0M2 = ( RM0M3 ( KK) -ALF a ( KK ) *RM 0M4 { KK ) ) * r FAC ( KK 
SUM RSHR 3( KK )-ALF A ( KK ) 4^RSHR4( KK ) ) *= TF AC ( KK 
~SUM6^'(6LSM (KK )-ALFA(KK )*ELSM4(K< ) ) ^^TFAC(KK 
SUM 7 = I 6SHR3(KK)-ALFA(KK)^6SHK4(KK))=*'rFAC(KK 
■ EDAE2=EDAE2+SUM6 


EDAE3=EDAE3“SUMT 

“EDA2PI =EDA2Pl +SUM2 

119 EDA3P l=EDA3P 1-SUM3 

RTdRK = 6'DA2PT-S ■ ' ’ 

AVTK=|RrORK + DTORK)/2, _ _ _ 

dtORK=RTORK 

IFt 1 .NE. I ) GO rO 13A 

AVTK=d. 

134 T INT=TI NT-»-AVTK»DLTr/ANFRQ 

F E£= 360. 4c' SF A C L OA T ( I - 1 ) / S f P 

I F{ FEE-360. ) 140, 140, 117 

117 FF = FF+rr • “ ' ' ■ 

F£E=FEE-360. 

i 40 TF ( r. N E'. 1) ■ GO 'to 4 4 • 

WR1TE(6,215)T,FEE,FF 

215. F OR MA TTI H "7 • t I ME " I N' SE C •'7 f 9 :37 3 X 7* F E E= « V F9 
__ WRITE( 6t 130)EDA2PI,E0A62,EDA3P1,E0AE3 
130 FORMAT UH ,»ELAS RMOM • , D 1 2 , 6, 3X , » T IP 

^RSH£AR = » ,012.6,3X,»riP SHEAR*,El?.6/I 

WRltE(6, 4l)riNT 

1/S0RT(L^)4cI|^TEGRAUM*dlTT) = 

4 4 D O ~T n “Fl = I ,1 01,50 " ■ 

SUM =0. 

DO Ti 2 Kk=l,2 ■ ' “ 

MODE = BPT (LL,KK )4'TFAC(KK ) 

112 SUM=$UM+MODE ' " ^ 

1 F( I .EQ . 1) GO TO 1 14 

I F n V AR7 eq7 1 ) ‘GO Tto" 1 1 1 

114 WRITF(6,115) SUM,LL 

111 CONTINUE - - - 

U5 FORMATIIH , * ADA • , E 1 2. 6, 3X , • STATI ON ' , I 3) 
IF(IVAR.nFI) GO TO 122 
6,121) SUM,TI NT,T,I 

122 continue " ‘ 

121_ FORMAT (IH , * END I S P • , F 12 . 6, 3X , • T INT ' , E 12 .6, 

120 WRITE I 6, 568) 

56^ FORMAT( IH /) 

"GO TO 100 " ” “ ~ 

U3 SAM=2 

WRftET67ll6) SAM 


FORMAT! IH 
CONTTNUE 
GO TO 42 
RETURN ' 


DUMPED SAM=*,I3) 


END 
IN = 


316, CAROS our = 



jUL>t\U«JtANu Lr»Jt\ V M » IM I j » c P • fUfj'i 

’ INTtLo^K N r P ( N ) r Q ( M ) 

01 Mc.:>lS I Ol\J A ( Ni » ,\i ) f rj ( M ) r C ( • J ) 

CPS^i •»(lU*)+=*^~7 

101 cP5=dPS/lu. 

lF(£Pb*LT. 60 To 102 

00 6b K=lrN 

PIVOf-u. 

DO 5 I=K»isJ 
00 b J=K»N 

lF(A6S(A(If J) ) .L£,AOS(PIVOD ) GO TO b 
PI\/01=/a(I fJ) 

P(K)=1 

0(K)-J 

5 COi>iTli'JO£ 

iFtAdS(PlvOT) •i_E,EP5)G0 TO 101 
.. 1F(P(K)-K)b»10f6 

6 00 7 
L=P(\) 

2=A(L»J) 

A(L» J)~A(K» J) 

r AiKfO)=Z 

10 IF(0(K}-K)llrlb»ll 

11 UO 12 1=1 nM 
■ L=Oti\) 

2=A(i»L) 

A ( X M- ) =A { I r K ) 

12 A(ln\) =Z 
1!> DO 2 b o=l»N 

lF(u-;<.)2UrlGr20 
16 rt(JJ=i ./PIVOT 
c t o ) — - 1 • 

60 TO 22 

20 ■ B( J)=“rt(K rJ) /PIVOT 
C(J)=A(J»K) 

22 A(K»J>=0. 

A(J»i^)=0* * 

25 CONfXiJU^I 

00 OU l = lrlM 
DO bJ %J=liN 

30 A(i» J)=A(I#J)+C(1)^B( J) 

35 CONI l^^tjE 

00 50 M=l#iM 
K=N-y.+l 

IF IP (A i ~K )40r4b»40 
40 DO 43 i = lr.M 
L=P(j\) 

2=A(i^l) 

. A(If L)=A(IrK) 

A(XfA)=2 
43 CONliMUE 
4b lF(6tA)“j<)46»5J^4G 
46 DO 4b J=ifN 
L=DC^) 

Z=A(l_»j) 

A(Lr J)-A(K» J) 

A(Krj)^Z 
46 CONTI- 4j£ 

50 CONTi-vuE 
RErUKN' 

102 I'jRXTE^O'lub) 

103 F0r<>;/\ri2lHU GJ:< COULD NOT uO IT) 
RETURN 7 

FN') 



CHAPTER 4 


Simulation of the Motion of The Central Rigid 
Body and its Elastic Appendages 

4*1 Introduction 

In the previous chapters, the problem of determining the modal 
shapes and frequencies of the rotating structure was examined, and 
applications were studied in which these modes are utilized. 

In the present chapter, equations of motion are written for the 
generalized coordinates representing the flexible structure and 
for the angular rates of the central rigid body. A simulation of 
the spacecraft motion is then possible. Various cases of simulation 
are examined, and the effect of modal truncation and of nonlinear 
terms is discussed. 

4.2 Modal Equations of Motion: equatorial vibrations (Case "E’l, 

for equatorial) 

-> 

4.2.1 Constancy of H. 

In what follows, it is assumed that the motion of the center of 
mass of the spacecraft is negligible (or that only antisymmetric 
motions of the booms are considered) and that the "limited approach" 
is taken^^”^^, i.e. the motion of the spacecraft’s center of mass in 
inertial space can be determined independently of the attitude. 

If over the time of interest, i.e. a few teas of spin periods 
or so, the torque- impulse due to all environmental attitude per- 
turbing torques (gravity-gradient, solar pressure, magnetic, etc.) can 
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be considered as negligible, then very sensibly the moment of momentum 
H about the center of mass remains constant: 

H = H(o) = constant vector (4.2-1) 

in which H is the value of H at t = 0. 


4.2,2 Representation of the elastic appendages 

Consider a particle of a boom, having non-dimensional abscissa 
located along axis +x in its undeflected position. Its elastic dis- 
placement, n = , is represented in terms of the modes $ , (? = 



(4.2-2) 


in which the q. are non-dimensional amplitudes, dependent on the non- 
— 2tt 

dimensional time t = to t, with o) = — the angular spin rate of the 

satellite in its nominal motion. N is some positive integer, which 

specifies the number of terms after which the series is truncated. 

We recall that the 'J’.(^) are orthogonal modes, normalized to 

J 

unit deflection at the boom^s tip, so that 




boom 


m 


l,j def 


j ^AOdK > 0 

boom 


j k 
j = k 


(4.2-3) 

(4.2-4) 


= ^ + 

both m. m . are assumed to be known quantities, determined as in 

1 >3 2 , J 

Chapter 2. 





4-3 


4«2.3 Kinetic energy contained in the elastic structure 

The total kinetic energy, T, is made out of two parts: one 

is independent of the generalized coordinates and the other one, T^, 
depends on the q. and appears as the integral of a density » More 

specifically (Fig. 4.1) 


->2 

V 

^ all particles 2 


( 4 . 2 - 6 ) 


V = 0 ) A(r + ^ ) + 6 

m m, 0 


(4.2-7) 


in which w is the instantaneous rotation, r^ ^ is the vector coordinate 
to m in its reference position and ^ is the elastic displacement from 


Computing T, 


^ 2 ntvQull / 

A ) . 5 f 2 (J A 6 j. S f 0 j 


Now, for small linear displacements of the elastic parts 


->■ -> 
6 = w 1 


-W(0 

^ A t Z 

o)A 0 = Q 


Z I i V. ) 

Hl fcirt57(c 2 2 


(4.2-9) 
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If w , 10 ,W are assumed, to be of first order of smallness, (4.2—9) 
X y 

is rewritten 


f) ^ 


2 


(4.2-10) 


in which w is the (constant) nominal value of the spin rate, 
s 


00 



0 ' 

-00 Xi 

I y J 


Then 


j_ s 

tlfiSTfc 

Furthermore, 


)AZ kl)^ 0 -W xl .0(£^J(4.2-11) 


5 A. fAfirs 


i ifVl 0 


2 . 1 ^ w 

2 


(4,2-12) 


and 


J. 5 I. 

2 AU ^ 

£i.fA?(?r 5 


r 2. rz. W 


Finally, 


(4.2-13) 


1 1 A S ),l 

Z Ml ^ 


=: 0 


(4.2-14) 
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Introducing expressions (4.8-10) through (4.2-14) in Equations 


(4.2-8), we obtain 


T = T, + - 


2 2 


^ds(w 0)^ + +2 (jj^Xj^W) + O(e^) (4.2-15) 


Since the element of curvilinear abscissa, ds, is related to dx by 
ds^ = dx^ + dy^ = dx^(l + 


or 


ds = dx(l + 


and 


dx = 


dx(l + +. . .) 


A /I l/9w. 2 
ds(l - 2 ( 3 ^) • • •) 


Therefore, consistent with the order of magnitudes retained explicitly, 
(43-15) can be rewritten with x instead of s as the integration variable, 

T = To + i [ P dx(w^co + (ji^ + 2 0 ) X vv) + O(e^) 

2 j S z 


The "flexible body" part of Tp , however, has to include a correction 
term, since for terms involving 


f£ 1 

9 

upper 1 

X ds = 

x^ dx + 2 

J 



(I 


“PP^’^2,3w,2 

dx 
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i.e« with an integrand of zero-th order of magnitude, we can write 

^ I ^ IbUsStrL P \ 

0 . iL N 


(IW (4.2. 


17) 


Neglecting terms of 3^^. order of smallness, 

^ o 2 ^ 


and using (4. 2-17) , 




O 


ii 


(4.,2-/5) 


Now, if the integrand is 


we obtain 
2 


Xj = (Xo + x)^ = H- 2XoX + Xp 


Til 

i/rt >/o 


J X J X oU + i f x C^i: 

O o 2 '^0 

or neglecting terms of third order of smallness 

1^' 1 Ji -. Jt - ( ^ )‘ J. 

■ {' - Jj/ !»-*)( 


7L 


Finally, 

■i , . ft 


1 

“ ^ ‘'o / 
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Therefore, T© is rewritten, with the inertia dyadic of the 

rigidified, undeflected total reference body as 



( 4 . 2 > 20 ) 


Collecting (4.2-16) and(4.2-20) 



w' -2 

+ O(t^) ( 4 . 2 - 21 ). 


4.2.4 Potential energy of the elastic structure 

For pure flexure in the (x,y) plane, the potential energy is 
given by 



rl w 2 

( — 2^ + O(e^) 

9x 


where it is legitimate to use x, instead of s, as the integration 
variable, to the order of the terms explicitly retained. 

4.2.5 Equations of motion for the elastic modes, equatorial vibrations. 

4. 2.5.1 Equation for the jth coordinate, q 

3 

At this yoint we introduce the modal representation (4.1—2). For 
the sake of simplicity, let the bars on t = [ry- ], be dropped. Further- 

-L/Wg 

more, let the energies be non-dimensional ized by 

and the lengths by I 

Although ( 0 ^, non-dimensional value of the nominal satellite spin-rate, 
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is 1, we shall for clarity retain it in the equations. 

In non-dimensional form, with ’ designating derivatives with 
respect to t, and = ? + ?o> 


T i ^ (T ^ 1 oT. j - — f ^ 


L_L_( 

— 2 

with X = ~ — m , as in Chapter 2. 
El s 


Now let, with t the non-dimensional time, 

n = 2 

J-i '■* <» 


Then 


= A?' 

= I Z 'I ^:K 

u> t. “j ^ 


•2 

n 

= Z 2 

k=' 


N 

9n 

= 2. 1 

9:^ 

J=' ' 


n , . 

r ^ o. a 

ri t-/ 'K 


u 

1 

N' 

Z 

£1 ■ i 

if; 



kzi 


^ c 


n 

1 

Z 

k=i 

1j1. 

i 

/fk 
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Thus, the Lagrangian function is 


U M 


2 ‘2. -Lj’’ k=i 1 - Vk <t ^ (j Ik J y s 


Now, for any i = 1, 2..., 


d ^91!^ V 9o^ _ Q 

dt » 

3 qi aq^ 


(4.2-2J) 


We recall that, as in (4.2-3), (4.2-4), (4.2-5) 


j = 0 (j ?^=k) 

boom 


d). E m_ . 

j l,j 

boom 


boom 


the <f>. having been previously normalized to unit deflection at the 
boom^s tip. Using these relations, and (4.2-25) for i = j, after de- 
f ining 


^jk ^kj def 2 


d3>. d4> 

[(1-C2) + 25<,(1-C)]^^ d? (4.2-26) 


boom 


^jk ^kj def 


d^(f) d^<t> 

2 .^2 


"^boom 


(4.2-27) 


we obtain, for the j th modal coordinate 


ofa, : 


I 


-2 

• id, 4 


ck ’’J ' '<( 


« A e 

+ ^4;.a +il 1. <] iJ. 
tn k r K. I« ^,L ■ 


(4.2-2S) 



4-10 


4. 2.5*2 Evaluation of the coefficient of q. 

J 


We nov/ evaluate the coefficient of q. in (4.2-2g), say c , 

J J 


N’ ^ 

c = 1 Z X 1 + a.i, - OTV . w; 

Tv krj kr| ^ 


i T -- 0 

^ <1, 1, — '11V, . W 


(4.2-2i/) 


From Equation (2.2-23), written in terms of 


1 

% 

V 

o 

II 

+ 


\ 





Thus (4.2-25) takes the simple form 

•• . -2 

m- ,q, + m, ,a>.q. = oj 

1,3 J l,j 3 3 z 


or 


m. 


“ ^ -2 2,1 - 

q. + u). q. ^ (jo 

j 3 3 z 


(4.2-JO) 


A few remarks should be made regarding (4.2-30) First of all, 
the modal equations for the coordinate reduce to a harmonic 
motion, in the case where to is constant. Second, as has been seen 
in Chapter 3, the "driving amplitude", measured by the non-dimension- 
al ratio 


m 


m 


lA 


is strongly a function of ^o, 3-nd to a lesser extent of X, Etkin's 
number. Thirdly, it should be noted that it is only because, for 
the sake of consistency, the difference of an integral in s (curvi- 
linear abscissa) and x was carefully considered when the integrand 
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was a quantity of zerot^"^ order (as detailed above) , that term 

- m- . 0 )^ q. 

1,J s 

coulti. finally be cancelled in Equation (4.2-30). Failure to make 
this distinction leads to having this extra term still present in 
the final equation and in order to "fall back” on (4. 2 “10), one 
has to introduce, rather belatedly, an additional term due to a 
"rotational potential”^ \ Finally, if linear distributed 
damping is introduced, Equation (4.2-30) takes the form 

q +2 P d i- io a. ^ ^ (4.2-3/) 

whose derivatives are taken with respect to non-dimensional time. 


4.3 Modal Equations of Motion: meridional vibrations (Case "M" , 

for meridional) 

Without repeating in the same detail the explanations of Sec- 
tion 4.2, we now derive the modal equations in the case of motions 
parallel to axis-z (meridional vibrations) . Only the relevant 
differences are underlined. 

4,3.1 Constancy of H 

In the absence of attitude perturbing torques, the torque-free 
motion has the integral 

H = Ho (4.3-1) 

where Hq is the value of the moment of momentum at t = 0. 


4.3.2 Representation of the elastic appendages 

The displacement n(x) = of an element of boom located at 

^ = - of axis +x in terms of the modes ^.(?) for meridional motions Js 
Hf 3 


x: f.. I.''-* 


(4.3-2) 



4-12 


Again q, are non-dimensional amplitudes, function^ of the inon- 

dimensional time t = V7 — • N is positive integer specifying 

1/ w 

s 

the number of terms after which the series is truncated. 

The "meridional" inodes are orthogonal 




and we have defined 


m, . - j, <l>.(^)dC > 0 

l,j ^ boom j 


(4.3-3) 


(4.3-4) 


, e,<l>.(Od5 with e = K + So (4.3-5) 

boom 1 j I 

These quantities are known as functions of X, Etkin’s number, and 

Co = , hub non-dimensional radius. 

X/ 

4.3.3 Kinetic energy contained in the elastic structure 
With the same notations as in Section 4.2.3, 

T = To (rigid part + flexible part) + 



Now 

? = w t 

z 


(U 


A? 


(0 w 

y 

-(A) W 

z 

. 0 . 


2 i A?)^- X I (^-3-6) 

2 2 J 


and 
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Let X. = X + Xo . Then 


0 ) A r 


m, 0 


0) X- 

z 1 


-w X- 

y 1. 


The next terms are 


J.5 I 2 6j^ 4)^ a^,i<r) + (4,3-7) 


2 /Ju 
£1. 


Z 

AU £z. 2 


— w C ~ 


X or\ W 
2 


(4.3-8) 




2 (Hi 




_ J 2 9 (Uy W 

2 J 


(4.3-9) 


1 1 Z'ty^ (oJ /\ i ) . I 

z 

£i. 


^ 0 


(4.3-10) 


Introducing expressions (4.3-6) through (4.3-10) into (4.2-8), and 
since we can substitute dx for ds when the integrand is of first order 


of smallness, or smaller. 


I . 'MiiSfi i 
2 2 




1 r /i 7 2_,,2' 


dx(w^ - 2w w x,v7 - 2m x^w - ^[(^ -x^) 
^ s x 1 V 2 


+2xo ()t-x) ] (~) ) + O(e^) 


(4.3-12) 
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4.3.4 Potential energy of the elastic structure 

For pure flexure in the (x,z) plane, the potential energy is 


V 

2 


+ 0(c3) 


9x' 


Again, it is legitimate to use x, instead of s, as the integration, 
to the order of the terms explicitly retained. 

4.3.5 Equations of motion: elastic modes, meridional vibrations 

5. 3.5.1 Equation for the jth coordinate, q^ 

The kinetic energy, T, and potential energy, V, are non-dimen- 
sionalized by the quantity Note that, although non- 

dimensional value of the nominal satellite spin rate, is 1, it is 
retained as '* in the equations. Let ~ 2 . 


Te — -- (y‘^20^ u, {y-2^iy - 


V 


/ _ 




---- ^ f 

9 -k ''L 


< t («| 2 a. y ^ 

With the same substitutions as in (4.2), x^^e obtain the Lagrangian 


X 




^ r 1 { y 


1 1 i ii -2<-i ‘j, ( <!, t -zij I f f 

V M ‘ ' y j j ‘h i 


2 j-., k., ^ > 'k M " /■'''<»<» '"j y -</ 

- = t(.-s)! U1, ^ U% SJt. I ■‘i 


Now define 


^1 

^jk \j def 2 J 


boom 


d<}> . d(|>, 

[a-^y + 2Co(i-e)]5^^ de 
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^jk ^kj def J 


d^(f) d^<l>. 

boom d? 


For modal coordinate , the Lagrangian equation is, with Equations 
(4.3-3) through (4.3-5), 


N i _ _ 

m . q. + = b.. + 1,^1 "■ (4.3-14) 

1,3 3 X k=l jk k=l jk 2,3 7 s x 


4. 3.5. 2 Evaluation of the coefficient of q^ 

From Equations } , Section 2.23, we obtain 


= b.. + a.- = 0 
X 3k 3k 


3^k 


7 b . . + a , . = w f m. . . j =k 
1 33 33 3 1,3 


Thus Equation (4.3-14) can be rewritten in the form 


_ ^9 • A _ _ 

q, + 0 )^ q. = -• (u) - oj (jo ) (4.3-15) 

3 j 3 n»i i y s X 

^ > J 


_ th 

in which o). is the j eigenfrequency of meridional vibrations, a 


function of X and 


Again, if linear distributed damping is introduced, the equation 


of motion becomes 


^ 9 • — i — 

a + 2vw Q + q = (w — co oj ) (4.3-16) 

3 3 3 3 3 ^l,j ^ ^ ^ 

m . 

—=^ thus appears as a "driving amplitude". As seen in Chapter 2, it 

m, . 

1,3 

is also strongly dependent on Co, 3,nd to a lesser extent on Etkin’s 


number X . 
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4,4 Equations for the rates; equatorial vibrations (Case "E") 

The equation for the time derivatives of the rates are now de- 
rived from the constancy of the moment of momentum for the torque- 
free motion, as given in (4,2-1). 

Since, about the center of mass, 


H = H(o) = 


r A r dm 


H = r A r dm 


Computing, with the same notations as in 4.2 and 4.3, 
m, 0 


r=6 + 0 ) A(r 


^ 


. . -t, 


r = + 2(ii A 6 -h (jj A r .-f-(jjA6+{i)A((oAr ^ ^ ^ ' 

m,0 s m,0 


Let H be divided between a part "relating to r , and a 

Til 5 

II 


relating to 5 , H 




and 


“ Jboom 


->■ . r ^ \ >1 1 

r A{o) A r „ + m A(m A r . ) } ds 


m, 0 


m. O' 


II 


boom 


[r A {6 -i2o) A6 + (c A6 + {joA(o)A?)} 


(4 . 4-1) 
(4,4-2) 


6 A (mA r ) + 6 (mA((oA r ))]ds 
m, 0 0 

Note that as has been seen in Section (4.3) and (4.3) 
ds = dx(l + ^ higher order terms 

in which the elastic displacement (along -Hy) is 


(4,4-3) 


6 = w(x)T 


part 


(4.4-4) 
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Therefore, if o) , (o (normalized to (n , nominal value of the satellite 
X y s 

spin rate) and n(x) = are considered to be of first order of 

smallness (0(e)), the equations for the rates deduced from (4.4-1) 
should be written with 


[ . , . ]ds = 


[‘•‘]dx , 


boom*" ■* "'boom 

for integrands of zeroth order of smallness, or smaller, if only 
quantities of first-order of smallness, or larger, are retained . 
Thus, neglecting terms of order 3 of smallness, or smaller (with 


[—Ids =1, [.»«]dx for an integrand of first order of smallness, 

boom ^ boom 

or smaller ) , 


~ W d Tc. 

rft c/ 


4i>, -J, 


CO w- W (0 

Jw J 


(Ojg, w) d.X 


(H ) = - f W (Lx 

vovwn 


Now, neglecting terms of order 2 of smallness, or smaller (with 


[- ]ds = 


boom 

or smaller) 


[•••]dx for an integrand of first order of smallness 


boom 




in the analysis. 

f * • * ]ds = 


boom 


boom 


g ^ 2 

[**•][! + t( 7T~) ]dx, in a manner similar 

2 oX 


to the one used in Section 4.2 and 4.3, if the integrand is of zeroth 
order of smallness, and if quantities of second order of smallness, or 
larger, are to be retained in the analysis. 

With this qualification in mind, the various terms in the integrand 
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of (4 
point 


.4-2) are computed without eliminating smaller terms at this 


VI 


KKi-. 


0 
0 

x^ w 

r' * "V 

- ^ 

0 


2 (0 ^ S u: 2. 

^ Wj 


/* • 


rv J A ^ r 

~r 

i A (uJ A it ^ - 


W W 

-- Xj co^ w 

^ w w i 

"I 


W 

-X, 

dO ^ w ^ J 


- 5Ci W 


i 

O 
O 


{ T . ^ \ \ 


A ) = 

0 

-V . / . v.."^ 

(w^’4 w.^) 


f? A (t*^ a(w aS )- 


VUJ vx 

1^-3C, (w/.4U}.J 




(Hj;) 


/ 5C, W ol-K + 0(t^j 


= [ + / A 'A, + 

f[l^ ^ (^j-I^) 


J 


CO 6 i) 1 /" 

X jf J : 
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if it can be assumed that x, y, z are principal axes of inertia of the 

total, rigidified spacecraft, of total moments of inertia I , I , I 

X y z 

about the corresponding axes. 

To summarize, we have, to order e, the following equations for 
the rates, in case E, 


I - I 
z y 

— 0) m 

I y s 


I - I 
X z 


(jl) 03 
S X 


0) + ^ 
z I 

z 


wx 


I - I 

dx ^ 03 03 0, to 0(e) 

1 ^ X V ^ ’ 


( 4 . 4 - 6 ) 


boom 


Let the time, t, be non-dimensionalized as t = vi — (from now on, 

1/03 

_ s 

• will designate derivatives with respect to t); the lengths are non- 


dimensionalized by , length of the boom, and 6= ? , Co = n 

w , _ ^x _ ly ^ ^ 

h= j - T“ > k_ We obtain 


X, 


X I 


y 1 


£U _ 
y. ^ 




"z = - 


Z 

J - ky 

kK-f 




J 




X, 




J f I, 


( 4 . 4 - 7 ) 


'■to 
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Using the modal expansion in terms of (O » having eigenf requencies 


io with and to. ^ functions of X and Co > 
j,E’ j J,E 


-xLl. I//'’*'-** 

Now, from Equation (4.2-3i), with (o^ = (o^ ^ 

a. - d) ' 

b M h 


i 


Substituting 


r- - ::;r^ -“,1j 


H^-7) 


AV> 

'4 


T~^ 

Thus the normalized moment of inertia is apparently reduced from the 
value 1, due to the flexibility of the boom, by an amount equal to 


y 


iii 




(4.4-9) 


or writing, with = moment of inertia, about z, of the central hub, 


and 


'i 3 


(4.4-10) 




(4,4-11) 


the non-dimensional inertia correction becomes, in Equation (4,4-8), 

with 

d 0-x -p 5 

zh 

, Z 
f' i 

- ^ ^ 

^-1 m, i 


(4. 4-/2) 
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and (4.4-8) is rewritten (Rate equations for case E ) 



The stability of the motion, in the presence of equatorial vi- 
brations, as studied in Chapter 5, can be done on the basis of 
equations 


(4.2-25) for the modal coordinates 
(4.4-13) for the angular rates 

with N = 1, 2 or 3, depending on the number of inodes retained in 
the analysis. 

4.5 Equations for the rates: meridional vibrations (Case "M") 

4.5.1 Equations for the rates, boom along the direction 

Without repeating the development of Section 4.1.1, the compo- 
nents of expression (4.4-2) are rederived for an elastic 

displacement 6 parallel to axis-z. 

i- W(y)4 


K ^ I Z 


0 

-X; W 


k)/ Iv' V'/ 
^ « 
U}j w w 
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ft K(j ^5 ). 




w 




w 

1 

✓ 







OC:.^ a). 

ftK(wfv(^AS)) 

- 



- Co^(i)._^ W‘ 


0 



63^ tO^ V 

) = 
'K,o7 

w ijO-^ 
0 




^ 0 




^ n(u3 h (to y 






— 

*1 ( 

"v t «z )w 




0 


Neglecting terms of order 

3 of 

smallness , 


-2 . ^ 


)w 


<« 0 >z 


boom 


X, WfO^ f<^i: 


” j, 

^ boom 


X, (w - vv) ( 10^^ + ) f ci)c 


(HrO 


30^ ^ U3^ iV - v^ ._ 6,' W j dx 


ll'z L 

^ boom 

Neglecting terms of order 2 of smallness. 




<Vy= -J 


boom 


X, ( W ^ W j> dx 


with the same assumption as in 4.4 

To summarize, we have to order e, after non-diraensionalization 

I I 

of time by 1 /m , and with k = k - 7 ^ , 

s X I y I 

rt r» 


to - _ 
X 

# 


l-iy 




k 




*' X y v»x>- 


O 
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In the second of equations (4.5-1), 

i 


I r 


I, h 


N 


J -^2 

Using the modal expression for n^ = being the 

J * j 3 ^ 3 

jth modal shape having associated frequency 
} ici ^ ii \w J ' * ^=1 h 


VUCTfly 


Since, from Equation (4.3-16), with m - 1, 

3 


we obtain in (4.5-1) 

-k 


~ W: ^ ^ f ( 3 ^ - ) ~ ^ 




J 4 


JL I _ Ay - 

6)- ^ 

•it \/ 


i l~ ^ - 

1 '-j u _ ^ r 

1,-1 ^ 

/ 

^ Jf 2. 

1 ^'-’-.4 ) 

^ hI 6 kj 

i 

1 

«v 

1 *' 

Ay 



l-^f& 

— ^ '^2 ; / 

^ 

<?. -2k'0S- 

4 J 

ii 

(4.5-2) 


4). . 0 


Investigation of the stability of the motion in the presence 
of meridional vibrations, as studied in Chapter 5, will be carried 
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out on the basis of equations 

(4.3-16) for the modal coordinates 
(4.5-2) for the angular rates 

with N = 1, 2 or 3, depending on the number of modes retained in the 
analysis. Since so far we have been considering a single boom lo- 
cated along the +x axis, it is of importance to generalize the analysis 
to multi-booms configurations. This is done in the following section. 


4.6 Generalization to Multiple-Boom Geometry 

The equations for the rates and modal coordinates were given, 
for equatorial vibrations, by Equations (4.4-13) and (4.2-25), respec- 
tively, and for meridional vibrations by Equations (4.5-2) and (4.3-16) 
respectively, in the case of a single boom located along the +x axis . 

In the present section, we proceed to generalize the developments to 
the case of multiple-boom arrangements located in plane (x,y) (A 
plane containing axes x^, y^, two principal axes of inertia of the 
ellipsoid in inertia of the rigidified, total spacecraft) (Fig. 4.1). 

In order to allow for various possibilities, the following de- 
finitions and notations are used 
I I 

- k , = ^ ~ , k , = - are ratios, smaller than one for quasi- 

X def I y def I 

rigid body stability, which relate to principal moments of 

inertia I , I , I of the total, rigidified structure. 

X y z j to 


given the Etkin'' s number, and non-dimensional radius Cq 

for boom "k”, the notation : 
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- <I>. used for the modal shape corresponding 

3 , k k 0 , k 

to these values of X and (there is no necessity to distin- 
guish between for the equatorial vibrations as opposed to 
for meridional vibration, since they are the same) 




luf 


and all , are normalized to a unit deflection at the tip. 
j » k 


- m. , a function of A , ^ , for given j, i^ the jth eigen- 

J ^ iC 1C V 9 ^ 

frequency for equatorial vibrations .whereas ai. , is the jth 

1 > k 

eigenf requency for meridional Ylbraitions. For the same pair 
(X, , ), we have from Equation 

1C ^ ^ 1C 

. + 1 = (joJ , (all j , k) 

j 3 ,K. 

- q. , is the jth modal coordinate (of type E, or M depending on 

J >k 

which equations contain it) for boom k. 

- is the angle between the boom's undeflected position (an axis 

normal to Z ^ z , thus contained in plane x , y ) and axis x 
P P P p 


of the ellipsoid of inertia. 



4.6.1 2 pairs of booms at right angle, along two principal axes of 

inertia 

In this case, we assume that booms (+x, -x) are aligned on x , 

P 

principal axis of inertia, and that booms (+y, -y) are normal to 
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(4-x, -x) , thus aligned on principal axis of inertia (Fig* 4-1). 

In order to generalize the previously obtained equations for 
the modal coordinates and angular rates, we observe that in these 
equations, 

(x,y,z) are a r.h.s. system, with 

+y in case E 

BOOM ALONG +x, deflection q along 

+z in case M 


Now consider the boom along -x . Equations analogous to these derived 
for the +x boom will apply, substituting 




for 

the expression 



axis 

x 

axis “X 



axis 

y 

axis -y 



axis 

z 

axis -z 




^x 

along z q ^ along z 


4ince 

(-X, -y 

» 

) is a direct system. The 

deflection along -y 

(i.e. 

in case 

E), 

will be measured, for the 

sake of convenience, along 


axis +y, in the same manner as q is measured. Therefore, in the 
analogous equations, written for case E, substitute 
for the expression 

q along +y -q_ along +y (4.6-1) 

~nX X 

Similarly, the substitutions needed are, in the following cases: 
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boom along +y axis Substitute 
for 

axis X 
axis y 
axis z 

V 

q . along z 
^+x 

boom along -y axis Substitute 
for 

axis X 

axis y 

axis z 

q along z 
^+x 

q_l_^along y 


the expression 
axis y 
axis -X 

(4.6-2) 

axis z 
q along z 

-y 

-q along +x 

y 

the expression 

axis -y 
axis X 

axis z (4.6—3) 

q along z 

-y 

q along +z 

-y 


Effecting these substitutions in Equations (4.4-13) and (4.2 iJ)> 
obtain 

Equatorial vibrations (case E) 

It should be recalled that , refers to "E" type, jth eigenfrequency 

J 

of boom "k". Although this is not done explicitly, the’v' could be 
subscripted to account for different damping ratios in the various 
booms . 


Rates 

ft • 

Booms -X, +y, -y: the equations for m m in (4.4-13) remain 

^ y 

unchanged . 

The equations for , in (4.4-13), r<iad 
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•X boom: 


(l. 

+Y boom :. 


r. 


2 1 _ 

i - 1 /hi t 


r « 

'-X 




2 (2i>(3 q fto 


4..J 


r /- -£l 

ur. 

-y boom : > 


^■Z ( ^- — ^ 

/y- r^y/[_y 

Modal Coordinates : ' 

-X boom: 


y 

b- 

z 



). 

V-Y ' 



(4.6-4) 


±L % {^)>u i q ) 

C ly r-' <l''1 Jr1 Jry %-f 


— z 


'j,-x J'-’' ■J’-* 

4-y boom : 

^ i,y ^ 1^-1 

-V boom : 

— z 




‘‘h 

.nil 


UJ 
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Meridional vibrations (case M) 




j>-y 


Again, m, , refers to ”M" type, jth eigenf requency of boom "k” , 
i > k 

and although this is not explicitly done, the v's could ke subscripted 
to account for different damping ration in the various booms. 

Rates 

Booms -X, +y, -y: the equation for remains unchanged, in (4,5^2) 


^ . (4.5-2) 

The equations for u) , oj read, with b, = 

X y 


'"k 

k k 


-X boom: 


to - _ 

X 


f-k> 


to 
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/ -ff 

K <s^' %ir 


'y>h 
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Modal Coordinates: 


-X boom : 

** ■— e> — i 

1),-. " 

H~y boom' 


Uy <»•? %j,y 


— 1 


-y boom: 

4«6,1,1 The Four Different booms 






Let 


. (('■), r . 


k def 


The equations of motion become 






Z 


Equatorial vibrations (case E) : 
Rates ^ 3 j ^ ^ 

i'- 


^2, 


id s, _ 
t. 


k 


OJ 


J 




(d 


Av 


'fL 
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1 ifi hvis: ^ << ]-i 

^ ^ t^' W. , ft. 


i Is 


2V(a^ ^ 


, -2 

+ tO 


2 ^ . * 3 

> ^ ]-2U.[2l>iJ ^ '*'“N| \- ^ , jjl'W <i + io^ a ll 


J-'I 

(4.£-«) 


Modal coordinates: 


For j = 

1,2 N; 

; k = 

1 

+ 


»f 

+ 2)) (0 c» 

— 2. , 

*^1 •> = •“ 

'"‘ 24 ^ 



'4,fe 

^ ' j. 

'*‘^1 i 

For j = 

1 , 2 , ...N; 

k = - 

-X, +y 




■ 0 

* V 1u- 



(4.6-9) 


Meridional vibrations (case M) 

Rates ; With k taking the values indicated; j = 1,2,...N; 

c (/- 1 1 k 1 (ii-i i 

* 4=' ^ ^ ix ix 


L ^ "a 


hf y.y M./ I (4.6-10) 
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Modal coordinates 


For j = 1, 2,...N; k = +x, -x, +y, ~y; 

IT A . ^ fc -Sy) 


%2,^« ({,) -00 J -- 

V /I 


(4.6-11) 








^ “j ) 


‘4'*^ 


4.6.x. 8 Aligned booms identical; different booms along (+x, +y). 
Equations (4.6-8) and (4.6-9), or (4.6-10) and (4.6-11) re 

O 

simplified, in view of the relations 

i *',4/ Ay. 

and similar ones for subscripts y, -y. {4,^- 12 ^ 

Equatorial case : 


Rates : With j = 1, 2,...N; 

- ^ 

£_ vO 

J 


0), . ^ 


Ik 

00 . - _ o3 

% a 


(4.6-13) 




Modal coordinates: 


Vj)] 


Same as (4.6-9), with Equations (4.6-12). 


(4.6-14) 
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Meridional case : 

Rates ; With j = 1, 2 N; 

rj 


w ((-.A- 1 . 

V t ^ 




J- 






.-k, 




2 5 "4.y ) 



; 1 Jf ('-(J ■* )<j - f 

<i'i \ mi ,.<,■. - 
f /- i_ L 1 - 5 . f iiiL + 4- 4 ^ ) 

/ (> ” 


/-w.'H () 

Jrl ^,1 M-i 


' jl 1 1 jj,xl 

~ i/-'* 5 t 


(4.6-15) 


< 0.2 =o 


Modal coordinates : 

Same as (4.6-11), with Equations (4.6-12) 


(4.6-16) 


4. 6. 1.3 Identical booms along x, -x, y, -y 

In this case, we can use in common for all booms, the notations 


W. 




1 ’ '>y 

Thus Equations (4.6—15) and (4.6—16) are simplified as follows. 

Equatorial case : 

Rates: With j = 1, 2,...N; 




to 




J 


a 


J- 


to 

k 


to. 


(/- Ai 


J.' 

5 




(4.6-17) 
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Modal coordinates : With j = 1, 2,....N; 

<rn; 


Q d 


hL 
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Meridional case: 


Rates ; With j = 1, 2 N; 


j ( u I. ii i 4 j. -s (hli - ^ ^ 2 2k) 

3 ( /- A i 2 ) . - “4 — + 4 ^ ) 

+ -L ll /m , / f/- 4)/) (’ 4 , )-^;)u -a )l 

1 JW X 7^,-x/J 
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Modes : With j = 1, 2 N; k = +x, -x, +y, -y 
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An alternate form of the equations for the rates has been used 
in the computer programs described in Chapter 5. 

Let 


K 


1-* 


Ur 


i 


b 






Assume furthermore that the motion is antisymmetric* i.e. q 
n , = -q : Then 

^ -L + 2 ra 

J, T)c -b 2 ^ t _ 


li, t.(r + 4 I 4 4 r il 


I 


T J_ + 2 TA 

il. » Km 

I ;+4ta 

■‘z 


and the rates, as given in (4.6-17) and (4.6-19), respectively, can 
be rewritten: 


Case E: 


U) 
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15 - 
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/- 
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. 0 (4.6-22) 

If furthermore, the transverse moments of inertia of the hub are 
equal, i.e» 

the above equations for the rates simplify to 


Case E: 


Oj t 




(O 


Case M: 


X 


f 2 TA _ 

: At 

-L i2r& 

(4.6-23) 

'* 2.rl 
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4,6.2 booms in x, y plane, necessarily along principal 

axes of inertia. 

The booms are all contained in plane x, y, with x, y as two 

^ A A r 

transverse axes of inertia of the rigidified structure, and are 
normal to , satellite spin axis. With the notations introduced 
in the beginning of Section 4.6, angle between the axis of 

the boom and axis x. 

Let q , be the ith modal coordinate of boom k. The equations 

j 

for the modal coordinates and the rates, written in Sections 4,2 to 
4.5 for the "+x boom", along a principal axis of inertia, will be 
modified as follows; 


4. 6. 2.1 +x boom "k ; angle C. with x 


Equatorial case . 
Modal coordinate: 




% I A - 

n + 2 y to n d t u> q - - — 


Rates 




iO - 

X ^ 


“y = - 


I- -fty - 
iO 


lx 

'^X ” ^ rr 




(/ - (i3 1 ^ 

^ ' I. 


(ffVa 


31 tf 
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(4.6-26) 


1 — - 1 n . . (j . f to ^ « q p 
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Meridional case ' 


a , +2>>(0 d i -fto Cool--w^‘^L 
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sxnce 


w in (4.3-16) becomes - u> 

y X 


+ to <^d 


0 ) in 
X 


(4.3-16) becomes 



Rates : 

In the case of a "+x" boom, Equation (4.5-1) shows that the 


vibrations parallel to the z-axis generate a torque along the direc- 
tion normal to "+x" , i.e. "4y" having projections: 


The equations for the rates will read; before dividing by 


respectively, 




T 

J ^ 


- 2>>i3 


‘ V'A (“-!*} z, ^ (■“< ‘-i(. 4“-Sn-Zi«ii. 


- 0 



Now define the following coefficients 
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Let 




System (4.6-28) is rewritten 
^ ^ 22 . 

.Let 



Then the equations for the rates are 



» 

W . C 
z 


(4.6-29) 


4. 6. 2. 2 General case; "B" booms, making angles C^(h = 1,...,B) 

with X . 

P 

Equations (4.6-28) will, in the general case of B booms, at 


angles C. (k = 1,...,B), have r.h. sides with sums over k, in 

fC 
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addition to the summation over j . Important note : all modal dis- 

placements are referred to the +z axis (case M) or to the normal 
to the boom (in case E) such that (x^, z) is a direct 

system . 


Equatorial case : 
Modal coordinates: 


— 1 


Cl t 2 W . , o -r ^ \ ■ W r 

in which expression (4.6-31) is substituted. 

Rates : 


- 


— — 


- — 


Ax 

u> 


1 

i-y 


_ J_ 
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(4,6-31) 
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Meridional case: 


Modal coordinates : 


h 




a 4-lHu; <i +*2^ * (-'®„ '■'^y " “x 

V V U.K ^ V > N 


- ^ 


(4.6-32) 


in which expression (4.6-33) is substituted 


Defining 
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The above formulation is the most general that will be considered 
in this work. 


4,6-3 System considering meridional and equatorial vibrations 

simultaneously. (B booms in X , y plane, not necessarily 
along a principal axis of inertia) . 

In the Lagrangian formulation, for an elastic displacement 

V 1 H- f f 1 oU t ou^) 

' (la 2 J V j 1 

E M 

and the kinetic energy is, if , hw =* *r~ , etc., 

k k 

I 




T, T I (7 ! n] 

a /!u -I 4o«. \k 

V ‘ '' 'Jt ‘ -* 




‘H 

It can readily be seen that when are expanded in their modes 

?o)» with associated frequencies o). , , O), , , the corresponding 

J >b 1 >k 


modal equations for 


M 


>k ,.J< 

E M 

are uncoupled. For ”E”, only u)^ will appear in the r.h. side, and 
this quantity is a function of the q. , only. For "M” , only , w , 

J^k * y 

^x* % will appear in the r.h. sides, and these quantities are 

functions of the q . Hence, in the total system, 

J > 

M 

• • 

- the two first equations of (4.6-33) are those for m , aj 

X y 
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- the last equation of (4.6-31) is the one for w 

- the modal coordinate equations for q. , are given by 

j 

M 

(4.6-32) 

- the modal coordinate equations for q. , are given by 

1 

(4.6-30) ^ 


4.7 Conclusion 

The equations of motion of the spinning spacecraft having 
flexible appendages have been derived in a rather general case, 
using the modes of the rotating structure at the nominal spin rate, 
and for a central hub of non-zero radius. They x-jere found to be 
in agreement with some other published results in the limit case 
of a central body of zero radius, and can be used with profit in the 
numerical simulation of flexible spacecraft motions. 


I 
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CHAPTER 5 

Simulation of the Satellite Attitude 
Motion and Stability Studies 


5.1 Motivation 

In the present chapter, we present a simulation study of the evolu- 
tion with time of the satellite attitude, from which stability charts 
can be obtained for use by the satellite designer. Of particular 
interest is the "nutational divergence” phenomenon, in which the 
satellite, although stable if it were ”quasi-rigid , exhibits a stead- 
ily increasing nutation angle. Its spin axis thus drifts ax^ay from 
the invariant angular momentum vector, on which it is assumed to be 
aligned initially. This instability is due only to the dissipative 

motion of the elastic appendage. 

To this effect, a set of computer programs, "FLEXAT”, has been 
developed which numerically integrates the equations of motion and 
prints or graphically outputs the variables of interest. This pro- 
gram quite markedly differs from earlier versions we have used in the 
x^ork, as will be explained later. The version given here accommodates 
thx'ee modes of the rotating structure and a diss^etric central body, 
and since it permits an easy visualization of the qualitative features 
of the attitude motion, it should appeal to the satellite project en- 
gineer . 

5.2 A Package for the Simulation of the Spacecraft with Flexible 

Appendages. 

5.2.1 Generalities 

FLEXAT is a set of programs, written in FORTRAN V, which were 
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. mostly run on tho UNIVAC 1108 at Carnegie— Mellon University. It is 
composed of the following parts: 

a) A short "MAIN” program calling on the relevant SUBROUTINES. 

b) A subroutine CASEM 2 called upon to study the stability 
of the meridional vibrations. This subroutine internally 
calls on its own subroutine RATES , which computes the angu- 


c) 


lar rates m , m , o) , 

X y z 

A subroutine CASES 2 called upon to simulate the equatorial 
vibrations. Again, this subroutine internally calls on its 


own subroutine RATES , which computes the angular rates to^, 11)^. 
In particular, this subroutine can be used to simulate the 
nutational divergence occurring when the GMI (greater moment 
of inertia) rule is violated, for the rigidified body. 


d) A subroutine SEARCH (NDS) called by the MAIN program and 
yielding the eigenfrequencies (up to j=3, if required) of 
the rotating structure, corresponding to the specified values 
of A, ?o* This subroutine, for the essential part, is the 
same as that described in Section 2.4. 

e) A subroutine PLOT, called Internally in either CASEM 2 or 
CASEE 2, giving a graphical output of the evolution with 

time of the satellite nutation angle, over a number of satellite 
spin periods (generally taken to be 10 to 20) . 


Each of these parts is now discussed in more detail. 
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5.2.2 Program MAIN 

In this program, COMMON, DIMENSION etc. are given. Then the 
"unchanging parameters" are specified by cards. The listing given 
at the end of this chapter, for example, specifies 

NSKP : skip the printing of 60% of the results is 

desired (NSKP = 1) else -NSKP = 1; all results 
plotted in both cases 

N0RU=NSUP=3: include 3 modal coordinates for each boom. 

XN0(1)=O.O5: the "x-boom" and the "-x-boom" have modal de- 

flections (1st mode) equal to + 0.05 times the 
length of the boom 

XN0(2) .. .YN0(3): thel"x-boom" and thei" :/-boom" have zero 

modal deflections, for the 2nd and 3rd modes. 

NU(1) NU(3)=0.05: same damping ratio on the 3 modes 


CASE = 'M' ! 

meridional vibrations 

SIO = 0 

value of 

LAM = 10 

value of X 

MGIV 

a switch. If equal to 1, the eigenfrequencies o)^ 

and m . , m are given as data (they are assumed 

1 >3 '2,3 

to be known from a previous study, or from a table). 
If equal to 0, the and the other quantities will 

be obtained "on line" by calling SEARCH (1) (in 


case E) or SEARCH (0) (in case M) 

GAM 

: r in the developments of Chapter 4. 

PKX,PKY 

: ratios K =I ~ u/ ^ u u 

px Zjhub*^ /:,hub py z,hub'' /,hub 

These measure dissymmetry of the ellipsoid of 


inert ta of the central body. 

PREC: 

: the integration interval in time is equal to 


— — or — (with j = PREC) whichever the smaller 


It has been found sufficient to take PREC = 1. 
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MAXP : maximum number of such periods (defined under PREC) 

to be considered. 

MODES ; 3 (should be the same as N0RU, NSUP) . Three modes 

are retained. 


5.2.3 Subroutine SEARCH (NDS) 

This subroutine has already been described in Chapter 2. It 
obtains in the relevant case (E or M) for j = 1,2,... NSUP. Note 
that 

a) NDS is an argument given in MAIN (0 for case M;1 for case E) 

b) SEARCH is bypassed if MG IV = 1, i.e. if the eigenfrequencies 
in the case of interest are externally given, other than 
completed on line. 


5.2.4 Subroutines CASEM2, CASEE2 

This subroutine, fed with the w , , m . » . values obtained from 

data or computed in SEARCH, proceeds to integrate equations (2.2-8) or 
(2.3-5), as the case may be, if MG IV = 0 , and bypasses the procedure if 
MG IV = 1 . 

It then proceeds to compute the quantity 
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the same geometric and structural properties (thus same 

to be along the principal axes of inertia of the rigidified structure 

= 0, /c for the x-booms, 1 ^ y-boom, in Chapter 4). 

The ellipsoid of inertia need not be of revolution (K ^ K ). Thus 

px'^ py 

the relevant equations have been written as equations 
(4.6-20) and (4.6-24) for program CASEM2 

(4.6-18) and (4.6-23) for program CASEE2 


Different assumptions (booms of different length, structural proper- 
ties) could easily be considered by the user, for any special applica- 
tion, after a rather simple rewrite of the equations, as given in 
Chapter 4, or a suitable distinction between etc, 

rather than the common “r"... adopted here. 

The method of integration is RUNGE-KUTTA with fixed step , 

the latter being computed in the program as some function of the spin 
period or of the vibration period of the jth mode, as precised in 5.2.2. 
under "PREC". 


The output consists of a .print of the case data, of the ' quantities 

NSUP ^ [jr ] 

Z ; A, vi, V 2 , V 3 ; m , m . , m (j=l, . . .NSUP) ; H initial= 

3-1 j 1,3 z,3 

(assxmiing Ho and z are initially aligned) ; then tables giving 


‘x,l 


Ho Ho Ho Ho^„ 

y,l ^x,2 ^y,2 ^x,3 y,3 x 


e 

(angle of nutation, 
degrees) 


ST^P 
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There exists an option to skip the print inf» of the first 60% of 
the results over the time interval considered, which makes sense if 
one is only interested at looking at the long-term behavior. 

5.2.5 Subroutine PL0T 

The PLOT routine graphically presents the results of the above 
computation. PLOT is internal to CASEM2 or CASEE2, as the case may 
be. 


5.3 Results from simulation study, using FLEXAT 
5.3.1 Comparison between the present and some previous results 
As compared to the approach previously taken by J. Rakowski 


and the present author 


, the equations used in the present 


O Q 

simulation do not include "extra" non-linear terms such as qS q^, 

mm... Including these terms, although they appear in the deriva- 
X y 

tions of Chapter4, did not seem. • fully consistent with writing the con- 
tributions to the kinetic and elastic energy with some terms of order 
3 of smallness neglected (such would be the case, for instance, if 
[{••♦}dx = [{♦•»}ds , with the integrand of first order of smallness). 


However, strictly for the sake of comparison, the stability 
boundaries, derived as explained in 5.3.2, were compared in a large 
number of cases using, on one hand, the equation with the extra non- 
linear terms, and on the other hand the equations obtained in Chapter 4. 

In no cases were the differences of much significance. All were 


well within the sampling interval (K + .016). 

P 
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5.3.2 Parametrization of the stability chart 

r 5— 1 5—2] 

Following the notation adopted earlier^ * , it is proposed 


to define a stability chart as follows, in the symmetric case 


(K = K = K )(See Fig. 5.1) 
px py p 


abscissa: K = P ^ a measure of the asymmetry of 

^ ^z ,hub 

the ellipsoid of inertia of the central body. 


-ordinate : 




a measure of the relative importance 


Z ,b 


of the inertia of a boom (- if ?o = 0) > and the inertia 

of the hub. All things being equal, small booms of small 
mass will give small values of V . 


- parameter, of the plane: 

= fixed non-dimensional radius of the hub (referred to 
the booms length) 

- parameters of the curves: 

- 2 ^ 

X = — 0 ) ^ ( ) , a ratio of certrifuqal to elastic 

El s CO ° 

cant 

forces, large for high spin rates or very 
flexible booms (E, I small; large) 


Thus X = constant curves will be drawn on the (K^,r) plane, 

for = constant, corresponding to the observed limit of stability, 

i.e. a point, at given T, , X, such that any slight increase in 

K causes stability of the observed motion, the nutation angle tend- 
P 

ing asymptotically to zero; whereas to the left of it (decreasing 

K ), the motion is observed to be unstable, the nutation angle 
P 

steadily increasing with time. 

In the asymmetric case, one more degree of freedom exists, and 
the chart will draw X = constant curves, corresponding to the observed 
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limit of stability, for given s K , in a (F , K ) plane of 

py px 

representation. 

5.3.3 The GMI rule 

As described in [5-3] , a rigid body undergoing a torque-free 
motion about its center of mass, but having interh.al energy dissi- 
pation, has a stable spinning motion only about its maximum axis of 
inertia, i.e. if 

I I 

and > 1 (5.3-1) 

X y 

If one of these ratios was one, there would be no preferred axis of rota- 
tion about which the satellite would spin after an initial nutation 
has been removed by energy dissipation. Condition (5.3-1) is commonly 
referred to as the GMI rule (or greatest moment of inertia rule). 

In the stability chart, planes described above, condition 
(5-1) will be represented, in the S3unmetrical case 



by a locus of equation 

2FA > ^ - 1 (5.3-2) 

K. 

p 

or 

4 +-0 > f ' (5.3-3) 

These curves will, whatever the value of tend to the common point 

r 0 K ^ 1 
P 

which they should not include. This corresponds to the case where 
the satellite has no flexible appendages (p^-^ ^ 0) and a spherical 
ellipsoid of inertia. The curves are shifted to the left as in- 
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creases (Fig, 5.1). Their X parameter is X = 0. 

Conclusion 

For the stability of the satellite with perfectly rigid appen- 
dages, and of the satellite with flexible appendages in the presence of 
equatorial vibrations (as explained in iT.3,4 ), the greatest moment 

of inertia rule 

y ^ X j 

should be satisfied for the total , rigidified satellite. On the 
(r , K ) stability charts, the design point 

for given 

should be to the right (i.e. in the region not including the origin) of 
the Quasi-Rigid (QR) locus given by Equation (5-3) . 

5.3,4 Stability with equatorial vibrations 

Stability in the presence of equatorial vibrations, was found 
to be equivalent to quasi-rigid body stability. The stability condi- 
tion for case E is thus the same as the Q.R. body condition given in 

Equation (5.3-3). This result. is in good agreement with 

F 5-41 

Hughes and Fung analysis in the case where = 0. Two examples . 

are given in Fig. 5-2 and 5-3. 

5. 3. 5.1 Stability charts (case M) , using three-mode analysis 

Using the FLEXAT program with subroutine CASEM2, and retaining the 
three modes in the simulation, figures such as 5 t 4 to 5.7 can be pro- 
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duced. Each of them corresponds to the same value of T- 10 and 
Eq = 0.1. For X = 100, two values of K are considered corres- 

p 

ponding to a slightly tiastable or a slightly stable condition (Fig. 5.4, 
5,6). The same applies to a higher X case (X = 1. ,000) (Fig. 5.6, 5.7) 
The final results of the three-mode stability analysis in the 
presence of meridional vibrations are summarized on charts 5-8, 5-9, 
5-10 for values of X = 0 (Quasi-rigid body case) to X = 10,000, and 
for Eo = 0> 0.1, 0.25. 

IMPORTANT NOTE : When using program FLEXAT, with subroutines SEARUI 

and CASEM2, for X K 5,000, the values of the relevant frequency and 
modal quantities: 


<0 






> ’^'2 


(5.3-4) 


should be given as input data, using option MGIV - 1, or described 
in Section 5.2.2. Quantities (5.3-4) cannot be obtained on line 
using program SEARCH DP, for such high values of X. They have been 
obtained using a multiple precision version (OS-MP or NP-package) of 
SEARCH, which is rather time-consuming and should be run only to set 


up tables such as in Section 2.8, for interpolation purposes. 


5. 3. 5. 2 Effect of higher modes, and of modal truncation 

As the tables in Section 2.8 show^ the effect of higher order 

modes (j = 2,3 ) on the motion parameters is as follows: 

2 / 

^.1 S' ’’J 

For small values of X, the changes of this sum by increasing NSUP from 
1 to 2,3 is at most 2.5% for Eo = 0, and 9% for = 0.25. 
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For large values of X (I = 5,000), the corresponding changes 

are 0.03% for Ko = 0, and 0.5% for Co = 0.25. 
m , 

(amplitude in r.h. side of jth modal equation) 

"^l,j 

It can be seen that this ratio is at most 25% (for j = 2) of the 
value corresponding to j - 1» when j is increased to 2,3 . 


c) m^ . (amplitude of some terms in the r.h. side of the rate 
equations) . 

The same comments apply to m. . 

> J 

To assess the effect of higher modes qualitatively, it should be 
remembered that, when non-dimensionalized by co , 






and the forcing frequency (precision frequency in body-fixed axes) 


on the terras would be, for |q [ , |q 1 « 1, 

ft. 


as opposed to 




i 
■i 


< 


W 




/ _ t- ZTA 

< 


for a quasi-rigid body 


mj 

Note that is always smaller than A . Typically, for 


4,0,, . 0.^4.? 

llL ^ 0 . 4 t ^ X r 0 Or 43d % - /OO 3 437 |o^X - \j0Q0 .. 
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Therefore, in an approximate sense, it can be said that angular 
rates will not appreciably excite modes 2,3,..* which are larger 

r 

than uj| by a factor of several units at least* 

With these observations in mind, we now discuss the conclusions 
of a detailed study of the effect of modal truncation on the stability 
charts (f, K^; constant X, 

Xt was indeed observed in the simulation that higher modes never 
developed to amplitudes of more than a few % of the amplitudes of 
the first mode, assuming i.e. which can be considered as "normal for 
the initial deflection, namely close to the shape of the first mode 
$ 

Within the accuracy retained in establishing the stability 
charts (K + 0.015), no noticeable difference could be reported between 
the stability chart determined here on the basis of three modal coor- 
dinates for each boom, and that we obtained on the basis of a 

single modal coordinate .* Se • ing times, however, were larger. 

The results of the 3-mode analysis, using program FLEXAT, are 
summarized in .. Figures . 5.8, 5.9, 5.10. 

5. 3.5. 3 Effect of some higher order terms 

As was mentioned in 5.1,. there was a lack of consistency in re- 
taining some non-linear terms of order 2 in the equations and neglecting 
some others. Equations (4.6-20) and (4.6-24) x^ere used in the present 
stability simulation. It should be noted that little difference re- 
sulted in the stability charts. The angles of nutation, however, are 
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computed here by 

■/h ^ + H ^ 

0 = nutation angle = arcsin ( — ^^77 = 0(e) 

tot 

and since they involve quantities of first order of smallness, should 
be accurate, whereas the use of formula 

will see 0 critically effected by terms of O(e^), none of which should 
then have been neglected. 

5. 3.5.4 Parametric studies for I 9 ^ I (Ellipsoid of inertia not 

X y 

of revolution) 

With the particular geometry considered here, 

I < I < I implies that 
y X z 

or ^ ^ 

A set of parameters is chosen, namely 

Co > ^ number of modes. 

In the (K , K ) plane, the bisectrix of the first quadrant, K = K , 
^ py px P" 

will correspond to the symmetric case, 

K = K = K 
px py p 

and the limit of stability K ,such that K > K will, ensure stability 

P« P 

of the motion, was found previously. Furthermore, in order to satisfy 
the GMI rule, we must have ^ 
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In 

program 

or 

and 


order to determine the parameter region to be studied with 
FLEXAT, it is useful to note that 


- I < - I.yA ( 

' ' <1 


k 


/■j 


k 


h* 




k^, 


>1 


Similarly, from 


• ^ _11 < \ 






K 






This is most conveniently represented on a ^^^py^ plane. 


(Fig. 5.11). Thus, if 

X, H 1 


X. J. 


k 




•M 


the admissible domain of study is bounded by 

-X, < I • X, <i+2rA 

X*-X<l 

^px 

In particular, for a constant ratio of (or » the limits 

py xK 


are shown by circleSon Fig. 5.11. 
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5.4 Conclusions 

A program has been developed for stability studies and simulation 
of the nutational motion of a spinning satellite with flexible 
appendages. The results of this program can be used with profit in 
the preliminary attitude design, to ascertain stability, determine 
the importance of structural damping and study the rate at which nuta- 
tion is generated or removed from the system. 



r 
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FIG, 5-1. QUASI-PJCGID BODY STABILITY 
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FIG. 5-10. Stability Diagram. Case M. 3 Modes. 
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A r-j D the i-j in a t r 0 f-j a f : o l i i h d E g r e r b 

L» J NS I (Ui Oi-HJi:i) 
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SFT PARAHETEin^ YO COtiTKCt SiIUIlAtIOM 
SET Nf)PT=l FOR KEVtK5LD INTEGRATION 
hi 0 S 3 0 i R t C T 1 nr! S ‘‘-I I T CFi 
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KO-0* 
d 0 - G • 

i F ( I^OPT . GT . G J C-O To t 2 
t 3 ? p ( 1 ) = J » 

B n a F 3 2 P ( ] ) 
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1 2 I- 3 ( > ) = I . 

L 0 IT B 3 < 1) 

! 3 C. f) a Q B 
^ 0 ~ 0 » 

c. 0 T 0 3 ... 

2 A 0 = 0 * 

0 0 = 0 • 

1 F U'HJPT • C-.T . G ) (lO 1 0 1 

L q 3 p « J J = 1 . - 

C 0 = E 4 3 P ( J ) 
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A 0 a r H 1 p ( n , - 

lb b 0 = 0 * 

t-O-G*.. 

3 A = AO 
t’ a B 0 
C s c; C 
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N N = N 
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b B ( f ) = « A 
L ( ) ) = »^ * b 
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Hu ] =1 . + TM; ( HOOL 1 «•• H 0 ( ; t () q p ) 

i F ( N 0 b • E 0 . 1 ) G 0 TO q D 
Hut=MMl“l. 

q n P { 1 } - < < 1 * - S T S 1 + 2 f * s i 0 * n • - S I ) ) / 2 * * r. -• ( S I ♦ b ! 0 ) « A - U 1 » E ) * L . A M * H 

I p ( M 0 P T . G T • t .n P < 1 » = • ( 1 - S i * 1 + 2 * * S 1 * 1 1 ♦ + 5 I 0 ) > / 2 , * B 
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B = H ( Z 1 /2 , 
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Go TO b 
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fc- 3 t 7 ) * I 3 < r j ) ••• < K < 1 < k < 2 ) ♦ K ( :U + K ( 'n ) / 6 • 

E 3 1 P ( 7. ) = t :u P < N ) + I L . ( i ) + V • « L ( 2 ) + 2 . » L ( 3 ) + L ( J ) / 6 . 

C 3 ? t 2 ) = F: 3 2 Wv > + t M ( \ ) + 2 • ’ f*M 2 > + 2 . « M ( 3 ) + M ( '< ) )//'.* 

E 3 3 P ^ 2 ) = t 3 3 P ( I j ^ ^ ( P ( ) ) + 2 * ^ P ( 2 ^ + 2 • * P ( 3 ) + p ( M ) ) / fr> • 

1 3 M P t Z ) =L A H * ( { I S Ki -M . } ♦ « 2 - ( S I + S I 0 ) » * 2 ) * C 3 2 P f 7 ) / 2 . 

1 ~ ( 5 1 + S 1 0 > ♦ f 3 i P < 7. ) + f'l L' I * E ^ ( 2 n 

i r ( MOP T . r, r . p ) L 3 ^ p 1 7 ) = i a m * ( ( ( s i 0 + 1 • * « * 2 - ( 1 » - s u s i 0 ) : 

] ’ r 3 2 p < 2 ) / 2 * + ( l « - S I + s 1 0 ) * E 3 1 P ( 7 ) H u 1 , E 3 ( 7 ) ) 

L S F-- 3 < N -*■ 1 ) 

A = f: 3 1 p ( K *M ) 

B = p 3 2 P ( N 4 1 ) 
c C r. 3 3 FM r; + 1 ) 

Eo = £ 

A 0 » A 
B 0 = & 

0 0 = C 
Ns;!.;+ I 

1 p ( t-: • 1. T * I N T P R ) a T , ) ^ 

t- M 0 f.‘, 3 = E 3 2 P ( I f-n E R j 

I p ( F'.’ (i f- T » fi T • Cl > E !•: 0 M 3 r E 3 ( I l■•l T E K ) 

Do 3G I^UtlNTF.P 
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9 Si=;rif]*H 
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E q /-i P ( 7 ) = L A F't * 1 ( 4 3 I C + i • ) * * 2 ~ 1 + 5 I 0 ) * * 2 ) « f ‘i 2 P ( 7 ) / 2 . 

1 “ { S 1 ^ ‘i 1 0 ) * E 'I 1 P * Z ) + lUi 1 » K *4 t 2 ) ) 
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^OsF 

BO = A 
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Do 2 16 I - 2 t J 0 1 
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2 I a s H = s H + I n f-’ T I I , iu> n t ) » (} p T ( 1 , m o p e ) + 'i p t < ! - i ♦ m o n e > ^ h p t ( i - 1 i ho lu- > ) / 2 « h 

Ami t H 0 1:) E ) = 5 m 


I 

c 

I 




^ 'i Co ij T I f'i u E 

EMU or POUE 5HAPt: AND MOOE PARAmETEP r A |. C Ul. A T I 0 N 

All VAKiApi.Erj are jvior*-p i mens i ofJAt 
N P I s the step f) U M h y 

Hpj IS THE internally calculated PKiNT inters/ A t, 

NMaX is The HUMPER oF segment^ jN T^E smallest PpRIDO 
Nx DET EF^ 1 NFS HoW MuCh Op ThE SIMULATION IS PRiNiED OUT 

0 A M A = R U 0 ^ L * «• 3 / 2 - H U B - I 

Co« = SUM Q\JEH M MOUL5 OF M2 SQuAEKD r.VER Ml , 


3 H □ Continue 

K I'l A X = 7 S 

I R = Li IM r R E c ) 

If ( pRE C*RE . 1 ) Up I =Np ! • U< 

UpRQ=f!U ( PREC > . - - 

PPAC“ 1 • /f LOaT ( KhaX ) 

P P R = 0 . * ( ;i . T A lM I • > I 

I p { R P R * G T » 1 * ) F’ E K = p E, R / 1*' U I P R E c J 

h=pek*pfac 

M K r: M A X P K M A X 

1 F ( P K E C . N E . 1 ) M K - M A X P ♦ K M A X • I R. . _ 

nx=*^*nk 

N P I ( M K - N X ) / 6 D 
C 0 f\ = u * 

Do 20? Jvi= 1 » houf;s 
A M 2 I ( ..MV ) = A M 2 ( J iV ) / A i'l I ( j ) 

A j M ( JVO -7 . » aH 1 UV. ) 

A7i-, ( Jvf ) = 7 • * AM? ' > 

20/ CoR = COR + Ail2 ( Jo J ( jL' 1 

D p L T A = S I 0 * S I 0 + S I 0 + o 3 3 3 3 3 3 3 
DnElT=OELTA 
U C 0 R = C 0 R 

WRIT £16,971 DCOr,E)DELT,NU 

9 7 FoRMAKir! , * C 0 E- * » F 1 1 , S / I H ,*oElTA=.,P 9.S/IH , » N U 1 = * , p 1 n • / 1 H , 

X ' |jU 7= • » P I 0 • M , / 1 H , * iNjU3= • , F' 1 n . H / // / 1 

Do 2 03 Jl;- I, MO Ups 
0 A M 1 ^ J 1 " A M 3 { J *■'' 1 
0 A!12 { J''-' 1 = AM2 ( Jl^‘ ) 

UMIJ ( J'-' 1 = MU 1 J ''«■ > . 

U f.| 2 ( J L' 1 = 7 . hi LI ( J ii 1 ** H u { j '•/ 1 
R I T F! ( G , ■'I Li 0 ) J G , 0 1 1 U ( IT. J I J - J , I) A M 1 I J '-'O » j V- , 0 A M 2 ( si •'’ ) 

Apn FOitMATFIH — 1|*F)>SX, 

X • !1 2 * ♦ 1 1 , ' . E 1 I • 5 / 1 

2 03 COiiTlI'-UE 

l>0 'lOG J'A =:I , MO Dps , , 

H Q P \\ ( J W ) - M U ( J\^ 1 * M u ( J VV ] 

G A LI A 2 " 2 • « gad ii 

2 }<; t si * + 2 • « G A M A ? * 0 E L T A 

3 6 A I = i . / P K X 

Y 1 s 1 • / P iC Y 

t 


U E r I E! E initial C 0 N & i T i 0 N S 
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0 ri X n - 0 • 

om Y 0=0 . 

0 M Z 0 = I «» 

0 n r fj A = 1 * 

OmxdO=u • 

J M Y 0 0 « 0 . 

0 M z n 0 = u « 

ThI^sE l.C* ALr'irj Z-aXiS on H-t^ECTOI^ 

t)Q HQH JVj=I,MOOF.S 

0 M X 0 = 0 r-i >; 0 + G A f I A ** 0 n Z 0 ♦ X ''I 0 ( J VV ) / I XI + 2 . ® G A A ® D F L T A ) 

M n H U f.t Y 0 = 0 M Y G + G A M A Z » 0 Z 0 ^ Y I‘) 0 ( -J / ( Y 1 ^ 2 » e- G A M A * D F L T A ) 

D() m05 

X M [) 0 ( J W ) =: Q • 

Yni.iO ( JtV ) = 0. 

V { t , J VM =: X N 0 ( J ft > 

V ( 2 , J* > sYhO ( Jft ) 

V { ,;i , JW) =XNDO { ) 

^ 0 B V ( M ♦ J •’'■)= Y N D 0 ( J > 0 

V ( b . 1 » = 0 F: X 0 

V ( A , n =OMYO . .. . ... . - ■ ■ 

V ( 7 , I ) = 0 z 0 

V ( .A . 1 > = 0HXD0 . ■ 

V ( 9 , n = omyoo 

V ( 1 n , 1 ) = OltZDO . 

Ij Q 9 J = 1 , I 0 

Do 9H jrisl iflOOEs ... - 

V G 0 V ( J I . = V ( vl » v O - ) 

N=t 
P = ! 

R A i = 1 .......... 

b U C 0 N T 1 f O E 

K ij N G r K U T T A I fj T E K A T 1 o N ' ^ 

Do ‘407 1 ,H0DFS 

A G X ( J '■' * ” ** 0 M Z ^ J •' J * D V ( 3 » j v'si 4 - 1) v '' .1 » 0 P A= < J 'V ) - A 2 1 ( J .v 

X ( 0 E G A « D V ' b » I I - 1:» V r 1 1 M 

Ac Y ( J'7 > =-*U02 ^ J-' ) *0V ( ‘5 » J *0V < 2 » J'’’'- J ‘*0^3 I J'/‘ ) -AM? I f f * 

X ( u M E G A * L> V u t n + lUM y » 1 n 

A K 4 1 * ) » 0 ' ' ) ~ H ** D '</ ( 3 » J /i 4 

Ak ( 2 > I > -OV J = M«0Y' *; -4 , Jiv ) 

AK I 3 » I • ) =H^ aCX ( JW j 

H n 7 A K ( H » I t J l( •' } = H * A C Y 4 J vV ) 

A K ( b t I ♦ 1 > = H s* i) * I ) 

A K ( A » J » 1 ) = H * O ''/ ( *7 » 1 ) 

K ( 7 » I » 1 ) = H » D V ( 1 0 » I ) 

I = J + 1 

I F ( I » G T . 3 ) G 0 T" 0 G Q 
Z=l-1 

Do 9 i J= J . 7 

DO 9 I JW-IiMODKS 

9 I D V I J I J V>( ) = V ( J » J + A K t J . Z t vJ W ) / 2 . 

C A 1. L rates 
G o TO GO 

6 n 1 F 1 I ® f > T . ) GO To ‘7 0 

U 0 1 0 = 1 * 7 

i.»o 10 J.v=-i » modes 
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0 V ( .i » J i'i ) = V { J » J ) + A K f J , 3 t J ) 

'^ALI. f-^ATCS 

‘JQ TO bO 

'^n f^O ). 1 J==] »7 . . .. 

i)n ti ovi a i , M0pf:s 

t-' C t J * -OV ) - K ( J J 1 t J W ) + 2 . e ( y\ K ( J I ? » j iv M A K ( J » 3 » J IV ) ) + A K < J » i j W ) 

1 V ( j . j 1 a V ( J , J vy M- c C ( j , j ) / 6 . 

UO BV J=^I,7 

'^0 BV JiVa I , Moor.s 
0 V ( J J JW ) = V ( J t J-- ) 

Call katls 


Kt.M'^GP KUTTa FlNlSiiLO 

no A calculate output vAkjABLLS 

1 F ( N P * N f * J ) C 0 T 0 J 

CoHpUTt COUPOtJi:iiT 5 or H-VECTOk |M BoDY-FlXLO AXES 
Am X = V f ^ » 1 ) ® ( X 1 ■»■ n A M A 2 r L T A ) 

A M Y a V < o , } ) » ( Y I G A • I A 2 « n r L T A I 

L>o MSI Kf= 1 , USUp 

An t a A M X + A N 2 ( K F ^ 4 ( V 1 H f ) - V ( 7 , M # V ( r , < F ) ) * f, A M A 2 . 

5 S I ^ M V “ A Y ~ A >-1 2 ( K F ^ ? V ( 3 , f;, P ) + V { 7 I T' ) » V ( ? I K F ) ) * G A P A 2 

AM7 = 2'^I*V(7f]l, ) 

fi S 0 = A M X f A l-t X + A H Y t \ M Y ^ A f-, Z *• A M Z 

f i S Q = S Q R T I H S t} ) 

U H X Y - A '“5 X * A M X + A 0 Y A M y 

Ohx Y-SQirr ( oux Y ) 

i-* I = 3 • 1 ^ 1 S 9 2 ■ 

1 ea • /p i 

C A - D E G » A 5 I 1-1 ( 0 H X Y / H S fj ) 

A F ( N S F. P a F' Q * n I i', 0 To M - 7 

I F ( N * 0 1 « 1 • 7i N |) • U a L T « i\i X ) G 0 T 0 Ml 

3^7 c o N T I n U L 

i F I N • nL • 1 ) GO To 29 
-“R J TF ( 6 , 33 ) HSU' 

3 F0R0AT(/1h jMIMTIAl h=',F12^/, ) 

I tF ( 6 , 9M ) 

M f- QRMAT ( //// 1 M ♦* Qx (-lOOL'i OY MOr.Li OX M0 Df:2 qY MonE2 

X X r i 0 D 1 3 0 Y r-f 0 1'J C 3 0 M L G ^ X C A N * t / > 

79 CoNT INUE 

1 F f ! I P * N F . 1 I GO TO Hi 


lYP I T F 0 U T }-' ij T V A K i A a L £ s 


b 2 2 


M I 


R 1 T E { 6 , 2 2 V V n M > < V ( 2 . I > 1 V ( I » 2 ) ♦ 7 f 2 ’ 2 ) * V <1 , 3 ) » V ^ 2 t 3 ) 1 V t S 

1 C A t TJ 

E 0 P f -1 A T ( 1 H * 6 I 0 M ^ S , » X ) , 0 ! 2 , 6 , i X » p I I , , 1 X , I S ) 

!'i p _ pi p + I 

i F < nP * CQ * ,jP I J Up a I 

Call plot 

iF{n*GT*MKJ RET'jRfJ 
Go TO B B 


1 J t. 


^iUPROUT I ivF. rates 

This KOuTiMF calculates the OfRIVATtVES of the angular rates 
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<: 0 - ( H 0 1'J V , D V » C 0 » ! 1 A M 2 , X J , Y l » D C I. T A i A !1 ^ 2 . I t U N 2 » 0 M f Ci A 
CoMMOf^ / /: KN /OMU 
C 0 M U N / T .< K C I: / (U» D E *i 
C o < • H 0 \' / T I U R D / I'J S U I* 

';f' NS 1 Oli OM'..l{.U 

i>0 ! f 0 L K P R t; C I S I ON V ( i D , 3 > , 0 V t 1 n , 3 ) » C R . A H 2 ( 3 ) , 

Ij i J 2 ( 3 )^ , G A i > A » 0 II L T A . X 1 I ^ 1 « ' 0 1 1 3 ) » A • ! X t » A U X 2 

\)0 7 11 J 1 - 1 »OSOp 

7 ? I .lu { J n -oMu ( j n . 

A u X l = - I . / < X 1 + G A h A 2 « ( OKI. T A •“ {• 0 R ) ) 

Dv {« , 1 ):^AUX \* ( 1. *-Y 1 +GA0A2* ( 0 f L T A ~C -M? ) I ^ 0 0 II G A * D V (A . 1 ) 

A(jj*2 */(Y1+SAM/\2*(0EI, TA”C OK) I 

o V ( 9 , 1 ) = A U A 2 ( X 1 - I * - G A A 2 « ( 0 F L T A - c n R > ) ♦ 0 M E G A » D V ( S t 1 ) 

0 V ( ) 0 t 1 ) ♦ 

On \Z J= It Moots 

o v’ ( H » 1 J = 0 V ( y I P + A u X 1 » 1, ^ II A 2 ♦ A 0 7 I J ) * t r V ( 2 « 0 ) * I 0 r G A 0 M t G ft - 

X HU ( J ) ♦ HU ( j ) > - U N 2 { vJ ) •» 0 V ( 4 ♦ . I n 

Uv { 9 , I ) =1>V < V , 1 i -aUX2«GaMa:>*aM2 ( j > « I V ( 1 . j ) • < 0NrGA*0HEGft» 

V H U I U ^ * i l 0 ( J ) ) - U N 2 { J ) !k n V I 3 » . I ) ) 

12 CoHTlMUt 

1 F ( T " tE « M 1 RK T Oi^lJ 

V ( y , 1 ) = 0 V ( 0 , n 

V { V , 1 ) == ['» V ( 9 , n 

V { 1 0 , 1 ) = n V ( UJ » I ) 

ReTupo 


) 


2 


3 




Suf.ROUTiNK PUOY 

This sour hut ink pi. OtS uOTATlOfj AMGLk vs u 
C n f'5 M 0 i 1 / 0 n II / G A M M I P K X t P k’ Y » X •'’i 0 * Y ,\j 0 , M U 
C o H 0 i-J / L E M / 1„ /i M » S 1 0 
C 0 hMON / T HR Ei ,:/ modes 

t 0 'i H 0 M / K I V K / M K » H j C A 
C 0 M M 0 N / 5 I X / M A X P , i-' R E C 

0 0 U B I- E P I' K E ! !> i 0 N S I n 
KpAl, MAX »U AM, link 

U I H E N S I 0 ' J S A V E I I U 0 ) , L 1 14 E ( 1 1 □ 1 .APIS) . A N I G ) 

Data b l a m k * s t a k , ih) t / i h > i h * ♦ I u * / 

1 F ( N • M E. • 1 ) GO ^0 2 

= ( OK +srj J / I f.u) 

L>o 1 J I - 1 , 1 on 

S A V K I ‘H ) - rj ♦ 

-) 1 = i? 

Oa/=:CU . . . 

4 1 =0 

5 I 0 R ~ S I 0 

I F ( H/N I 01) i • UK .0 ) GO To 3 


J I = J ] + 1 


^AVE^'^^ )“CA 
J F ( A n s ( c A ) . r; T . I i A X ) 
i F ni * 0£ • \'K ) PET UE4M 
vx R I T K I f. , q ) 

B 0 H M A T ( 1 t-l i , » P U-* T 0 r 
'*■ R I T F. ( *» » S ) L A l l * S I (> f\ 
F ormat t i h , * i. amri^ a = 
X i H , * P K X . F • 4 , / 1 H 
X IH ,*PKKr.“*i 


M /\ X ™ A B S ( C A ) 

I'J U T A T I 0 N A 11 G I . K 1 N i>C G i? K i: S \> S' H F Q P * ) 

. G A M '' J P K X » [* K Y , P P ;• C , M A X P , M 0 OE S 

• . K 6 • 0 » / I H , • s I - ; E p 0 = ’ . F 4 . 2 , / 1 H f * G a M A - 

, f !•' K Y = * , F A o 4 / 

12/)H , • !|AxP- » » 1 3/ I H , MinoK5= • » 1 2//// ) 


.F70,/ 
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M 1 s H A X / 0 , 

L> 0 /> H a ! , G 

A N t I M a - A ) ♦ ( 6 0 • 10 . * I ! ) 

^ A p ( « I 1 I s *. ,u.H 1 i ) 

I \ ^Q 

l'»R 1 TE «■ 6 ,7 ) A0,JiliAP 

^ E 0 ! ? ! 1 A T ( I H . 1 X ^ b U‘ 6 * ;? . q X ) , 3 )( , I I t Z X » r, ( S X f P Fj • 2 ) } 

0 0 fl J I - 1 , 1 J 0 
L I N p ( J i ) = fi L A 0 K 

^ J F ( ( ^J t + 'I ) / 1 0 « I Op f: u» . ( j I + q 1 ) L I ,i F. ( J 1 ) ta s T A R 
i‘/ R I T P. ( 6 ♦ 9 ) L. I N t: 

F n M A T ( 1 H I 3 2 X * ) I 0 A 1 > 

• 0 0 U) J 3 = 1 t 11 0 

Jn L imE ( ••M ) --.star 

Wr t TE * 4 , 9 ] U I IJE 

0 0 I I J i - 1 » 1 1 0 
1) L I me { J 1) =UL ANK 

On K I =1,1 00 
J I = 5 A V E ( K; i ) / A 1 •*• 5 A • b 

L 1 M F { b 6 ) s s T A I’J 

1 F { K 3 / 1 G * 1 0 ♦ N E - K 1 ) <:;0 T 0 12 

L I E I B b ) =: s T A R 

12 L I rjE I J 1 } -'OUT 

‘■V R I T f t R) LI N K 

1 F f K 1 2 1 OM CteNE ->K I 1 GO TO 1 S 

1 F ( j 1 i. GE , BO . ANO , J 1 • LF . 5"^ ) GO TO I 5 

0 1 = rn * K 1 

'‘V F? I T I- 1 A t I q ) Ml 
1 P 0 R M >\ T ( I H + t A » X > T B ) 

lb L f E ^ J 1 1 aFJi. A NK 

13 U I F ( b b ) - L A 0 K 

K F T U R f'i 

Lnd . 

' L /. T ] 0 i i ; M 0 0 ; A G N 0 B T I C 5 * 
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CHAPTER 6 


Other Topics 

The present chapter contains a short note on the use of the sta- 
IjiXity charts in deployment dynamics, and two bibliographical reports 
on passive nutation damping devices. 


6.1 Stability charts and deployment dynamics 

6.11 Dynamic parameters during deployment 

A deployment phase such as the one for IMP-I may be summarized as 

follows, if H, A-, Kot Wg designate the angular momentum, length of 

X . — or + 

booms, non-dimensional radius of the hub (■^)» spin rate 


STATE 

** 1 " 

(for t<0) 


STATE 

(for 

0<t<tj) 


”k" orbital periods 



We define a "state" as a set of values H. 


If non-dimensional variables 


H 

o 


def 


H(£-0) 



m 

So 


are used, let 
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V _ HCanv state. 
h(any state)^|^ 


h(£=0) = 1 


0} (0 ) = (D 

s def so 


(1 + 

I ,0) 
zh So 


2 1 "s 

= (1 + Ar(5 + c + ^))-^ 

o o J U) 

S 


__ p£^ 


A = + 3 


I , = moment of inertia of central hub about "z". 

z,h 


p = linear density of boom 


In view of these definitions, an extension maneuver at t corresponds to 


h(t +0) = h(t - O) 

£(t + 0) = UX - 0)+ AZ 


0) (t - 0) 
s 


(1 + 4rA)t-o 
(1 + 4rA)t+o 


in which A^ is specified. 


A respin maneuver at t will give 


2 1 0 ) (t+o)-w (t- 0 ) 

h(t + o) = h(t - o) + (1 + + Cq + 3 ^^t+ 0 ~" 


2 

(1 + 4r(c; + ?. + 5))^+o = + «. + 5 ))t-( 


in which 6m (t) = w (t + 0) - m (t -0) is specified, 
s s s 
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For a satellite of given hub (xq, I , specified) 


Co - ^ 

zh 


zh 

1 / ^ 

£_ .V3 ^cP 1 




= i 1 


£ 


1/ 3 


( 2 ) 


with S a fixed non-dimensional number def 


X P 


1/3 


1/3 


zh 


Now, substituting (2) for Co Equation (1) 


u = (u [1 + 4r(s r + S 

S Sq 


2 - 2/3 _ - 1/3 _ 

f + 3 )] 'h 


If £ is specified in any state. 


P ^ 


z,h 

can be computed. 

To that state there corresponds an Etkin's number 


„4 2 4 -f 1/3 - ^/2 + h)]”^h' 

X = x,[i + 4(sr + sr 3 

El S EJ;^ 


The quantity is specified for a given design. Let R be the non-dimensional 

EX 

quantity 


R ^ 

i 


{-) 

Wm ^ 




£I 


H £ i*® 


Then 


l. i C \ r"‘ ,r 




(2) and (3) thus give C^, ^ during the "states" of deployment as functions of 
boom's length and angular momentum. In these relations, end S are fixed 

for any given design. 
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6.1.2 Stability during deployment 

The determination of the stability during deployment will thus 
proceed as follows: 

a) R/m in computed (a fixed quantity), then R for co^ given. 

b) S, a fixed quantity, is computed; 

given the state ^o c^nd H for some t: 

c) compute h; 

d) compute F ( ); 

e) compute A (Co) 3 

using the relevant formulae for either respin or 
extension maneuver 

f) -compute X from (3); 

g) determine the stability of the corresponding (K^, T) point 

on the stability chart corresponding to the computed 
values of X and Co » using program FLEXAT of Chapter 5. 



6.2 A SURVEY OF PASSIVE 
NUTATION DAMPING TECHNIQUES 


Prepared by 


Wllliean 0. Keks 



1 


I . I n. t r o d u c t i o n ^ 

In this paper, several methods of passive nu- 
tation daiV^ping ai^o survoyod. In a review of rigid 
■body dynamics, conditions of stability are pre&icnted. 
Bal 1 , p e ]•] d ul uvu , a n d f 1 ui d damp er s a,r e s ur vc y c d , amo ng 
othex's, along v/itli effects of magnetic and gravita-* 
tiona'j. torxiucs and, stx'uctural hyste>rf;sis energy dis- 
sipation, Fincilly, a few active and semipassive 
system.8 are mentioned in the way of comparison. 



A,B,C 

X,Y,Z 

in 

.H 

V 

r 

0 


r 


62 J, 


T 


ji 


i,3,k 


m 

s 


Mo Til 0 11 1 s of i n e r t i a ab 6'u t x , y , ax e g 

Products of inertia for xy, xz, yz planes 

Body-fixed axes, z along spin axis 

I n e r t j; a 1 ax e s , Z al o n g 11 

Total anp;ular velocity 

C 0 uip o n c n t s o f al o ng x , y , z ax e s 

d( )/dt 

Angular momentum 

Gene v a 1 i z e d coo r d i n a. t e s 

Moment in direction of q. 

X 

Euler’s angles 
Precession rate 
Spin rate 
Nutation angle. 

Magnetic or structural hysteresis factor 
(p"^‘ r^') = component of in xy i>lane 

(C ^ A)/A 

r ~ fol’cing frequency 

Unit vectors aAong x,y,% 

?\lass of main body 
Dajiiper mass 
Piadius of gyration 


Other symbols are defined througbout the 
text as needed/ * ;■ 
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II. R e V i o w o f J L ;l g i. d 1 ’ o d y Dy n arn i c s [3 

A, Uclinxtions 
1. Eulex’^s Angles 


If X,Y,Z is fixer] in space and x^y,?! is tlie body 
fixed system, v/c dex?iiic the Eulox'' Angl.es and <f in 

Fig. H--1 . The spin axis is along z, and: 

~ precession I'atc 
Q = nutation angle 
^ ~ spin x:*ate 

The unit vectors lie along x,yjZ. 

\Ve have : 


x' 

Y 

f< 

Z, 

. 


cos ^ ’“Sin <p cos 0 sin yO 
sin^*{-Gin cj> cos O cos ) 
(sin O sin Cj>^ 


(-sin cosji' “Sjn ^ cos cos 
(“■sin^^ sinj;^ j-cos cos ^ cos (/') 
( s i 11 6? c o s p ) , 


(sin O sill f ) 
(-sin O cos ) 
(cos 0) 



' " 


X 




j 






2. Angular Velocity 

If (?) =: OJ /i -f cu *5". -H c/y'k is the total angular vel- 
X y z 

oeity of the body, then: 



xl 


(sin <9 sin (f ) 

(0) 

( cos {?? ) 




( sill cos p ) 

(0) 

(-sin 

— 

c.(t> r. 


' “ ( cos ) 

" ( i ) 

(u) 

Note 

: th, 

at 

r is not the 

spin 

rate , 


In most eases 
of mass is ignored. 


the linear velocity of the center 
for d a !'.!p e r an o,l y s i s . 



d 


3 , An K ul ai* M o iri o n t urn 


II = 


A 

:i;iA 

-13 J) 

-iii.E 

A A /• 

,iiD 

• JJB 

A A 

-jlcF 

A A ^ 

klE" 

A , 

^1 1' 

A A 

kkC 


cv 


Tb. e a]) o V e :i. s the angii], r^r mo va e n t u.iii j. n x , y , o r 
•most cases, v/e can ignore the external torques pro due* 
ed by electromagnetic f id, els and gravitational gracli- 
.ents. Thus Tl is constant in j.nertial fixitice (X,Y,Z), 
and tlius avc can align the Z axis along H. If x,y,z 
are aligned along the principle axes of the body, 
D;ir,F'= 0, ami: H ^ Apr r Bqj -i- Crk 

where p,q,r are the for al igmnent ' v/ith the 

y /j ; . 

pr i n o ip 1 o ax e s » 

4. Kinetic Ener.gv 


The kinetic energy of the body is: 

o p P 

rp _ iif \ /n, "* f .T^vj-1 • r* //} > ^ - D 


T rr yo{A(^ -*r ) 

X y '/j 

an d f 0 1 ’ t ] 1 0 p x' i n c i p 1 c ax e s : 

T l4(Ap^ + Bq^ + Cr^) 


E60 w-' , Fcg; oj 
X ’A y 55 


5 . F. ul or * s l.i] q uat ion: 


Hex’e Euler's Equations ajt^e presoyited only for a 


pr inciplc 

air X s 

X , y , 55 : 




Ap -!• 

qr(C--B) 



= 

Bq X- 

pr(A-C) 



k. = 

Cr X- 

pq('B~A) 



vdiex'o E^, 

I i ^ , a 

nd are the 

Kf 

c >r. 1 0 X' Vx c J. m o rn o n t s 

about the 

coi-r ospomling 

j j T i n c i p I e ax e s 

; hero they will 

usually 

be 7.ero, 







r o j n s o t ] !i 1 1 i p s o I cl 


l^cr a rigid body, T constant, c^nd thus: 

-- 

6?H-I „ 22 „ o 

n r' “ u " 

This iimst be the coiirponcnt of £c> along II and Ze 
sides of t)ie energy relationship are divided by 


If both 
T, we 


get: 


2 

V 

2T/A 




2 




r 


2 


2T/B ' 2T/C 


This Is the Po inset ellipsoid. If a plane is placed per 

, 'i 

pendicvilar to H a cUstance Q frein the center of tJiis 


ellipsoid, see the Poinsot ellipsoid rolls on the plane 
(called the invariant plane), without slipping. The cou'- 
tact point is the tip of Cd (J-'ig.lI-2) . The curve traced 
out by the contact point on the plane :i s the herpolliocle , 
and that on the ellipsoid is the polhode. 

B o dy a n d S p a e Cones ( /u:, i s y n ]m e t r 1 c Bo dy ) 


From the above we see that ^ sweeps out a surface 
in both the x,y,K and X,Y,Z fravnos^ If 6 0 and we 

have an axisyirtmetr ic body (A-D) then these arc both 


r -j. gh t c i I ' 0 ul ar o o n c s c , Fr o m t h e r c 1 at i o ns b e t we on p , q , r 
and ^ substituted Into the Euler moiDcnt equations, 

WG li.avc: : 

™ ... c y 

/ CK-Z)qos 0 


( a ) C rA : ('> a n d.' ' ^ ar 'e ' op p o s it e i h " s ig n , ah d t h i s i s 


k n 0 \7n a s r c t r o gr ad e ]:j r e c e s s ion. 

( b ) fX A : an <\ ^ h a. v e t h e s a m e s i g n , 


and this is IvnoTm 


as direct or posxgrade preocs.sion. 



^ — I/ines or nodes 

Euler's angles. 



plan 


[“?■: rojnsot ellipsoid. 


7 . 



(a) Retrograde process Jon; C>A. 


n 



(b) Direct precession; C<A, 

-3: Precession o.C body cone rolling’ on 

space cono.cc^ is along the line of contract 



8 


The body ooug rolling on the .space cone for each 
of tijcse cases is illustratecj :ln Fig« 1 1-3. The angle 
bo tree n '^p ami is )( : 

tan >r - (p^-j-q^ ) 'i^/r 

v.li ere W i s tli e c o ir.p o n e n t o f ^ 3. y i n g i n. t h e zc^y pi a n e , 
The ajigle between and If is6> : 

. tan 6> == {A/C)/{(<'^,^yr) 

By .subst :i tut ing _a u: [(C-A)/a3^' ~ Ar into the 
liluler equations, v/e have: 

p + J1 q " 0 — --- p rdtq 

q - -Up = 0 
2 

Thus p + ^ p - 0 


and 


P 

q = P 


pQ cos-d,t -i- (Pq/h.) sinllt 


sin -At “• (p^/a) cos~at 

^ ^ p . p 1 / 

These last imply tha.t<^,^ ^ (p -i-fi rotates about the z 
axis at therateD-. 


By using a complex analysis, Ames and Murnaghah 
show tha,t I'll : 

CA» f id ~ scfl'® 

T [AfA~"c'a ® 

8. A Note on Unsymmetr ica.1 Bodies 


' The relations for ^ are. 'given . by Thomson for the 

c a. s e“A^lT^ C aii d ' ] I 2 T H . a b d dy s p inn iiTg ’ ah out 1 IT^' ax. i s 

of least inertia [33f0 : 


- 

2Tcl 

a{a 

C)J 
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>- • 2T cV^^ 
]j(ir~ c'j-j 

'[c"(A ‘. C)J 


sn f(t t^) 


dn f(t t^) 


'ivhero f = 


C)(2TA •- 


and the modulus ot the cllxptla fiin-ctioiis is: 

.), r (A. rs ) (id' - 2TC )1</. 

_{D - C) {2'J’A - H^)j 

This results in spin about the z axis with a super j.inpo sod 

wobble, with a. (9 and O . : 

max min 

cos^Pj.,^ = C(2TB - ld)/(B - chd' 

= C(2TA -• )d)/(A - C)I:I^ 

hlXil 


B* Miscellaneous Concepts 
■ 1. Stability of a Rigid Body 

For. a rigid body, T is a constant. If we let the 
initial condition bo: 

p r.:: -I 6 

q,r small 

Vihe3;e <S- is small, v/e can differentiate the Euler equa- 
tions and substitute for p,q,r and p,q,r. Then: 
a Pj^qCA - B)(A - C)/J3C 0 

r + p^^r(A " B)(A C)/BC - 0 

These arc stable only if (A - D) and (A - C) are of the 
same sign. OBuis they arc unstable only if A is the in- 
termediate rotationaJ. inertia. 

2. Energy aaid St ability 

In a real space Ci a.i. t , xhf;!rc is always an energY loss 



:io 


cl u e t. o i’ 1 c xur e of n o i ) r :i. g i cl p a r t s , m a g ji c t i c hy store s i s , 
etc. Thus V;G have T-^O, 

For an axfsyminetr ic hody, ive have: 

2T = H- Cr^ 

I] - A. C r 

Si'noe Cr ~ IT cos <9 : 

n^' - 2TA = cos% II^(C - A)/C 

or T = 11^' [l - cos% (C - A)/c] /2A 
Since there are no ex'ternal torcjiies, II Is constant j and 

= [h^(C '= A)/Ac'](siD ^ cosO )S 

= (H^A/C) ( sin d" COG O- ) S 

Thus, for decreasing T, ^decreases only if C>A, 
and the satellite is, spinning avboiit its axis of muximiirti 
inertia. This is the stable condition. B'or a prolate 
body, there must be cin energy input for stability, which 
i r,ip 1 i e s an a c t i v e n u t a t ion c o n 1 1 ' o 1 . 


The change In energy recjulred to stabilize a pro- 
cessing hody can easily be found. The desired 'energy 
state is: 


T^ - ^Cj;^ 


2 


where the subscript f denotes final condition. Since If 
is constant f3b] : 

= A>/- . = cb/- = h/- 

Tlien = (h/0)w,y + 

Tlrua at = |T - T^j = j>Al(l „ A/C)t^{ 
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For an oblate body (A-^C), this is tlie precessi onal en- 
ergy, the amount to be remove cl; for a prolate body, it 
is the amount to be added. . 


Also = (j:" tan 6'^ = (C/A) r tan, 


III, Pels s 1 ve Da.mp e r s 


Unless stated otherv/ise , the satellite will be 
assumed axisyimnetr ic about the z (spin) axis, A. = ,B 
and oblate (A-^^C) for the beloiv, 

A, Ball- type [4, 24, 36^. 

1 * Mounted i n the Mer i d i o.n Plan e 

This type was first used in Telstar and later in 
ESPO II; These consist of a ball alIov;ed to roll in- 
side a circular cross section cui'ved tube vdiicli is 
filled vrith a gas. Two a,re used, diametrically opposed, 
to maintain symmetry, and mounted in a plane through the 
spin axis. Energy dissapoAion comes about through vis- 
cous friction betv/een the ball and ga,s , rolling friction. 
bet\veen tlie ball and tube wall, and coll.isioji of the 
ball \vith the tube end, the latter only at large nuta- 
tion angles. . - 

Such a system is shovm in Fig, III-l, According 
to Yu, Die rotational motion of the ball (of radius a) 
is given: 

(2/5)ma^(^R/r) ha - N . 

ivh ere k is the friction force at the c o n t ac t point , and 
N t h e 1 ' 0 1 1 :i n g f r i c t i o n t o r c i u e . N i s a] ) p r o x :l i n a t e 1 y ctn 
order of magnitude smaller than the viscous ternic 
lecting M and assuming B small, Die motion of the ball 
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is describofi by: 

5 + (ric/7ni)£2 + (5brV7!;)oi^''~ ( 5b/7R)cosJX t 
where c is the coerficient of viscous friction. The 
time average rate of energy dissipation is, for viscous 
frictioji: 

dT^/clt - T^/1/2 y/ = 

where = &(,! ~ ./L‘'/r^) [(1 -i- 4n^l.^'/P^] 

n ~ 5c/l4m 

and P -- (5br"^/7R) is the )iatural frequency, the 
square root of the coefficient. 


We end up with an expotential damping: 


and 


5cln -.M^/P^)\,. 4n^"/P^ 

7nraK^;t ( >, 1)^(1 - A 


If rolling friction dissipation is included: 
d^g/clt = SFRcy^qI/I//-;/- 

v/herc F is the rolling friction, paid: 

(Pq 

'Y' ^ Tin [a ■!• Pq//^-)/(1 -y 0-Q//-n)] 

Numerical computations show that is substantially less 
than one degree. Thus the viscous-only results. can be 
used if {9^ is somewhat greater than one degree. 


The darapi ng can be greatly educed bj^ designing 

a .r.esoiitant system, making P ~f>, , Then: 


31 


r c 


3b/7;^^- 
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The time co:.istant Is then: 

"^rer "* ^ *■' ~A )^ 

A resonant damper coaid not he used in Telstcir be^ 
cause A 'v-’as close to zero and room had to be made for 
an electronics package, preventing a smcill value' oi‘ b„ 

It is possible to conceive ol dampers using 
straight tubes or tubes concave outv/ard. It is easily 
seen, hov/ever, that the equilibrium iiosition for the 
ball during nutation v/ould be at the ends of the tubes, 
and the final spin axis v/oi.ild not coincide with that 
0 f t h 0 s a t e 1 1 i t e vj s t h out t h e h a 3. 1 s <> 

The pa,rametcrs for Tel star v/erc A/C - .95, 

^ 20-^180 rpm, R l5 ft, in t-_ 0.0021 slug, a = 0.24-2 in 
(tungsten for its largo density), c - 0.00193 Ib-sec/ft 
(neon for its high viscosity). The theoi’et ieal damping 
time v.-as caiculated to be a maximum of about three min- 
utes. 


Note that a gas of low viscosity should be used for 

a tuned (resonan.t) damper, £is n, proportional to c, ap- 

liears in the numcra.tor of the exin-esslon for 'Y' . 

res 


The problems in this £inal 3 ''sls are 
s u m e d s in a 1 1 0 a n d 1 3. 1 'i e ar i z a t J. o n o f t h o 


duo to the as- 
equations o G.T. 


Kossyk devised a ground test of a model supported at 
its center of gr.avity \Yhich shov/ed that the exTjerimen 


tal ''y was £ibout four times that calculated using the 
111 o a n V a'j. u e of the t r a n s v c r sc i n c r t i a i ii o m e n t s , an d n j. n o 


1 3 hi c s t ; 5. t using 1 1 J G li ! i, n i mum v.al u o . T ak j, n g t l i e s e f a c 



■J.5 


1 0 3 ' s :l n to a e c o im t , t h e f o r T e 1 n: t ar v/a. s c a 1 c ul a t c d to 

be no jiiorc than thirty rtiinntec, 

2. Mounted in a Plane Parallel to the Equatorial 

Two of this type were mounted in Ihh’Ij and one in 
the IIEOS spacecraft, which also used a liquid damper* 

If h is the distance from the damper pla.ne to 
the center of gravity, Routh criteria applied to the 
Euler equations indicate that b/ll ^ 1 - inli '/A is neces- 
sary for stability. Also, optimum damping (minimum T) 

. is given by a viscous friction coefficient of: 

- ri'U^r[l4n;lP(7' + l)VsA^ 

This results in: 

"^opt = 0/r)[5GAA/inli^(A.>. 

Experimental results agree v/ell with the theorcti-- 
cal. For. two dampeurs and ^ 0.61, h = 0.15m, Pt - 0.2m, 

r = 0.2 rcid/sec, anci 250 gm give a maxinmni of 120 see 

for reasonable^. The e.xpei’iinental result was 130 sec. 
•V/ith all part'iinetcr s equal, the efficiency ratio of the 
equatorial to meridian damper is [(1 -I >1 ) / ( 1 ^ ^ , 

B. TEAM Damper [24,25] 

The TEAM danipcu', used in Tiros, is essentially tlie 
same in concept as the mor id ian-mo unted ball damper. A 
small mass fitted with roJ.lcrs is allo^ved to run along 
a curved monorail (Fig. IlI-'2). The difference lies in 
that 1 1) er e is no f 1 ui d i n vo 1 v e d , so o nl y r o 1 1 i ng f r 1 c- 

t:ion exists. From the ball damper analysis, it can be 








17 


«ecn that this would behave v.'oll only a. t ssnall « 


F o r T j. r o 8 , t h e d. amp e r i iict. s s v/a s al.) o ut 0 « 0 0 1 of t li e 
total s a t e 1 1 :l t e rn a, s s , an d a s s ur e d. a 6^ of less t h aai 0 ^ 5 
d e gr 0 e s . The 1 1 me t o d a,mx) fro m 2,5 to 0,5 d e gr e e s v>a s 
about one minute. It v/as chosen becau8e tests shov/ed 
that the tube radius H of the ball damper ivould be 
g r e a t. er t h a.n the t ci- c Jc X' ad i u s of T ExlM * A1 s o ^ it v/a s 
found that the ball damixex,- I'equxi'ed an A/C not less 
than 0,375), Y/hex^e A/C for Tiros was 1,-1 5 

( - o,3i) , 


C, Peixciulum Dampei" 

1 , Six i IX Alii s P i VO t c d 


The motion .fo'r a. satellite with a ixenxluluri} pivoted 
on the spin axis, and moving in a plane perpendicular 
to the axis was described by Cai'’twx'ight , Mas.si.ngill , 
and TTnoblood [g] , The driving fx'equency of the pendu- 


lum is the fx'enticncy of the acceleration due to mxtation, 
Jl = 1 ' , Without fx-iet lon the jx^nrlxilum xvould oscillate 

in synchronism opposite at-Q , as in Fig,IlI™3a, ITov/- 
evex’, if the pivot exerts a f i-ict iona,l toxuque, the pen- 
dulum lags behind this position by an angle S (Fig. III*‘-31x ) 


The resulting iorqxie on the axisymn-0t3:'ic main body caixc- 
es the damping. As this lag angle ixici'cases, so does 
the damping, pi-'ocmeing the convex portion of Fig, I II- 3b, 
and called the '‘rmtation s3’nchronous*' mode. 


V.'h e n x:* e ac h. e s 

is n.o long ex.' ijx Fxync 


90 degrees, ho xv ever, the pendulum 
XV i t h , b u t i s dr i. v o n t o xvxxr d s y n- 


elironism with i'. TJiis xs a. dccreasing-x'ato. decay with 



la 


a r, up 0 r i. i: ■ p o .n c cT c o n v o r o r? t o n o :i. 1 1 rat i o n . 1 1 ayo ti 1 d b o d o - 

sir able to iiva'cc t]ie ti-ans it ion ..botraen tbc tr/o stkjOcs at 
a s s r, 1 al 1 ra 0 a s p o s s i b 1 g . 

. • J r t h c m a s s is a s s uia e d s a 11 so 1v i i a. t 1 1 1 e r o t ra t e s 
precisely at tlie nutatio3'i luate ' ( ?\ -’r . 1 )r - XL r in iner- 
tial space, a,nd -Cl ^ye Iirwe: 

O ~ (niih/A)(/^ -f 1 )r sipo.^ 

-i- ( c / ir: )<k ( h/.i ) ( "A 1 ) r s i )i o4 - - ( c / r 

\v]3ere is tbs betT/een the x a:;: is and dq.,, a.ssunied 

a.ppro:r:ii:iately equal to ^ , Also, is the frJctioji oo~ 
eiXicient ot the relative velocity betv.'eeji, the }pendulu:in 
and inain body. 


Conip liter finaXysis has shov/n that t)ie terir: can be 
neglected. To find the time • and nutation angle at tran- 
sition betivocn modes, i7e setl^ ~ Tf/2 . aiid integratt> tlie 
above equations. Tiius: 

t. = [ ( ) /cy ^] ( C/2 >) 

N umer i e al i n t o gr a,t Ion of ex ac t e quat 1 o n s show that t h e 
first equation overestimates 0^_ by as iimch s,s a factor 
of 2, andthe second umicrcst imates t..,, by as I’-iueh as a 
factor of 2„ Also, for these equa^tions to be valid, oL 
iiius t “be no .ai’ "Z ei' o ■ at' t j, me — t 0 j ■ thus : : — ; — ^ — 


mf 




C^^^^(A/C)/(1 A/C) 


is a necessary condition for tlrcir validity 


13 e c exu s e of i. t s no ji s y mm o t i^y , ■ 1 1 1 er e w i 3. 1 b e a vS iiiai. 1 
f i n al n u t a 1 1 o n an g 1 c xdi c n on 3. y o n, c d amp e r is u 53 c d : 




(mi?h/C) ( A ■!- 1)/'A 
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If c i s 1 ar g e , t de or e as os but (9.v. t n cr eas c s , 

If is an increasing fuucticn of velocity, there 
v.'ill be strong damping at the beginning. As the rel- 
ative velocity decreases, so does c^^, and the damper- 


inaiJi body ' system is decoupled enougli to dcla,y tran- 


sit ioj^. 


Another improvement ivould be to use two pendulums 
of diffei'cnt radii. Experimental results show these 
to act independently, the long one damping quickly at 
large $ (Pig . III-3d) , the short at small 6^. 

2, Pivoted Away From the Spin Axis 

The problem of a pendulum moving in a plane per- 
pendicular to the spin axis and pivoted at a point 
av;ay from the axis have been studied by Ilaseltine (^16, 
17] and Newkirk, Haseltino, and Pratt f23] . 


If is the rotation required' to reach a point on 
the body, tlie kinetic ene3rgy of the system is [23j: 

T = 'h [cf^- + 

where § r + • 


m Mm/(M m) 


a n cl - i* — is t h e - d i s t a n e e 

m 

the d a ! • : p c r m a s s . 


. 1 . j : \jin~ 


til e • ■ c e n te j. ** of ■“ gr a vi. t V" t o 


Using a set of modified Lagrangian equations; 
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Centr Ifugal 
force 


Offset 



( a ) 0 f f - d . e s i g e o q u i 1 r 1) v 1 u rii 
f o i” t V/ o of f sot p e 1 i cl 11 1 uni 



\ 



( b ) 0 ;l f - d c s :l, g n e a ii j. 1 J. b r i u cii f o r for r 
pendulums . 


Fi 


JIT-S 


Offset pendulum dampers. 
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in which all tlie h are c:NCGpt: 

Q J. ♦»•• 

L -C ' 

■ V i . 

TliG r 0 1 1 owing e q iiat i o n s x* o s ul t : 

I -^(c^^/a)(.? -- r) 

r r. >AOp(,^ r)/(AC ) 

H- p/(AC D^) .-[- pq 

q -^pr ^!- CC^p/(AC D^) 

+ DC (.^ xO/(AC ~ D^‘) 

p ~ (A - C)qr/(A c) 

-I- C^q/(A + C) - D(r^ - . q^)/(A + C) 

^ - 2 
■where A ~ A my 

C - fnh^ 

D - ivjby 

and y is tJie y coox'dinate of the mass. No shiall aiigle 
afifuiniptxons have heen made. If however, O and m are 
small, J ooiistaiii, and other limiting assTimptions are 
made'; 

S ~C^(S 5 )/c 

+ (CD/2AC)|(Uc‘”^^ + 

U i(C/A)lU = ~(D/A)S^e^^ 

v/hei'c S =: p -X 


U sin O (cos^ -I- i sin^) 
Thus [u( ~ sin <9 O, 
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Three dlxTcrcnt aolutionc^ were tried for this set of 
equations: 

(a) The stable solution in which the damper does not 
rotate r el active to the main body. Then |U/ ~ ^0 ~ 

constant. 


(b) ”Slov7 damping" in vdiioh the damper has a small 
oscillation a.bout a fixed point on the body, resulting 


in: 


0 - 


r - 


6 J.)C- 


■t/r 


2(C - A)AC^' [(c^ycp^i- (C ~ A),^jgVc''1 




CD ?o 


P 


(c) "Fast damr)ing" in wbicb, the damper rotates at 
the nutation' froqueney Cr/A, This solirtion is good 
only when 0 is not small: 

0 "^- - oj' ~[2Cp<C - A)/CA ~ 2?/#C#/A^]t 


The advantage of offsetting the pivot point from 
the axis is that it v/ould appear that the pendulvun 
will align itself radially outward from the i>ivot* 


Then a couiiterniass corild be 
um ^ s i i i ^ -t o p e s e:c v e . t h e - 


mounted from the equilibri-- 
s y in?n e t r y - o f . - 1 h e - s at e 1 1 1 1 c , 


with no residual wobble. An alternate is to employ t^v'o 
d 1 am 0 1 r i c al 1 y 0 pp o s e d p e n du 1 urn s ( F i g , 1 1 1 . 

r 



F or a ])ea cl u 1 u rn o f c s o t a tl j. s t im c c b ao rl o f ar ni 
long t h f- , t li. a f r e r? u en c y J. s : 


277 


v"r* I 


hC 

/A 


t£LnO 


For resonance: 


I ^ [h + (A + DhOj/Z'-y b/A^ 

as suiting ” i"* for small ^ , If ?* is only slightly 
greater than. zero, V can be largo, A solution is to 
uso a pendulum of x'adius of gyration s. Then: 

^ = [b -I- (A + l)h<5?/A^(l + 


ITaseltine [igI has shown that, when: 

(.//b)[mh^/(C - A) - % 

The angle between the two dampers will not be 180 de- 
grees in the steady sta,te (Pige IIl-'-5a) . Then: 

cos ^ " lesser of one oi* b(C ‘-.•A)/2.nfi ^ and the 

apparent wobble aingle is approximately ( 2riih sin )/(0 - A) 


-Ilaseltino also studied the motion with four identi- 
cal pendulums moiuited 90 degrees apart. Again, experi- 
mental ' result s showed possible equilibrium positions 
X'esi.’lting in <x x'esiduaul wobble (.Pig . III-5b ) , 
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D. I.i liquid Da rnp 01 ’ n 
1. Spin, Concentric 

The use o.r an annuius pan'tial’Jy filled \7ith a, dense, 
high Yisco.sit 5 ^ lieprid, usually niercury, has proven very 
popular; It vaas first used in Syncoi?} and the Explorer 
series [24-]., The basic tJjcory vaas laid out by Carrier 
and. tiiles [s] for laniilar flov/. The equations of motion 
for the body are similar to those for the pcndulian, since 
both systems are circularly constrained-. The dimensions 
of the system a.re given in Fig,' III-6, For small ^ , it 
v/as assumed tha^t the liquid wris in contact v/itii the en- 
tire outer surface of tJie annuluSo The r*ate.at v.iiich 
enei'g 3 ^ is dissipated throughout the fluid is, if is 
the density: 

T X v)^dV 

Vvhere v is the fluid velocity, -v is the kinematic vis- 
cosity, and dV is a differential element voliune. As*^ 
Sliming the irro tationtil component of velocity cannot 
con t r i b u t G to t h e into gr al : 

T n..,^|G| //yA| ^ 

vdiere n.. - (1 

, 9 V, 

Cr = {xfir/y)- 
A = (n^. '■ n)G h- (n... + n) 
n =. 1 + 2 

T h X s r e s u 1 1 s J. n a time c o n s t, an t o 1 d e c a,y f o r 0 of; 

r := rA/oid)df ]f\ 

~ n)^|G|^ + - n^)j(;| -i (n.,.. + n)^J 

3K'V( A 4 l)r,,.d'|G| 


mirt r 



2G 


ThlB is at a minimum in the nciriibborhoocl of n.,,. == n* 

Then:, 

-- 1 ) iGl 

and , if a..,/n<-l: ■ 

nhVd +> 1 )^ 1 / 2 ’^"! Glj 

is the resoiiajit condition. The vai'iatlon in thickjiess of 
a.,:, has been, a s s n.nic d s 111 al 1 , 


For 1 ar g e ^ , t ii e fluid c 0 ni}:) X 0 1 e 1 y f i 1 1 s the cr o s s 
section of the anuuXun^ over .ah angle @ (Fig , ) . 

The energy dissix>n,t ion is then: 





The time conEitant for large is: 

r = K\!y\B^ 78 \f>\ (Ip)-i' (fi -I- 


For R = 10 cm» h - 10 cm, d = 0.25 cm, a,„. r-. 0^05 cm, 

A ~ 1.3 kg ~ m^, }\ i"~ 1/3, ay ~ l2 rad/sec, ^ - l3.6 gni/cc, 

""•32 ' 

and ^ " lO cm "/sec give a dfimping time of 14 sec for 

sma,ll , If a resoUvant dajmper were designed, a,, Vvould 

be 0,637 cm and T - 0.00044 sec. 


The large 6 result for the above i:>arameters and 
Oq - 1/6, 0 ~ and (a + d) “ V 2 cm givoE’ T " 200 sec. 
However, the Ucjniolcls n limber is past critical for 
t h c s e j- • a n cl t It e - i n c r e I'l s e d"~ f il c t i o i i ■■'wg 1 1 1 0. ' r c d u 0 e~ 0 "to 
about 70 see.. 


Tlio above would indicate that it would be desirable 
to d c s ;i. g n t h 0 d £unp c r f o i' r c s 0 n an c e . .ITo v;’ ever, a s t u dy 
by Fltx gibbon ajid Smitli [o5] show that significant en- 
c r g y (T El n be s 1 o .t' e d 3 rj, the s u :r f e c e iva %• e s o n t } 1 c f i n 1 d 







near reGonance, v/ilh the result that energy is ti’ad- 
etl back and foi'th betv.een liquid and rigid body. This 
can result it the dcuiiper n^ass is as little as 2-/f ol 
the main body, resulting in a history oi B as sho>vn. 
in Fig.TJI-S [sl] . This can be ovcrcoine by da-mping 
tJie ’ivave 1110 1 ion by the use of baffles., filling the 
void with a light liquid such as alcohol, or usln.g o- 
nough damping fluid so that the void is small and the 
v/aves impact the inner surface of the damper. Also, 
damper masses are usually much small ei" than 2^^ of the 
main body pveight . 

The advantcige of this conf iguration is that it 
assures symmetry in the steady state, v/ith no appar-^ 
ent residua,! v/obble, as is the case with single, and 
some multiple, pendulums. A comparison of a fluid 
damper euid single spin ao^is pivoted pendulum damper 
of equal mass from experimental results is shov/n in 
Fig.III~9l'C)] . 

The IJEOS used a spin axis concentric mercury 
and alcohol damper for small 0 , less than half a ' 
degree, and one equatorial ball damper for fast damp- 
ing at larger 6 . 

Unsymmetr ioally Mounted 

Ayache and Lynch anali^cd toroidal and I'ect angu- 
lar dampers of circular cross section and a U-shaped 
resoria.nt damper mounted in planes parallel to the spin 
axis [ 2 I in terms of a frictional coupling factor f 

JJS 

inversely proportional to the time constant. Only the 
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V e suit r. a.r o pr e s e n t. o d Ii c r c , This i s 1 o r s iii a. 1 X o ii 1 y , 
t }i c s p a c c c r a 1 1 ii c 1 y d e s p u r» * '')i lo g t o ;C the s t al > i. .U. a— 

tioii is due to a llyivlieel cn the spin a>:is^ 

For a to3;oidaX da;aper as described In Fig . Ill'-lOa: 
where il - 2Jj ) 


^0 ) 




jQjJj^ “ Bessel functions of oracr zero 

and one, respectively o • 

V? 

This is plotted vs. e.^{Sl/y) ' for v^iX’iou.s bubble size; 
in Pigo 1II*-101), where a^. is the tube inside radius. 


The rectangular daviiper lias a frictioiial cotipling 
o 

factor (i -- \v ) that of the toroidal, where ( Pig . Ill-lOo ) 

\v' ( a b ) / ( i.1 ^• b ) 

This means a greater time constant. 


For t he U- s h ap e d d. amp or ( F i g , 1 1 1 1 1) 


^Dg ^ iiTiag 


2.a-^ (JVVv./A^)(i --11)1 


vhcro--K:.- - 2(-A7 /'A.J (r/Xr)?^-./^ 
• t. s 

A,j - (v^ + 5 )/Jq(/3^j) 
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Fig. Toroiclal and rectangular liquid dampers 

mounted along transverse a.xis. 





velocity of the tii.bc wall 

c„ ~ ^ ' ; 

N " ■'( 2c..ya,^^) v.r.-.. clr. 
r.,,. -- ro.dial distance from tube center 

Disk Type [26*] 

In this, a disk is mounted on a ball and socket 
at the center of gravity. Fox’’ best results , the fric 
tion should be small* To my knowledge, this type has 
only been used in a, test model by Perkel* 


When the entire body is spinning smoothly and 
then disturbed, the disk damps dov/n more quickly than 
the main body. For si-n.il 1 friction, the dfimpor plane 
is perpendicular to the precession cone axis. Up to 
a point, gr^e a t or f r i e 1 1 o n c an 3 e s fas t ei' darnT) i ng , Th o 
limit is when stiction occurs, freezin.g the dauiipcr. 


Th e damp 1 ng i s exp o n e n t i al 

e -- 



*^3 n 
e~- “A 




where Aj^ =- polar and transverse moments of 

inertia of the damp or 

“ initial angle betvyecn II and disk axis. 
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If the an " 1 e l> e t v; e e ii t l'.i o d i s K. a 1 1 d h o cl y o-x e s is small, 
r may he approximated by H/C. ,v 

The stiction problem can be overcome by using a 
1 ub r i c a t e d b o ar i n g , T h o v i s c o u a ;C r 1 c t i. o n c o ]i s t an t 
for minimum "T is: 

: JCq = CpO, - CAyCjjA)( A -i. l)r/A* 

V\ ltass-“ Spring Systems 
1* Periiendicul.ar to Spin teis 

V/adleigh, Gallov'ay, 'and Mathur have treated a 
spring™ mass systeirt mounted on and pe.rpeiidlcula.r to 
the spin axis [35] « If K is the spring constant, 
c the da,ni})ing, and cv the nature.! frequency: 

r 0 

p -Q q ^ 0 

q -Xip - 2(c/c^) (c«^^^mh/A)X ~ (Kb/A)x 0 
: x + 2(c/c X “ ( ^')x + hrp + hej - 0 

O 1*1 

~ (Tv/m)'^ and » 2/(K,m) ^“ 

If it is assumed tba.t the sinusoidal character of 
the spinniiig body is not affected: 

p - pQ cxp(-d‘\.,t/2) CO silt 

where i\.j. is the Jlayleigh dissipa.tion function: 

p/cpeJiPCTr -r i ) ( 1 “-I- - 


P,,/2 




i of spring (laiiipers 

and perpendicular to 
cis. 
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T i 1 e 111 ax :i ni i irn f imp 1 :l t u do o f t h c in a s s o s c 1 1 1 it i o n Is: 

F i n a 1 1 y » t .Ii e n i..i t <:i t j, o n a n g 1 c is e xp r c s s e d : 

^ 

Becaiiso this system ivill he sliglitly asymmetric, this 
converges to an apparent wobble angle of q^/(r ) , 

In a laboratory test, witli ^ 3cp.s, c/e^ = 0*5, 
cc^^/r 1»1, mh^/ji = 0, 00135, and an initial spin of 
3 cps, all nutation damped out in C sec. See Fig.JIT”!?. 

. Parallel to Si) in Axis 

Such a damper is inherently unbalanced. The nu- 
tation angle is a decreasing exponent iotl ivith a super- 
i np'i o s e d o o n v e r g o n t o s c i 1 1 a t i o n * ;Vgai n , h o we v c r , t b e 
apparent residual wobble is small [2l] . The damper 
is not oj). the spin axis. 

G. Spherical 


A II end Ilium pivoted in a ball and socket ami immersed 
in a fluid was mounted on the despun portion, of OSO [g] . 
However, it will work for a single body satellite for 
C/A I, 



If p.., ^ (G/A)^7 and c is the. damp in g constant o:C 
the fluid: 


i./ 


U 


For resouanoo: 

(C./A.)f= I,Y4!i(Ti:jj)’'V(ra/'* + 2,4G7rii/s^)';'‘ 

where ^ is the pendulum length, is the criamctral 
moment of inortifi oi the pendulum wii^e, and s is the 
transverse radius ox gyration oi the bob* 


For OSO, there Vv^as no evidence of nutation for 
8000 orbital passes. 


Mas s- 1)1" uiii Sy s t cm 


T h i s 5, s an o t. h e r 
and cojif:;ists of tTvo m 
v/rapped around a dxum 
main body concentric 
Bpr ing- damper system . 


system devised by Fvcrhel [26] , 
asses strung on wires vdiich are 
, The drum is connected to the 
with the spin axis by a tortional 
V/hen nutatioii occurs, there is 


a restoring torque due to the relative deflection of 
tlie 'wires in addition to energy dissipation in the 


dampers / Fig . 1 1 l-l 3 , 


Exper ivncntca,! v/ork cn a lab model Indicated this 
system was capable of damping tiio nuta.tion of a pi'o- 
late bodi’’. Of coujrsc, if the cables wore long enough, 
the actual polar inertia niori?ent could be greater than 
that of the pi'olate main body alone, possibly greater 
than the transverse moment of inertia. 



37 



Fig. II 1-13: Mass-clruni nutation damper ajid 
si)in rate control. 
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Ad other xM)f.;sih Hi ty along these linos would he 
to difjpevjse \vith the drum, mounting masses on damped 
springs o n 1 1 1 e o u t s i d o o f th e s p ac e c j:‘ a x t , o jrp o s i t e 
each otlier. In this case there w'ould he no direct 
coupling cl the motion of t)ie tv .'0 dampers* 

I , Mag n 0 1 i c Dam p i n g 

One method is used to align' the spin axis of a 
spacecraft aJong the local external magnetic field* 

A strong, permanent magnet is mounted in the spa.ee- 
craft along the axis* This method vras used in 
TPtANSIT IB c?viid 2Ae The spin had to he reduced to 
h e 1 o \Y 0 * i r s , 0 1 h er w 1 s e , t h e oh 3. at c s p ac c or af t 

could have overcome the magnetic torque arid assumed 
an attitude fixed in space [l3] * 


Energy dissipation also comes about through ed- 
dy currents and magnetic hysteresis* If a rod is 
rotating about a trams verse a.xis perpendicular to the 
external field, the cormponent of the field along the 
rod is a. function of tirae, and thus there must he an 
induced cuin'cnt. This eddy exarrent causes heat to 
he radiated due to the resistance of the members* 

For a spaoecraft of i^clar monicnt C, n number of perm- 
eaJjle rods of volume V and dicniicter D, and spinniiig 
perp Oil (1 i c u j- cxr ■ t o ■ ~t'h e ■ ■ j. 0 0 a 1 -f i c X cl — i nit. i al 1: y '■ a, t ■ 

6?; exp t] 

■ x lO'^'* erg-sco 



wliex'e o-‘ “ separation eifcct ilue to dis.tcinec bet’.vcen 

c 

rods {cr'^ - 1 tor c^) 
o 

- r e s i s 1 1 V i t y o f r ad ( o 1 im~ c m ) 

/ ? s 

(n ■ j ~ avero.'>;e of sauare of maxi mum flux 

density over oiitlre length of rod for 

2 

one orhlt (gauss )* 

. Hystei'csis damping j.s due to the friction hctv/eeii 
the magnetic domains in the spacecraft This results 
in a linea.r damping. 

Note that in all of the above, there t?.re external 
torques, and angular momentum is not conserved, .Since 
there ar*e ejiex’gy losses, ho v: eves’, they can be applied, 
to nutation damping. The latter 'two methods will gen>- 
erally cax.iso energy loss no matter v/hat the oriorita*- 
tion of tlie satellite is intended to be, fixed in space. 

In general, the mag)ictic torques are disturbances 
that must be overcorae by other nutation dam-pers, and 
til us are benef'Joicil only for sx>in removal and ctiign- 
m e n t with the 1 o c al magne tic field. 

Gravity Gradient 


As in t h e fib o v e , t h i s 
damping only v;hcn the spin 
i s { al \:i\y s a }:> r o I at e b o fly 


can. be used for nutation 
is very lov/, and the spin 
in this ease) orienteu 


t o Y,*a.i’ d e ar t h , 1.^ o r t h is type 

i s u s u a 1 1. y d. e s i r e d al o a g t h :i 


of s p a c 0 c i ' a f t , no s ]) i. n 
s axis. In satellites 



uot meant to be gravity gradient stabilized, it is a 
cl i s 1 1 1 r b an c e to be o v e r o o m c b y t c nut a t i o )i d am};) er - 

[20,32] . 

, According to' Thomson, the torque on a satellite 
with spi)i axis' perxiGiidicular to the orbital x>lane is: 

L = ~ C)&^ 

where B - deviation of spin axis from normal to or- 
e 

bital plane tov/ards earth (small) 

-- orbital angular velocity 

and L is about the axis tangent to the oi'bit. Condi” 
tions for stability are defined in terms of: 

ti/2<a;.f y2[-5A - (1 - A )^ ] + A - 1 - 1) A. 

-'/2(^q/«iJ^(A 1)2 

= 4A^ + 5(fj/«'*)A(A 1 )■^(f j/a'., )^'( A -i- 1)^ 

" Spin relOvtive to the tangent to the orbit 
For stability: b^/ 2 cq.,,^ ^ 0 

/oo.f > 0 . 
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K. Sti nclura.i Di.ss:i, pat 


ion 


No Structure is perfectly rigid, and the accelera- 
tions on a preoessing spacecraft will cause energy loss 
througli mechanical hysteresis. Usually, however . part 
p„i t., 01 tlie spaeeorait can bo consiciered rigid with 
energy disetnatlon only from the relatively flexible 
parts, suoli as anteiu-iie or solar panels. Two examples 
have beo]i v’-orked by Thomson [31 ^ 33], 


We have already shown that, for no external torques 

foi on i S3''n!;iie tr i c both' T'bp' t 

uuu., . ine cneroy loss per cycle of 

stress per unit volume is: 

^^cr-VsE 

where W is Young -s modulus, o- the normal stress, and 
the hysteresis factor. Integrating this over the vdiole 
structure, for period of stress oscillation t(j: 




r<3/2EtQ)dV = T 


■ Consider i.is an arbitrary point on the spacecraft 
at coordinates (x.x), we can compute the acceleration 
at thao point,, wliloh is the excitation. If e is com- 
|j 3.1“ 3, 1 1 v e 1 y E m 1 1, ; 

CQj — ( P f., i ]r| Q g j |-J ^ I { S '\ Ji O 0 0 8^)"^ 

-l i^Y' -I- )'ic 


W -- ^f.s-iad? (cosp t - sin^,j') 
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Note that (x, 7^) does define an arbitrary poir3t 
j, n t i 1 0 s p cl c e 0 r ait, not res t r 1 c t e d to o n e p X an c , b o - 
. c a XI s e of t h. e a.x 1 s 3^ 1:. tm c t r 3’' , If the r e 3 . a 1 1 v e mo 1 1 o 1 1 0 f 

the ):)oi]Xts on the spacecraft can be considered small, 
the accelex’at ion at a point is: 

a “ 03 ^ X ( cS ) •{- X-^' 

^ ^ /I 

v/here X- - xi + . 


For an example, let us assume that the elastic 
pai't of a spaxccraft serves onlj^ the fxmetion of en- 
ergy dissipation, and the deflections caxxs.e xxo changes 
in the inertia* The satellite in. Fig* 111-14 consists 
of two disks, each of inertia, and j> and mass , 
connected by a flexible tube of radius and length 
2 -d * The gyroscopic inomexit reqx.)ired b^*- each dis.k is: 

Lg = C^{(^ -r (/' COS& ) ^/ siTUi? “ sin<^^ eos^ 

and C = 

A £? 2(.Aj I- Mp^') 

Then L = siiiC* oosl? 

The moment distribution is linear: 


L 


z 





and thus the maximum stress is: 
er = I.^xj/I 


xvhere I is the cross-sectional moment 
the tu 3 ie.' Substi txit ing these into the 
tion equation: 


of inertia of 
e n ex.'' g 3^ d i s a i p a-' 


• ^ ~ (//24T7- li;)(Uj-<?^'3;..yi)^'(V/C)(G7A)‘6yQ^ si.n£> cos^& 

~ K 00 



Cent, force 



IT-14:: Structural eneri^y dissipation example 



1 1 1 - 3. 5 : Var i ci t ion i n rate of t um'o ling. 



Til is is s h 0 1.1 1 ii F i .XXI - 1 5 » V i s t Ii c v o 1 u.m e of t Ii e 

stressed material and ol>^ the initial angular velocity. 

If A C, such as for a missile: 

p /\ A 

a ^ ’x( 0/A) siiiO cos6> sin y Jv = 


If v;e consider the inertia cf the deflected mem- 
bers, resonance is observed. Pig.III-l6 shows a cylin 
dr'ical spa,oeo,ra.r t o,r radius R, v/ith four beams^ of 
length -4^. If the elastic deformations v/(^,t) are as- 
sumed small and in the z direction only: - 


•El 




m ^\v 


a. 


in 


V/ £ 


Vv'here m and I ax'c the mass and cross-sectional mor:icnts 
of inertia for the beams. Tliis gives: 


0 -- TC... sin^ cos^<9 / [( ^- -y" cos^)^ 
K.,._ = IG Ciiy(i^-.j/3jn + 


( 7/2 

2 


= (1 " c/A)«.'y/.a 

-H. ^ = first natural frequency of beams 


Also, and 
graplri cally 
ance effects 
shape as the 


/5| ai^e tabulated in [37] , This is shovm 
in Fig.IXI-lT, which clearly shows reson- 
. The envelope of this curve is the same 
c ur v 0 i n F i g . 1 1 1 - 1 5 , 




Fiir. TTI-IT: 


Res 0 n a n. c e e ff e c t g = 



SornipasBive and Active Svf;tGKis 




Such fiy stems vrill he me nt toned here only in p£iss- 
ing* A.ctivc systems have energy sources activated 
either on -board sensing equipment or ground coirmjaiKh 

Semi p ji s s i v e h a v e e n e x' gy so ur c e s that e i t Ii or r e m a 1 n 
const an t o r. r e a c t na t u r ally t o a. 1 1 i t u d e c h an g e s » 

0 sc 11 lasting }lass 

This was proposed by Kano and Sobala [l9] . Two 
ma, s s 0 s , d i c t r i c a 1 1 y o p p o s e d , ar o force d to oscillate 
back and forth along the spin axis at constcint frequen- 
cy. The spin axis (axis of syjmnetry) is noriiial to the 
orbital plaue. Tliis is capa.ble of maintaining attitude 
at very low spin r'ates. 

Dual Spin 


The general reasoning behind dual spin spacecra,ft 
with despun datiipax'’s v/as mo.st recently out lijied.hy 
Tojxkin rsij . Fig .lV-1 shows the position of &) rola- 
tive to II for oblate and ju'clate bodies. Vre have 


li 


and 


1 a>^ as the components of la parallel and normal to H, 
respectively.. ...Note.. J:hat:4:hG- in.-each-c-as-e- is-opuo-- 
site In sign. An liiternal torque to reduce nutation 
must- bo of 7.cro average value- (II concervod) and remove 
energy for oblate bodies or inject 'it for prolate. The 
first requirement means that the required torque is 
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no r nui 1 t o 1 T and s p i j'j n i n g vr :i 1 1 ) </> i n i n e r 1 1 a], s pace. 
Since povi-er is tlie scalar pu'odiict of torque and an- 
gular velocity, the torque roqui]'’cd nrast be cppositc 

(c> 5;'oi" oblate bodies a.n.d of the same sense as ec;„for 
c e 

prolate bodies. 

The torques produced by a damper dissipate ener- 
gy, thus the coi'in)OD.Gnt normal to IT is opposite . 

• There is also a, component along the spin axis, thus 
changing the spin of the body upon ^7hioh the damp oj; 
is mounted. If the clamper is despun, the motoi^ must 
comnensato for this speed differential. This is the 
source of energ}^ injection for the prolcite body. Sev- 
er cal refcx*enccs a.i-’e presented in the bibliography. 

C. Ma.gnetic 

i\s was shoxvh before, eddy currents iiiduced by thci 
cartlids magiietic field can cause toi^cjvios on a space- 
craft. This can be overcome by supplying a torqiiing 
coil whose axis is normal to the spin axis v.-ith a cur- 
rent. 180 degrees out of ’phase with the externally in- 
duced 1 j * spin axis can be oriented by an- 

otber coil v»'hose axis j.s parallel to the spin taxis, 

TJ)e current in this is switched on and off, and the 
tor-que being a sinusoid v/hilc on, to give Kero average 
torque on one transverse axis, and a resulta.nt torque 
on the otber. 
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D. Jcit Pulse 


Another iiiethort of supplying torque is to activate 
a single attitude motor aligned paxaliel vrit)i the spin 
axis. The pulsing is controlled by an on-bocird nuta.'- 
tion sensor, firing v^hen the motor is inside the body 
cone [15] . 
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^.3 EFFECTS OF A TOROIDAL 
LIQUID NUTATION DAMPER MOUNTED ON A TRANSVERSE 
AXIS OF AN AXISYMMETRIC SINGLE-SPIN SATELLITE 
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AIvSTUACT 


I n 1 1 1 i s paper, an at. t e m p t J, s m ctrl o (. o (J i s < m > \- o r- U t f' 
parameters relevant to the pcrtormaneo uT a torn i da I 
liquid nutation damper mounted wiUi its axis aloii'..': a, 
transverse axis of a singie-spin saielJ.ite* ama 1 ! iH'- 
1, tial values of t lie nuttit ion ain^lc; wera" a-'SUp;>* 

[5y rlescrlbinti dissipation of enei'i'.y hv the 1 , 1 , ui, d. , a 
time constant for the nutation .aiirtle is found os a ide. 
tiou of a Bessel funetion oT eoriiplox ariiUi'ient , 



NOMENCLATURE 


Small radius of torus 

Radial coordinate within torus 
Function of only 
Constant 

Magnitude of angular momentum of satellite 
Unit vectors along transverse axes of satellite 
Unit vector along spin etxis 

Moments of inertia with respect to transverse and 
spin ax©8, respectively 

Bessel function of order zero 
Kinetic energy of satellite 

XL 

Argument of Bessel function = 

Large radius of torus 

Surface area of control volume 

Time 

Function of time only 
Fluid velocity 

Fluid velocity at wall ( = a) 

Complex fluid velocity for two dampers 


Volume 



w 

r 

. V 

r • 
{') 

( )‘ 

< >1 
0 


Work done on fluid 
Angular velocity 
Precession angle 
Spin angle 
Fluid density 
Absolute viscosity 
Kinematic viscosity 
Shear stress in fluid 
Time derivative except for T 
Derivative for A and T 
Magnitude of complex ( ) 
Angle of nutation 


1 


The minimum energy condition for a spinning 
axisymmetrlc body is when the angular velocity is a- 
ligned'^lth the axis of maximum inertia. When there 
are no external torques, the magnitude of the angular 
momentum, 


H = +«y,/) + 


2 2 
tiU 

z z 


is constant. Values for w and u> exist when the nu- 

* y 

tat ion angle Is not zero. Here we have assumed the 

momenta® of the damper to be small. We also have 

2K = I . (cy ^ + I0J>^ 

t X ^y ' z z 

as twice the kinetic energy, again assuming the motion 
of the damper small. If the nutation angle decreases 
slowly, we have the angular velocity in the satellite 
frame given by 

sine sin^) + i^ijfsinO cos^) 

+ ^z^ f ^ f cos e ) , 

and the precession speed , 

f “ = n/I| . 

Combining these equations gives 

- 2K1^ = - l^)cos‘’e . 

If there is energy dissipation, as with a damper, M 
remains constant while K decreases. Thus (5), 

K = (H^/ljI )(Ij, - cose)©. 



2 . 


With a liquid damper, ener(i;y dissipation oecurs 
because of viscous erfocts, hein : represented hy the 
time rate of work done on tlie rinid by the wall of its 
container ( 2 ) , 

W = V dS 



- 


The damper is illustrated in l^^ig.!. From the 

above, it is seen that the velocity distr ibut ion of the 

liquid must be found. Assuming that. r , , the veloc- 

1 / 

ity of the tube wall relcative to tiie liquid is i^lvcri t)v 
Va = sin& sin^ . 

It is assumed that the velocity will be entirely tan:r;en- 
tial to the torus, and pressure variations due to cen- 
trifugal body forces are small. Thus the fluid momen- 
tum differential vector equation (for a paiticular ^ ) , 
pd\/dl = +y«v^v, 

becomes fdv/^t =yu. 7 "^v, 

or, introducing the kinematic viscosity v ^ ^ 



In ' cylindrical cocrdlnates, this is 
Ov ^ c)v V 

5t = s; 

Assuming that the solution is a product of a function 
A of o^y, and T of t only^ we have 
AT’ = T(a^AM’/a^, 

or a^AT* = T(A» + a^A" ) . 

If transients due to initial conditions are considered 
to damp out -quickly^ the time function will be in phase 
with v^. Thus we let T = sin . Then, 
a^A^cos^ = Vslnp(A’ 4 a^A”) 

Rearranging terms result in 

a^A" 4 A* - (jp/v)cot ^ a^^A = 0, 
which is rather difficult to solve. The coefficient in 
the third term Is a function of time; it cannot be aver- 
aged over one revolution for ^ to give a constant, for 
cot JJ? has values of infinity. 

A simpler model may be had by assuming we have 
another damper, identical to the first, mounted on the 
y axis. Then vie Cfxn represent tlie rotation of ttie 
transverse component of the angular velocity vector us- 
ing complex variables. Since 



4 . 


and e = cosp - j sin^ , 

the velocity component in the xy plane can he described 
by a phasor^ 

”®y = ’ 

^ ,2 Zs'h 

where C3>^y = ^ 

= ^ sin 

Is the phasor magnitude. 

Since we now have the terra e as the excitation, 
we let this be equal to Tg » Then, 

Tg’ = ’ 

and the equation for the fluid velocity position depend- 

function A becomes 

a^A” + A' + (jj;/i^)a^A = 0 , 

O o 2 

or + a^A* + A = 0 . 

This is a complex Bessel equation, the solution of w'hich 


is (1) 

A = CJQ(q) t 
where q a^(j^/^) 

and C is a constant, which may he found by examining the 


velocity at the wall, which is given by 

“.if 


V « = -r ,Cg; e 
a2 t xy 



5 . 


If we define A at = a as where 

then C = Y.(A^ ^J^(q ) 
t S.Y 0 

Thus we have 

Using this expression in the integral for work done 
on the fluid will result in a complex function. This 
would represent two energy flows for dissipation, ninety 
degrees out of phase with each other. Averaged over one 
revolution for jP , the work done could be represented bv 
the magnitude of the complex work. Also at this point we 
can say that this is (tvice the energy dissipation rate 
for a single damper. Thus for the single damper, 


W = 



At any instant, 

\-3f\ 2 2 . 2 ^ 

je =: cos j> H- sin ^ l 
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Making this substitution, ^ve have 


For a given value of 9 , we have 

and ^ ^ cos^ . 

Thus - 1^) 

= jiu ^ V<^t - Iz>> 

Also, again for a particular O with Q small, 


tan © := \m /j o> 

t xy' z z 

Rearranging terms gives 


OD 


xy 




Substituting the relation between cy and a? 

z T » 


C&; 


xy 


= )/(I. - I ) 


we liave 


If we assume no external torques and the Tnomenfiim 
of the damper relative to the satellite main body small, 
the Euler equations for the satellite are 


and 


° = 


0 = I C4^ 
z z 
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Thus we can take o; constant during the damping action, 

z 

and therefore and ^y.y also be held constant. 

Note that this requires that and be small cornpared 

to H , thus meaning 0 is small («^12°); therefore 
z 

sin 9 , 

cos© 5? 1 , 

and tan© =? © . 




8 . 


or © 


3^2 




- 'z 


MVM Jo' ' J 


The term in brackets is almost constant for small nuta- 
tion angles, and is thus the inverse of the time constant 
for a decreasing exponential solution. Thus 
9 = exp(-t/t^) 


where t 


-1 


2 r.Vl., f 


jo 




and 6q is the initial nutation angle. This may also be 
expressed using the approximation for angular momentum by 


-1 


= (2 r//rVit^lJoUa)|'') j 


a^.da^ 




10 . 


CONCUISION 

/ 

The above is valid for small initial nutation an~ 
glea for an axisymmetrlc single-spin satellite. Values 
for complex Bessel functions may be found in references 
3 and 4. However, these are good only for Bessel func- 
tions in which the magnitude of the complex argument q 
is less than ten. However, the best liquid for use in 
the damper is mercury because of its high density; its 
kinematic viscosity is (0.5)10 ft- /sec at 75 F. Since 
q is inversely proportional to the square root of i/ , its 
magnitude will be on the order of 10^ or lO^for reason- 
able values of f . Bessel functions for complex tirgu- 
ments of these magnitudes have not been tabulated, and 


must be calculated. 
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CHAPTER 7 
General Conclusions 

As a result of the present study, equations of motion and 
computer programs have been developed for analyzing the motion of a 
spin-stabilized spacecraft having long, flexible appendages. Sta- 
Ijility charts were derived, or can be redra\i7n with the desired accu- 
racy for any particular set of design parameters. Simulation graphs 
of variables of interest are readily obtainable on line using program 
FLEXAT. Finally, applications to actual satellites, such as UK-4 and 
IMP-I have been considered. 


) 



